Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,710,708 T→C 5.4% intergenic (‑233/+91) eptB ← / ← yhjX Kdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,7080TC 5.4% 227.2 / 2.7 75intergenic (‑233/+91)eptB/yhjXKdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter
Reads supporting (aligned to +/- strand):  ref base T (41/29);  new base C (2/2);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCGTCGTA  >  NC_000913/3710560‑3710856
                                                                                                                                                    |                                                                                                                                                    
taTAAGTAAAATCTACACGCTCGAATAAGTAGACCACAGCAAAGCGAAAAAGTGTGGGGAAAGAGTTTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                                                                                     <  1:1309479/149‑1 (MQ=255)
    aGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGcc                                                                                                                                                 <  1:373539/150‑1 (MQ=255)
     gTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCAAGAAGCACAAAATTACAAAAATGCAAGTATGAAAAGAGACAATAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAATGGcac                                                                                                                                                >  2:1301256/1‑148 (MQ=255)
      tAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCa                                                                                                                                               >  1:595969/1‑150 (MQ=255)
       aaaaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAATCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                                                                              >  2:1032056/1‑150 (MQ=255)
        aaaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaaa                                                                                                                                             >  2:414409/1‑150 (MQ=255)
         aaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttt                                                                                                                                            <  2:667038/150‑5 (MQ=255)
         aaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttt                                                                                                                                            <  2:367437/150‑5 (MQ=255)
             tACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTa                                                                                                                                        <  1:648722/150‑1 (MQ=255)
              acacGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                                                                       <  2:1412967/150‑1 (MQ=255)
              acacGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                                                                       <  1:610254/150‑1 (MQ=255)
                          aaGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGttt                                                                                                                           <  1:1173349/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                          >  2:181649/1‑150 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                          <  1:844671/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                          <  1:1619328/150‑1 (MQ=255)
                              aGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATAc                                                                                                                       <  1:391378/150‑1 (MQ=255)
                                     cgCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAACAGAGACTTTAGCGGTACAAAAAGCAGTAAAAAGGGCCAGGTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                                                                                <  2:307579/149‑1 (MQ=255)
                                           gggAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGc                                                                                                          <  2:656068/150‑1 (MQ=255)
                                               aaaaGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                                                                      <  1:181649/150‑1 (MQ=255)
                                                aaaGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTg                                                                                                     <  1:414409/150‑1 (MQ=255)
                                                     gtgGGGAAAGAGTGTTCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                                                                                <  1:1032056/150‑1 (MQ=255)
                                                     gtgGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                                                                                <  1:955202/150‑1 (MQ=255)
                                                      tgGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGtttcgccagc                                                                                               >  2:877504/1‑143 (MQ=255)
                                                      tgGGGAAAGAGTGTGAATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGAATTTAGCGTTACAAAAAGCATTAAAAAGGGCAAGTTAATTGTCCAAAATTTTACGAGTCATTGTTTATACATTTCGTTTATGAGtttcgacagc                                                                                               >  2:871896/1‑143 (MQ=255)
                                                           aaaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                                                                          <  2:46747/150‑1 (MQ=255)
                                                               agTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCt                                                                                      >  2:686248/1‑150 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                                                                                 >  2:995480/1‑150 (MQ=255)
                                                                         aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGt                                                                            <  1:1250542/150‑1 (MQ=255)
                                                                           gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAg                                                                          >  1:13381/1‑150 (MQ=255)
                                                                           gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAg                                                                          <  1:1432600/150‑1 (MQ=255)
                                                                              caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                                                                       >  1:84942/1‑150 (MQ=255)
                                                                                  tttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCgg                                                                   >  2:983467/1‑150 (MQ=255)
                                                                                   ttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc                                                                  <  1:941972/150‑1 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                               >  2:608373/1‑150 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                               >  2:1241266/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                       >  2:750133/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                       >  1:223967/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                       >  2:1323008/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGACAGTTAATTGGCCCAAATTTTACGATTCATTGTTTATACATATCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCAAGTAGCTCACGGCTGAGCATTaaa                                                       >  1:181970/1‑150 (MQ=255)
                                                                                                 gTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAggg                                                    <  1:1586922/150‑1 (MQ=255)
                                                                                                    tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGc                                                 >  1:1196383/1‑150 (MQ=255)
                                                                                                    tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGc                                                 >  2:46819/1‑150 (MQ=255)
                                                                                                      aaaaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCa                                               >  2:634228/1‑150 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                           >  2:172068/1‑150 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                           <  1:849282/150‑1 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                           <  1:184488/150‑1 (MQ=255)
                                                                                                             aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcca                                        >  1:1148565/1‑149 (MQ=255)
                                                                                                               tttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCa                                      <  2:1432459/150‑1 (MQ=255)
                                                                                                                 tAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACg                                    >  1:157752/1‑150 (MQ=255)
                                                                                                                      ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACa                               >  1:709288/1‑150 (MQ=255)
                                                                                                                       tACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAt                              >  1:549413/1‑150 (MQ=255)
                                                                                                                        aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAtt                             >  1:597429/1‑150 (MQ=255)
                                                                                                                         cAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAttt                            >  1:149364/1‑150 (MQ=255)
                                                                                                                            aaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCt                         <  2:1196383/150‑1 (MQ=255)
                                                                                                                              aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGc                       >  2:906852/1‑150 (MQ=255)
                                                                                                                               gCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGct                      >  2:935411/1‑150 (MQ=255)
                                                                                                                               gCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGct                      >  1:1526199/1‑150 (MQ=255)
                                                                                                                                  gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTg                   >  2:1481983/1‑150 (MQ=255)
                                                                                                                                   aaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTgg                  >  1:665676/7‑150 (MQ=255)
                                                                                                                                    aatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc                 >  2:2221/6‑150 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcc                                         <  1:1082464/126‑1 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcc                                         >  2:1082464/1‑126 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATTAGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt                >  2:219679/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGTAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt                <  1:952782/150‑1 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt                <  1:686248/150‑1 (MQ=255)
                                                                                                                                        aGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGAc             <  2:149364/150‑1 (MQ=255)
                                                                                                                                          ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGa           >  1:1329206/1‑150 (MQ=255)
                                                                                                                                          ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGa           >  1:1093880/1‑150 (MQ=255)
                                                                                                                                            ccAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAt         >  1:1524407/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg       >  2:84887/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg       >  1:253104/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg       >  1:1252202/1‑150 (MQ=255)
                                                                                                                                                ttAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcgtc     <  2:709288/150‑1 (MQ=255)
                                                                                                                                                 tAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcgtcg    <  2:1093880/150‑1 (MQ=255)
                                                                                                                                                  aaTTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGAACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATcgtcgt   >  2:642016/1‑150 (MQ=255)
                                                                                                                                                   aTTGGCCCAAATTTTTCCAGTCATTGTTTATACATTTCGTTGATCCATTTGATTCACACGTAAGCTGAAGCCAAGTAGCTCGCGGCTGACCATTAAAGGGAGCCATCCGCCTCACGTAACATTTTCTGCTATGGCTGACGAATCGTCGTa  >  2:1221250/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCGTCGTA  >  NC_000913/3710560‑3710856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑