Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,686 T→C 5.1% Q169R (CAG→CGG)  flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,6860TC 5.1% 398.7 / 8.6 118Q169R (CAG→CGG) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base T (60/52);  new base C (3/3);  total (63/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 118 total aligned reads displayed.

TTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCA  >  NC_000913/1965539‑1965833
                                                                                                                                                   |                                                                                                                                                    
tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                                                    >  1:1704579/1‑150 (MQ=255)
tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                                                    >  2:787737/1‑150 (MQ=255)
tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                                                    >  2:787690/1‑150 (MQ=255)
  tttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCgg                                                                                                                                                  <  2:510021/150‑1 (MQ=255)
     cAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                                                               >  2:100236/1‑150 (MQ=255)
         tGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGa                                                                                                                                           >  1:1041472/1‑150 (MQ=255)
         tGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGa                                                                                                                                           >  1:1041473/1‑150 (MQ=255)
         tGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGa                                                                                                                                           >  1:1547638/1‑150 (MQ=255)
         tGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGa                                                                                                                                           >  1:12604/1‑150 (MQ=255)
          ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                                                          <  2:1216564/150‑1 (MQ=255)
            cTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGc                                                                                                                                        >  2:460762/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  1:577263/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  1:401419/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  1:541550/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  2:828284/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  1:981156/1‑150 (MQ=255)
             tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                                                       >  1:1044689/1‑150 (MQ=255)
                atatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                                                    <  2:1201757/150‑1 (MQ=255)
                atatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                                                    <  2:1311007/150‑1 (MQ=255)
                     tGGAAAAAGACGACAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑aaaa                                                                                                                              <  1:1564697/150‑1 (MQ=255)
                       gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑aaaaaa                                                                                                                            <  2:600849/150‑1 (MQ=255)
                       gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑aaaaaa                                                                                                                            <  2:779281/150‑1 (MQ=255)
                       gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑aaaaaa                                                                                                                            >  2:1648119/1‑150 (MQ=255)
                           aaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAAcccc                                                                                                                        <  1:1122923/150‑1 (MQ=255)
                               cGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCa                                                                                                                    <  1:1423624/150‑1 (MQ=255)
                                  cAACTCCCACCAGTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 <  2:1352834/150‑1 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 <  2:1078953/150‑1 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 <  1:926914/150‑1 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 <  1:372336/150‑1 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 <  1:372333/150‑1 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 >  1:1131914/1‑150 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 >  2:1513854/1‑150 (MQ=255)
                                  cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGc                                                                                                                 >  1:1305408/1‑150 (MQ=255)
                                     aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGc                                                                                                              >  2:1482262/1‑150 (MQ=255)
                                     aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGc                                                                                                              >  2:36391/1‑150 (MQ=255)
                                     aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGc                                                                                                              <  2:1480624/150‑1 (MQ=255)
                                     aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGc                                                                                                              >  2:1482218/1‑150 (MQ=255)
                                         caccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAAc                                                                                                          >  1:1177783/1‑150 (MQ=255)
                                           ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCa                                                                                                        <  1:1696915/150‑1 (MQ=255)
                                             aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAgg                                                                                                      >  1:80403/1‑150 (MQ=255)
                                                  aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                                                 >  1:1651076/1‑150 (MQ=255)
                                                  aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                                                 >  1:1289743/1‑150 (MQ=255)
                                                  aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                                                 >  1:1239789/1‑150 (MQ=255)
                                                     gACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAACCCCGTCACGCTGCCAACCAGGATGGttt                                                                                               <  2:1026512/150‑1 (MQ=255)
                                                     gACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAACCCCGTCACGCTGCCAACCAGGATGGttt                                                                                               <  2:1026510/150‑1 (MQ=255)
                                                     gACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGtttt                                                                                              <  1:1682860/150‑1 (MQ=255)
                                                     gACAACAAGCACCACCAGCAGTGAGCATAGCCCTTCCAATTCCATAGCATTACCCATGGCGGTTATCGGAGACTCGGCCATCAAGCTCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAACCCCGTCACGCTGCCAACCATGATGGttt                                                                                               <  2:1026508/150‑1 (MQ=255)
                                                      acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                                                             >  2:786198/1‑150 (MQ=255)
                                                      acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                                                             >  2:194199/1‑150 (MQ=255)
                                                        acCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaa                                                                                           <  1:888710/150‑1 (MQ=255)
                                                         ccAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaa                                                                                          >  2:228669/1‑150 (MQ=255)
                                                         ccAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaa                                                                                          >  2:1101697/1‑150 (MQ=255)
                                                          cAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaaa                                                                                         <  2:1145438/150‑1 (MQ=255)
                                                          cAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaaa                                                                                         >  2:1551991/1‑150 (MQ=255)
                                                          cAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaaa                                                                                         <  2:1145436/150‑1 (MQ=255)
                                                             gCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                                                      <  1:12530/150‑1 (MQ=255)
                                                                 caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                                                  <  2:1088821/150‑1 (MQ=255)
                                                                 caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                                                  >  1:1159019/1‑150 (MQ=255)
                                                                 caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                                                  <  2:609971/150‑1 (MQ=255)
                                                                 caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                                                  >  2:744037/1‑150 (MQ=255)
                                                                 caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                                                  >  1:1159008/1‑150 (MQ=255)
                                                                   ccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTaa                                                                                <  1:497740/150‑1 (MQ=255)
                                                                    cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                                                               <  2:1217568/150‑1 (MQ=255)
                                                                     agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                                              >  1:1704932/1‑150 (MQ=255)
                                                                     agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                                              <  2:741480/150‑1 (MQ=255)
                                                                     agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                                              >  1:209408/1‑150 (MQ=255)
                                                                     agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                                              >  1:436350/1‑150 (MQ=255)
                                                                          tGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTc                                                                         >  2:329332/1‑150 (MQ=255)
                                                                                  gCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAg                                                                 >  2:153897/1‑150 (MQ=255)
                                                                                      tACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTg                                                             >  2:856374/1‑150 (MQ=255)
                                                                                                aTCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAAc                                                   >  2:627500/1‑150 (MQ=255)
                                                                                                aTCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAAc                                                   >  1:1136115/1‑150 (MQ=255)
                                                                                                 tCGCATTACCCATGGCGTTATGCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCATTGATGCCAGAT‑AAAAAACCCCGTCACGCTGCCAACCAGGTTGGTTTTCAAAATTGCTTAAAGCAACTCCGCGCCAGTCTGAGCCGAGAACa                                                  <  1:1594836/150‑1 (MQ=255)
                                                                                                  cGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaa                                                                                           >  1:828171/1‑108 (MQ=255)
                                                                                                  cGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaa                                                                                           <  2:828171/108‑1 (MQ=255)
                                                                                                  cGCATTACCCATGGCGGTAATCAGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACAt                                                 >  2:1399896/1‑150 (MQ=255)
                                                                                                         cccATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa                                          >  1:1156671/1‑150 (MQ=255)
                                                                                                         cccATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa                                          <  2:235263/150‑1 (MQ=255)
                                                                                                          ccATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTaa                                         >  2:114514/1‑150 (MQ=255)
                                                                                                               gCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGccc                                    <  2:1315921/150‑1 (MQ=255)
                                                                                                                        ggAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCggg                           >  1:644973/1‑150 (MQ=255)
                                                                                                                         gagaCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGt                          <  1:416181/150‑1 (MQ=255)
                                                                                                                         gagaCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGt                          <  2:820218/150‑1 (MQ=255)
                                                                                                                         gagaCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGt                          <  2:820188/150‑1 (MQ=255)
                                                                                                                         gagaCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGt                          <  1:1008338/150‑1 (MQ=255)
                                                                                                                          agaCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTGGCGGGtt                         <  2:1546016/150‑1 (MQ=255)
                                                                                                                            aCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGa                       <  2:338389/150‑1 (MQ=255)
                                                                                                                            aCTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGa                       <  1:1138020/150‑1 (MQ=255)
                                                                                                                              tCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCATGATGGTTTTCAAAATAGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGt                     >  1:255272/1‑150 (MQ=255)
                                                                                                                                   cATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCCGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGa                <  2:1289743/150‑1 (MQ=255)
                                                                                                                                   cATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGa                <  2:1239789/150‑1 (MQ=255)
                                                                                                                                   cATCAAGCGCATCATCCGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGa                <  2:1651076/150‑1 (MQ=255)
                                                                                                                                    aTCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGaa               >  1:639076/1‑150 (MQ=255)
                                                                                                                                         gcgcATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACtt          <  2:1222788/150‑1 (MQ=255)
                                                                                                                                          cgcATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTc         <  1:787737/150‑1 (MQ=255)
                                                                                                                                          cgcATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTc         <  1:787690/150‑1 (MQ=255)
                                                                                                                                           gcATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATTGTTTTCAAAATTGATTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGTTTGAGTTTGGAAAACTTCg        <  1:52902/150‑1 (MQ=255)
                                                                                                                                           gcATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCg        <  1:52886/150‑1 (MQ=255)
                                                                                                                                                 tcTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCa  >  1:356564/1‑150 (MQ=255)
                                                                                                                                                 tcTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCa  >  1:356453/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAG‑AAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCA  >  NC_000913/1965539‑1965833

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