Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,294,366 A→T intergenic (+227/‑115) nrfG → / → gltP heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,294,3660AT99.0% 376.0 / ‑3.4 97intergenic (+227/‑115)nrfG/gltPheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (42/54);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCCTGGCCTGGC  >  NC_000913/4294218‑4294511
                                                                                                                                                    |                                                                                                                                                 
cAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                                                                                                 >  1:1156206/1‑151 (MQ=255)
    attattGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                                                                                              <  1:943945/150‑1 (MQ=255)
     ttattGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                                                                                             <  2:106137/150‑1 (MQ=255)
          gCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                                                                                                        >  1:180396/1‑150 (MQ=255)
          gCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTa                                                                                                                                       <  1:254736/151‑1 (MQ=255)
          gCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTa                                                                                                                                       <  2:416475/151‑1 (MQ=255)
           cAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                                                                                                      <  2:1088014/151‑1 (MQ=255)
                tAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                                                                                                  >  1:373802/1‑150 (MQ=255)
                tAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                                                                                                  <  2:96206/150‑1 (MQ=255)
                tAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                                                                                                  >  2:783897/1‑150 (MQ=255)
                tAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                                                                                                  >  2:881638/1‑150 (MQ=255)
                   aCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                                                                                              >  2:1155641/1‑151 (MQ=255)
                    cAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                                                                                              >  1:202806/1‑150 (MQ=255)
                      aaTAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                                                                                            >  1:350810/1‑150 (MQ=255)
                         aaaaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                                                                                         <  2:223635/150‑1 (MQ=255)
                         aaaaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                                                                                                        >  1:691109/1‑151 (MQ=255)
                         aaaaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                                                                                                        <  1:999727/151‑1 (MQ=255)
                          aaaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                                                                                                       >  1:302489/1‑151 (MQ=255)
                            aaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCa                                                                                                                      >  2:246090/1‑150 (MQ=255)
                              tCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                                                                                                   >  1:974463/1‑151 (MQ=255)
                                tttCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGc                                                                                                                  >  1:668475/1‑150 (MQ=255)
                                   ccATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                                                                                              <  1:665317/151‑1 (MQ=255)
                                   ccATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAac                                                                                                               <  2:72155/150‑1 (MQ=255)
                                    cATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                             >  2:844946/1‑151 (MQ=255)
                                    cATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                             <  1:776805/151‑1 (MQ=255)
                                     aTAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                                                                                             >  1:639716/1‑150 (MQ=255)
                                      tAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                                                                                            >  1:1141576/1‑150 (MQ=255)
                                       aaCAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                                                                                          <  1:268043/151‑1 (MQ=255)
                                       aaCAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                                                                                          >  2:188653/1‑151 (MQ=255)
                                          aaaTGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                                                                                                       <  1:357683/151‑1 (MQ=255)
                                              ggTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                                                                                                   >  1:14749/1‑151 (MQ=255)
                                              ggTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                                                                                                    <  2:997586/150‑1 (MQ=255)
                                                ttATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                                                                                                   <  2:914225/149‑1 (MQ=255)
                                                ttATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                                                                                                 >  1:929922/1‑151 (MQ=255)
                                                ttATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                                                                                                 >  1:820724/1‑151 (MQ=255)
                                                  aTTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTgg                                                                                               <  1:557118/151‑1 (MQ=255)
                                                      aTTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                                                                                           >  1:1000335/1‑151 (MQ=255)
                                                        tAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                                                                                         >  2:109073/1‑151 (MQ=255)
                                                         aaTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                                                                                         <  2:1034492/150‑1 (MQ=255)
                                                         aaTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                                                                                         <  2:598646/150‑1 (MQ=255)
                                                         aaTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg                                                                                        >  2:153828/1‑151 (MQ=255)
                                                            ccTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTc                                                                                      <  1:691087/150‑1 (MQ=255)
                                                            ccTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCa                                                                                     <  2:925929/151‑1 (MQ=255)
                                                                  aTAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct                                                                               <  2:333034/151‑1 (MQ=255)
                                                                  aTAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct                                                                               >  2:1141719/1‑151 (MQ=255)
                                                                  aTAAACTCGCGTTATTAAGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct                                                                               <  1:1217654/151‑1 (MQ=255)
                                                                    aaaCTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg                                                                             <  2:178961/151‑1 (MQ=255)
                                                                    aaaCTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg                                                                             >  1:720076/1‑151 (MQ=255)
                                                                    aaaCTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg                                                                             <  1:241919/151‑1 (MQ=255)
                                                                            gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                                     >  2:443332/1‑151 (MQ=255)
                                                                               attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAAc                                                                  >  1:1052633/1‑151 (MQ=255)
                                                                                  atGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACggg                                                               >  1:1214550/1‑151 (MQ=255)
                                                                                   tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                              >  1:731548/1‑151 (MQ=255)
                                                                                    gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                              <  1:757712/150‑1 (MQ=255)
                                                                                     cATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGaa                                                            <  2:256550/151‑1 (MQ=255)
                                                                                         aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACgg                                                         <  1:273846/150‑1 (MQ=255)
                                                                                          aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgcg                                                       >  2:757938/1‑151 (MQ=255)
                                                                                            gCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCt                                                     <  2:570902/151‑1 (MQ=255)
                                                                                             cAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                                    >  2:1071513/1‑151 (MQ=255)
                                                                                                cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAt                                                 >  2:641407/1‑151 (MQ=255)
                                                                                                 gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                                <  2:574515/151‑1 (MQ=255)
                                                                                                 gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                                <  2:184460/151‑1 (MQ=255)
                                                                                                  aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTg                                               >  1:254175/1‑151 (MQ=255)
                                                                                                    aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAg                                             <  1:652283/151‑1 (MQ=255)
                                                                                                     aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                            >  1:318356/1‑151 (MQ=255)
                                                                                                     aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                            >  1:12159/1‑151 (MQ=255)
                                                                                                           gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGt                                        <  2:1204613/149‑1 (MQ=255)
                                                                                                             tgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcat                                     <  1:384149/150‑1 (MQ=255)
                                                                                                              gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattc                                   <  1:835961/151‑1 (MQ=255)
                                                                                                              gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattc                                   <  2:422466/151‑1 (MQ=255)
                                                                                                                tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattca                                  <  2:798504/150‑1 (MQ=255)
                                                                                                                tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattca                                  <  1:992511/150‑1 (MQ=255)
                                                                                                                  aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCatat                               >  1:1032947/1‑151 (MQ=255)
                                                                                                                  aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCatat                               <  1:963970/151‑1 (MQ=255)
                                                                                                                   aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATg                              >  2:1114909/1‑151 (MQ=255)
                                                                                                                    gggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGa                             <  1:188658/151‑1 (MQ=255)
                                                                                                                     ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaa                            <  2:1184159/151‑1 (MQ=255)
                                                                                                                        tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaaaa                          >  2:463568/1‑150 (MQ=255)
                                                                                                                        tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaaaa                          <  2:715051/150‑1 (MQ=255)
                                                                                                                        tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGaaaa                          <  1:597657/150‑1 (MQ=255)
                                                                                                                         gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAat                        <  1:164343/151‑1 (MQ=255)
                                                                                                                         gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAat                        <  2:1156175/151‑1 (MQ=255)
                                                                                                                          cgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAata                       <  1:671673/151‑1 (MQ=255)
                                                                                                                           gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAatat                      <  2:837405/151‑1 (MQ=255)
                                                                                                                             aaCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAAtata                     >  2:756613/1‑150 (MQ=255)
                                                                                                                              aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATaaa                   <  2:738865/151‑1 (MQ=255)
                                                                                                                              aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATaaa                   <  1:148205/151‑1 (MQ=255)
                                                                                                                              aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATaaa                   <  2:202801/151‑1 (MQ=255)
                                                                                                                                aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAAt                 >  1:449310/1‑151 (MQ=255)
                                                                                                                                  aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAAtt                <  2:1166933/150‑1 (MQ=255)
                                                                                                                                   ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCa              <  2:1148459/151‑1 (MQ=255)
                                                                                                                                     gcgcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgc            <  2:1177379/151‑1 (MQ=255)
                                                                                                                                             gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctggcctg    <  1:785592/151‑1 (MQ=255)
                                                                                                                                             gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAgcctggcctg    >  1:273864/1‑151 (MQ=255)
                                                                                                                                              aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGcctggcctgg   <  1:246095/151‑1 (MQ=255)
                                                                                                                                                 cTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCctggcctggc  >  2:740524/1‑149 (MQ=255)
                                                                                                                                                 cTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCctggcctggc  >  2:392526/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATATGAAAAATATAAAATTCAGCCTGGCCTGGC  >  NC_000913/4294218‑4294511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑