Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,190 T→A 5.4% L18I (TTA→ATA)  narJ → molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,1900TA 5.4% 761.3 / 22.9 219L18I (TTA→ATA) narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (101/106);  new base A (5/7);  total (106/113)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 221 total aligned reads displayed.

TGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGA  >  NC_000913/1285043‑1285338
                                                                                                                                                   |                                                                                                                                                     
   tGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                                                                                                                   <  4:1667409/149‑1 (MQ=255)
   tGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATaag                                                                                                                                                   >  1:517590/1‑147 (MQ=255)
     aTGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgc                                                                                                                                                 >  1:299891/1‑149 (MQ=255)
        gCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgcagc                                                                                                                                              <  3:258265/149‑1 (MQ=255)
              aCGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaac                                                                                                                                        >  3:362566/1‑149 (MQ=255)
              aCGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaac                                                                                                                                        >  4:1347856/1‑149 (MQ=255)
              aCGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTCATGGCAGCATcaac                                                                                                                                        <  1:206022/149‑1 (MQ=255)
              aCGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                            >  1:1444094/1‑145 (MQ=255)
              aCGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                            <  2:1444094/145‑1 (MQ=255)
                ggTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                      <  2:161344/149‑1 (MQ=255)
                ggTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                      <  4:186827/149‑1 (MQ=255)
                 gTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGCTGTGACTAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGTTTGTATCGCGTCTCTTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAag                                                                                                                                     >  2:1011754/1‑149 (MQ=255)
                 gTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAag                                                                                                                                     >  2:743474/1‑149 (MQ=255)
                  ttCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                                                                                                                    >  4:885679/1‑149 (MQ=255)
                  ttCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAAAGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                                                                                                                    >  3:784259/1‑149 (MQ=255)
                      gATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCTTGATCGGACTCGTGATTGTATCGCGTCTGCTTGAATTTCCCGATGCTGCCTTATGGCAGCATCAACAAGAGATg                                                                                                                                >  4:866777/1‑149 (MQ=255)
                      gATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATg                                                                                                                                >  3:1683264/1‑149 (MQ=255)
                        tACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAACCGGAGCCGCATCCATGATCGAACTCGTGATTGTCTCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGtt                                                                                                                              <  2:91732/149‑1 (MQ=255)
                         aCCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGttt                                                                                                                             <  4:1623770/149‑1 (MQ=255)
                          ccAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTg                                                                                                                            <  1:743474/149‑1 (MQ=255)
                           cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGTCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGa                                                                                                                           <  3:1650659/149‑1 (MQ=255)
                           cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGtttaa                                                                                                                           >  3:638504/1‑147 (MQ=255)
                           cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGa                                                                                                                           >  4:1220566/1‑149 (MQ=255)
                           cAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGa                                                                                                                           >  3:975087/1‑149 (MQ=255)
                                    tCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGc                                                                                                                  >  3:1172673/1‑149 (MQ=255)
                                          cAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGATCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTCCGTTATGGCAGCATGAAGAAGAGATGTTTGAGGCGATTGCCGCGTc                                                                                                            <  4:1036765/149‑1 (MQ=255)
                                          cAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTc                                                                                                            <  4:1469669/149‑1 (MQ=255)
                                             cAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGaa                                                                                                         <  3:1405467/149‑1 (MQ=255)
                                              aGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGaaa                                                                                                        >  2:771735/1‑149 (MQ=255)
                                                 cgtcgtATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAAt                                                                                                     <  2:675172/149‑1 (MQ=255)
                                                      tATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGcc                                                                                                >  4:106582/1‑149 (MQ=255)
                                                       aTCGATGCCATCGATGTGACCAGCAAAACGGAGCCTCATCCATGATCGAACTCGTGTTTGTATCGCGTCTCCTTGAAGATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGcct                                                                                               >  4:906027/1‑148 (MQ=255)
                                                       aTCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCa                                                                                               <  4:667748/149‑1 (MQ=255)
                                                        tCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                                    >  2:311510/1‑95 (MQ=255)
                                                          gATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCaaaa                                                                                            <  4:1659317/149‑1 (MQ=255)
                                                            tGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                          <  1:321278/149‑1 (MQ=255)
                                                             gCCATCTATGTAACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAAGAGCCGGTTGCTGCCTCAGGGCAGCATCAACAAGAGATGTTTGAGGCGTTTGCCTCATCGAAAAATCTGCCATAagag                                                                                         <  4:784259/149‑1 (MQ=255)
                                                             gCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAagag                                                                                         <  1:661775/149‑1 (MQ=255)
                                                              ccATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAgg                                                                                        >  4:1464052/1‑149 (MQ=255)
                                                              ccATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAgg                                                                                        <  1:1189624/149‑1 (MQ=255)
                                                                aTCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGttt                                                                                                                             <  1:1000682/110‑1 (MQ=255)
                                                                aTCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGttt                                                                                                                             >  2:1000682/1‑110 (MQ=255)
                                                                aTCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGAt                                                                                      >  3:292664/1‑149 (MQ=255)
                                                                 tCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGAt                                                                                      <  4:416381/148‑1 (MQ=255)
                                                                      gtgACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAt                                                                                <  1:926715/149‑1 (MQ=255)
                                                                        gACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgc                                                                              >  1:1097845/1‑149 (MQ=255)
                                                                             gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                        <  2:603821/149‑1 (MQ=255)
                                                                                       gCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCt                                                              <  1:83856/149‑1 (MQ=255)
                                                                                       gCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCt                                                              <  3:1574286/149‑1 (MQ=255)
                                                                                       gCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCAtcctctt                                                              <  3:885679/149‑7 (MQ=255)
                                                                                         cGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTgc                                                            <  2:941978/149‑1 (MQ=255)
                                                                                         cGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTgc                                                            <  4:117886/149‑1 (MQ=255)
                                                                                          gCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGaaaaa                                                                                                      <  1:1491130/107‑1 (MQ=255)
                                                                                          gCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTgcg                                                           >  3:1224574/1‑149 (MQ=255)
                                                                                          gCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTgcg                                                           >  3:1247539/1‑149 (MQ=255)
                                                                                                 tGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCCGCATCAACAAGAGATTTTTGAGACGATTGCCGCGTCGAAAGATCTGCCAAAAGAGGATGCCCATGTGCTGGGCATTTTCCTGCGCGATTTa                                                    >  4:15097/1‑149 (MQ=255)
                                                                                                 tGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                          >  1:989781/1‑112 (MQ=255)
                                                                                                  gATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTTGGCATTTTCCTGCGCGATTTaa                                                   <  1:1040257/149‑1 (MQ=255)
                                                                                                   aTCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                  >  4:1416585/1‑149 (MQ=255)
                                                                                                    tCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGGAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGGCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacg                                                 >  2:1124633/1‑149 (MQ=255)
                                                                                                       aaCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgacg                                              >  4:609930/1‑149 (MQ=255)
                                                                                                         cTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGAt                                            <  3:1113506/149‑1 (MQ=255)
                                                                                                          tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATg                                           <  2:1413702/149‑1 (MQ=255)
                                                                                                          tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATg                                           >  3:1059235/1‑149 (MQ=255)
                                                                                                          tCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATg                                           >  2:1475237/1‑149 (MQ=255)
                                                                                                            gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                        >  2:287769/1‑118 (MQ=255)
                                                                                                            gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGa                                         >  4:113385/1‑149 (MQ=255)
                                                                                                            gTGATTGTATCGCGTCTCCTTGAATATCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                        <  2:879770/149‑1 (MQ=255)
                                                                                                             tGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                        <  3:1347856/149‑1 (MQ=255)
                                                                                                              gATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATc                                       >  1:876652/1‑149 (MQ=255)
                                                                                                              gATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATc                                       <  1:39634/149‑1 (MQ=255)
                                                                                                               aTTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCACGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATcc                                      <  3:1026330/149‑1 (MQ=255)
                                                                                                                     tCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGa                                >  1:325639/1‑149 (MQ=255)
                                                                                                                      cgcgTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGAt                               <  3:609930/149‑1 (MQ=255)
                                                                                                                      cgcgTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGAt                               <  3:525795/149‑1 (MQ=255)
                                                                                                                      cgcgTCTCCTTGAATATCCGGAAGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGAt                               <  2:1186817/149‑1 (MQ=255)
                                                                                                                        cgTCTCGTTGAAGATCCGGATGCTGCGTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCACTGGGCATTTTCCTGCACGATTTAACGACGATGGATCCGCTAGATGc                             <  4:48893/149‑1 (MQ=255)
                                                                                                                         gTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGcc                            >  2:1718541/1‑149 (MQ=255)
                                                                                                                          tctcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                           <  1:771735/149‑1 (MQ=255)
                                                                                                                               ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                           <  3:1192807/144‑1 (MQ=255)
                                                                                                                               ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                      >  2:986671/1‑149 (MQ=255)
                                                                                                                               ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                      >  2:459175/1‑149 (MQ=255)
                                                                                                                               ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                      >  2:542066/1‑149 (MQ=255)
                                                                                                                               ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGgcg                      >  1:573216/1‑149 (MQ=255)
                                                                                                                                 gAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGAt                                                                                      <  2:187603/84‑1 (MQ=255)
                                                                                                                                 gAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGAt                                                                                      >  1:187603/1‑84 (MQ=255)
                                                                                                                                  aaTATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGccc                                                                                  <  2:1630133/87‑1 (MQ=255)
                                                                                                                                  aaTATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGccc                                                                                  >  1:1630133/1‑87 (MQ=255)
                                                                                                                                     atCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                          <  4:222144/76‑1 (MQ=255)
                                                                                                                                       ccGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCAGCGCATGGGCa                                                                     >  4:1435971/1‑94 (MQ=255)
                                                                                                                                        ttgATGCTGCCATAAGGCAGCATCCACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGa                                                                                       <  3:1435537/74‑1 (MQ=38)
                                                                                                                                         ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCg                                     >  4:207883/1‑124 (MQ=255)
                                                                                                                                         ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCg                                     <  3:207883/124‑1 (MQ=255)
                                                                                                                                         ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg            >  1:434175/1‑149 (MQ=255)
                                                                                                                                             gctgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACt        >  3:1440738/1‑149 (MQ=255)
                                                                                                                                              ctgcCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTg       >  3:722666/1‑149 (MQ=255)
                                                                                                                                                gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAac                                                  <  4:356337/104‑1 (MQ=255)
                                                                                                                                                gcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGtt     <  2:391635/149‑1 (MQ=255)
                                                                                                                                                 ccTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAag                                                                                           <  2:129533/63‑1 (MQ=255)
                                                                                                                                                 ccTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAag                                                                                           >  1:129533/1‑63 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGA  >  NC_000913/1285043‑1285338

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: