Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,470,108 A→C 5.8% intergenic (‑255/+37) chiA ← / ← tufA periplasmic endochitinase/translation elongation factor EF‑Tu 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,470,1080AC 5.8% 1771.8 / 75.4 509intergenic (‑255/+37)chiA/tufAperiplasmic endochitinase/translation elongation factor EF‑Tu 1
Reads supporting (aligned to +/- strand):  ref base A (242/236);  new base C (9/21);  total (252/257)
Fisher's exact test for biased strand distribution p-value = 3.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 519 total aligned reads displayed.

AAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTT  >  NC_000913/3469959‑3470257
                                                                                                                                                     |                                                                                                                                                     
aaGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGa                                                                                                                                                       <  1:5038469/150‑1 (MQ=255)
   aaaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgc                                                                                                                                                    <  2:1161216/150‑1 (MQ=255)
   aaaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgc                                                                                                                                                    >  1:2994097/1‑150 (MQ=255)
    aaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgcg                                                                                                                                                   <  1:8650796/150‑1 (MQ=255)
    aaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgcg                                                                                                                                                   <  1:8175502/150‑1 (MQ=255)
         gCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCtt                                                                                                                                              <  2:7064221/150‑1 (MQ=255)
          ccTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttt                                                                                                                                             <  2:2252838/150‑1 (MQ=255)
           cTGACAAATTAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                            <  2:8413179/150‑1 (MQ=255)
           cTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                            <  2:5264007/150‑1 (MQ=255)
           cTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                            <  2:5759050/150‑1 (MQ=255)
           cTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                            >  2:6132288/1‑150 (MQ=255)
            tGACAACTGAAGCATTTTAAAAACAGAACCATTCATATTTAAAATGTGAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttc                                                                                                                                           <  2:4612152/150‑1 (MQ=255)
            tGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttttt                                                                                                                                           <  2:8995262/150‑2 (MQ=255)
            tGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttc                                                                                                                                           >  2:7955265/1‑150 (MQ=255)
            tGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttc                                                                                                                                           >  1:4597667/1‑150 (MQ=255)
             gACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttca                                                                                                                                          >  2:2762059/1‑150 (MQ=255)
             gACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttca                                                                                                                                          <  1:5419035/150‑1 (MQ=255)
                aaaTGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTTAATAGATATTATAAATATGAATAATTTATTAGTACTGACAGTACGAATAAGATATGCCGACAACAAATGCAAAAATGGCGCCGAAGCGCCCTTttcaatt                                                                                                                                       >  1:6166861/1‑150 (MQ=255)
                   tGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaattcaa                                                                                                                                    <  2:1215167/150‑1 (MQ=255)
                   tGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaattcaa                                                                                                                                    <  2:9552967/150‑1 (MQ=255)
                         aTTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTaa                                                                                                                              <  1:2517204/150‑1 (MQ=255)
                           tttAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAAtt                                                                                                                            <  2:3541202/150‑1 (MQ=255)
                              aaaaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAAc                                                                                                                         <  1:2931154/150‑1 (MQ=255)
                                aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                       >  1:9396225/1‑150 (MQ=255)
                                aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                       >  2:3910192/1‑150 (MQ=255)
                                aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                       <  2:1427638/150‑1 (MQ=255)
                                aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                       >  1:2395560/1‑150 (MQ=255)
                                aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                       >  2:6432866/1‑150 (MQ=255)
                                     gAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAAt                                                                                                                  <  1:736281/150‑1 (MQ=255)
                                           ttCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGccc                                                                                                            <  2:3672420/150‑1 (MQ=255)
                                            tCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCa                                                                                                           <  2:4652282/150‑1 (MQ=255)
                                            tCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCa                                                                                                           <  2:506194/150‑1 (MQ=255)
                                            tCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCa                                                                                                           <  2:1352947/150‑1 (MQ=255)
                                             cATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACCAGTAATTAGCCCAg                                                                                                          >  2:8565645/1‑150 (MQ=255)
                                               tatTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGaa                                                                                                        <  1:5948271/150‑1 (MQ=255)
                                               tatTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGaa                                                                                                        >  1:754396/1‑150 (MQ=255)
                                                atTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAAc                                                                                                       <  1:4852935/150‑1 (MQ=255)
                                                 tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACt                                                                                                      >  2:4279487/1‑150 (MQ=255)
                                                 tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACt                                                                                                      >  1:2237921/1‑150 (MQ=255)
                                                 tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACt                                                                                                      >  2:6589872/1‑150 (MQ=255)
                                                   tAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACttt                                                                                                    >  1:6501497/1‑150 (MQ=255)
                                                    aaaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCTTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCAAGAACTTTa                                                                                                   >  2:3076467/1‑150 (MQ=255)
                                                    aaaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTa                                                                                                   >  2:5614630/1‑150 (MQ=255)
                                                    aaaaTGTTAAATTGAATTGATATATTAAATATGAATAATTTATTCGTACTGACAGAAAGAATAAGATATGCCGAAAACAAATGCAAAAAGTGCGACGAATCGACCTTATAAAATCAAAACTAATTAACGTGTAAATAGCCCATAACTTTa                                                                                                   >  2:1866040/1‑150 (MQ=255)
                                                         gTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaac                                                                                              <  2:1265674/150‑1 (MQ=255)
                                                                  aattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcgc                                                                                     >  2:2985020/1‑150 (MQ=255)
                                                                       atatTTTAAATTTTAATAATTTATTGTTTCTGACAGTCCTAATATTATATGCCGTAAACAAATGAAAAATGGGCGCCGAAGCGCCCTTTTAAATTCAAAACTAATTAAGGTGTAATTAGCCCATAACTTTAGCAACAACGCCCGCGCCaa                                                                                <  1:951481/150‑1 (MQ=255)
                                                                          ttttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGACGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACgg                                                                             >  1:5507698/1‑150 (MQ=255)
                                                                           tttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGt                                                                            >  1:3980173/1‑150 (MQ=255)
                                                                            ttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTa                                                                           >  1:6940853/1‑150 (MQ=255)
                                                                            ttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTa                                                                           <  2:2914126/150‑1 (MQ=255)
                                                                            ttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTa                                                                           <  2:3329115/150‑1 (MQ=255)
                                                                               aaTATTAATAATTTATTCGTTCTGACGGTACTAATAAGATAGGCCGTCAACAAATGAAAAAAGGGCGCCTAAGCGCCATTTTAAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACgg                                                                        <  1:3397582/150‑1 (MQ=38)
                                                                                 tatGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgcc                                                                      >  1:8042368/1‑150 (MQ=255)
                                                                                   tGAATAATTTTTTCGTTCTGACATTACTAATAATATATGCCGTCAACAAATGAAAAAAGGGCGCCGAAGCGCCCTTTTAAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgc                                                                    <  1:8814852/150‑1 (MQ=39)
                                                                                   tGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgc                                                                    <  1:7876542/150‑1 (MQ=255)
                                                                                   tGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgc                                                                    <  2:815815/150‑1 (MQ=255)
                                                                                       taatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTc                                                                >  1:4792515/1‑150 (MQ=255)
                                                                                       taatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTc                                                                >  2:3183999/1‑150 (MQ=255)
                                                                                        aatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCa                                                               <  2:2470499/150‑1 (MQ=255)
                                                                                        aatTTATTAGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCa                                                               <  2:5855846/150‑1 (MQ=255)
                                                                                          tttATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAACGGCAAAAAGGGCGCCGACGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACg                                                             <  1:7076680/150‑1 (MQ=255)
                                                                                            ttttCGTTATGACATTCCGATTAAGATATGCCTTCAACAACTGAAACATGGGCGCCGACGCGCCCTTTTCAATTAAAACATAATTACCTTGTAATTTGCCCTGACCTTTAGCACCAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGa                                                           <  1:963787/148‑1 (MQ=19)
                                                                                            tATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGa                                                           <  1:1293091/150‑1 (MQ=255)
                                                                                              ttCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGAtt                                                         >  2:196732/1‑150 (MQ=255)
                                                                                               tCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTg                                                        <  2:4964151/150‑1 (MQ=255)
                                                                                                  ttCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGACGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGa                                                     <  1:2745171/150‑1 (MQ=255)
                                                                                                  ttCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGa                                                     >  1:6103578/1‑150 (MQ=255)
                                                                                                   tCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAACCTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaa                                                    <  1:6753496/150‑1 (MQ=255)
                                                                                                      gACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACg                                                 >  2:8341722/1‑150 (MQ=255)
                                                                                                       aCAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGc                                                <  2:9283149/150‑1 (MQ=255)
                                                                                                       aCAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGc                                                >  2:4419520/1‑150 (MQ=255)
                                                                                                           tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGAc                                            >  1:8381557/1‑150 (MQ=255)
                                                                                                              gAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgt                                         >  1:1983173/1‑150 (MQ=255)
                                                                                                               aaTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtc                                        <  1:3569104/150‑1 (MQ=255)
                                                                                                                  aaGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCgtcgtc                                     <  1:1824551/150‑1 (MQ=255)
                                                                                                                    gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCa                                   <  2:4046902/150‑1 (MQ=255)
                                                                                                                      tatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                 <  2:4756643/150‑1 (MQ=255)
                                                                                                                       atGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcg                                <  2:7153774/150‑1 (MQ=255)
                                                                                                                       atGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcg                                >  1:5342480/1‑150 (MQ=255)
                                                                                                                        tGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTACTTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcgc                               <  2:3839024/150‑1 (MQ=255)
                                                                                                                        tGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcgc                               >  2:7553779/1‑150 (MQ=255)
                                                                                                                        tGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcgc                               <  2:8907737/150‑1 (MQ=255)
                                                                                                                           cGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGAt                            >  2:9170084/1‑150 (MQ=255)
                                                                                                                           cGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGAt                            >  1:6352524/1‑150 (MQ=255)
                                                                                                                             tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCg                          >  1:4397/1‑150 (MQ=255)
                                                                                                                               aacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCggg                        >  1:1104861/1‑150 (MQ=255)
                                                                                                                               aacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCggg                        <  2:8705047/150‑1 (MQ=255)
                                                                                                                                acaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGt                       <  1:1688960/150‑1 (MQ=255)
                                                                                                                                acaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGt                       >  1:9268521/1‑150 (MQ=255)
                                                                                                                                acaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGt                       >  2:6855239/1‑150 (MQ=255)
                                                                                                                                acaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGt                       <  1:3783389/150‑1 (MQ=255)
                                                                                                                                  aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTgg                     >  2:5534621/1‑150 (MQ=255)
                                                                                                                                       tgaaaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAg                >  2:2717104/3‑150 (MQ=255)
                                                                                                                                          aaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt             >  1:6367543/1‑150 (MQ=255)
                                                                                                                                               gcgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaacaa        >  2:6256913/1‑150 (MQ=255)
                                                                                                                                               gcgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaacaa        >  1:3672564/1‑150 (MQ=255)
                                                                                                                                                 gcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAAcc      >  1:9054285/1‑150 (MQ=255)
                                                                                                                                                 gcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAAcc      >  1:5724061/1‑150 (MQ=255)
                                                                                                                                                     aaGCGCCCTTTTCAATTCAAAACTACTTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCAttt  <  1:1560387/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTT  >  NC_000913/3469959‑3470257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: