Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,854,897 A→G 18.4% E326G (GAG→GGG)  emrD → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,854,8970AG18.4% 103.6 / 2.1 60E326G (GAG→GGG) emrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base G (5/9);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATC  >  NC_000913/3854663‑3855131
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
tGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCg                                                                                                                                                                                                                                             >  2:1121496/1‑239 (MQ=255)
  tttttttggccgtttcggtgggtgTTTTTTGGCGCATGGTTTTCCGGACGCCCCAATAAACGCTTCTCCACTTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGtt                                                                                                                                                                                                                                           <  1:22836/216‑1 (MQ=255)
     aTTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGTTGCCGGGATGCAGTTTCCGCTCGCGACCCGCGGCCGGATGGGGCCGTTcc                                                                                                                                                                                                                                         >  2:900745/1‑238 (MQ=255)
     aTTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCTGGGATGCTGTTTCCGCTGGCGACCCGCGTCCCGATGGAGCGGTTcc                                                                                                                                                                                                                                         >  2:1188621/1‑238 (MQ=255)
     aTTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGGGCCGGGATGCTGTTTCCCCTGGCGACCCACGGCCCGTTGGGGCCCTTcc                                                                                                                                                                                                                                         >  2:1321405/1‑238 (MQ=255)
     aTTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGGGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGGGCCGGGAGGCTGTTTCCGCTGGCGCCCAGCGGCGCGATGGAGcccgtc                                                                                                                                                                                                                                          >  2:400435/1‑233 (MQ=255)
                     cAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGGGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGGGCCGTTCCCCTTCCTGGCGGGCACg                                                                                                                                                                                                                       >  2:740344/1‑239 (MQ=255)
                          ttttttGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGGGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGCCCGCCGGCGCGATGGGGCCGTTCCCCTTCCTGGCGGGGCGGGCTgg                                                                                                                                                                                                                  >  2:1299640/1‑239 (MQ=255)
                                    cATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGGCTGGACGCTGCTCGTTCCCGCCCCGCTGTTCTTTTTCGGTGCCGGGGTGCTGTTTGCGCGGGCGCCCACGCGCGGCATGGGTGCGTTCCCCTTCCTGGCGGGGCCCGCTggggggcggggg                                                                                                                                                                                                        >  1:1021677/1‑229 (MQ=255)
                                     ttgttttcccggcgctCCCAATAAACCTTTCTCCCCGTTATGGGGGGCTGTGTTTTTCTGGCGCCTGCTGGGTGGTTTGCGGATGGGGTCCCCGGAGTGTTTTGGGGTTATGATGGTCGGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTcg                                                                                                                                                                                                       <  1:1013461/224‑1 (MQ=255)
                                      tGGTTTGCCGGACGCCCCAATAAACGGTTCTCCAGGTTAATGGGGCAGTCGGTTATCTGCGGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGAAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTcgg                                                                                                                                                                                                      <  2:816321/239‑1 (MQ=255)
                                                 aCGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACgg                                                                                                                                                                                                                      >  2:1220650/1‑212 (MQ=255)
                                                   gTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCCGGTCGGCGGGGTGCAAAAc                                                                                                                                                                                         >  2:275838/1‑239 (MQ=255)
                                                   gTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAACGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCAAGGCTGGGCCCCTGGCCGGCGGGCTGCcaaa                                                                                                                                                                                          >  2:451051/1‑238 (MQ=255)
                                                   gTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGGTTGGCGTGATGGATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCCCTTGCCACCAGCGGGCCGATGGAGCCGTTCCCCTTCCTTGCCGGGCCCGGGGGGCGCCCGGGTGGGGGGCGCGCaaaa                                                                                                                                                                                          >  1:785316/1‑239 (MQ=255)
                                                             aCGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGCCGGGGACGGCTGGCGCGCGGGTCGGGGGGCTGCAAAACATTGGTTCCg                                                                                                                                                                               >  2:923726/1‑239 (MQ=255)
                                                                           tAAGGTGGCAGTGGGTTAACGGCTGCCTGCTGGCTGGCTTGCTGATGTGGACCCCCGACTGGTTTGGCGTGATGAATTTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCg                                                                                                                                                               <  2:982795/239‑1 (MQ=255)
                                                                               tgtgGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGCGTGCTGGCGTGGc                                                                                                                                                              >  2:1181745/1‑237 (MQ=255)
                                                                               tgtgGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCGGTTTCCGCTGGCGACCGGCGGCGCGGGGGAGCCGTTCCCCTTCCTGGCGGGCAAGGCCGGCGCGCCGGGCGGGGGGCCGGAAAACATTGGGTCCgggggggggggggcggtc                                                                                                                                                              >  1:1325964/1‑223 (MQ=255)
                                                                                                 gctgcCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTTGCGGGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTCCGATGTTGCCGcca                                                                                                                                         >  2:578473/1‑237 (MQ=255)
                                                                                                 gctgcCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGGGGTCTGCCACACCTTGGTTCCGGCGTGCTGGGCTCGCTTTCTCCGATGTTGCCCCa                                                                                                                                          >  2:956633/1‑238 (MQ=255)
                                                                                                 gctgcCTGCGGGCTGGCTTGCTGATGTGGATCCCCGTCTGGGTTGCGGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCTCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCaa                                                                                                                                         <  2:111736/239‑1 (MQ=255)
                                                                                                  ctgcctgcGGGCTGGCTTGCTGATGGGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCaaa                                                                                                                                        <  2:108532/239‑1 (MQ=255)
                                                                                                     cctgcTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGGCGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGCCCCGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGGCCGGCCGGGCCCCTGGGCGGCGGGCTGCAAAAACTTGGTGCGCGCGGGGTGGCGGTCGTCTCCGGGATGTGGCGGCCAAccc                                                                                                                                     >  1:1022316/1‑238 (MQ=255)
                                                                                                             ctggctTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAggg                                                                                                                              <  2:680663/238‑1 (MQ=255)
                                                                                                               ggctTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGCCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGAGGTTGCCGCAAACCGGTCAGGGCa                                                                                                                            >  2:832369/1‑238 (MQ=255)
                                                                                                                 ctTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGAGCTGGCGTCCCTCTCTCCGAGTGTTCCCCAAACCGGTCAGGGCAccc                                                                                                                         >  2:952144/1‑239 (MQ=255)
                                                                                                                               tccccggtgggTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTtgat                                                                                                            <  2:955883/229‑1 (MQ=255)
                                                                                                                                  ccGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGCTGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGCCGGGCGCGCGGGCCGGGGGGCTGCAAACAATTGGTTCCGGGGTGCTGGCGCCCCCCTCTGCGATGTTTCGCCAAACCGGCCCGGGCAGCCCGGGGgggttggttg                                                                                                          >  1:1028284/1‑233 (MQ=255)
                                                                                                                                     aCTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGGCCGGCTGGCGCGCTGGGCGGCGGTCCCCCAAACCTTGGTTCCCGCGGGGCTGGGTCGCCCTCTGCGGGGGTGGCGCCAACCGGGCAGGGGCGCCCgggggggtggtgggcctta                                                                                                     >  1:1114245/1‑224 (MQ=255)
                                                                                                                                       tGGTTTGGCGTGATGAATGTCTGGACGCTCCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTaa                                                                                                    <  2:224468/238‑1 (MQ=255)
                                                                                                                                        ggTTTGGGGTGATGAAGGTCTGGACGCTGCTCGTTCCCGCCGCGCGTTTCTTTTTCGGTGCTGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTTGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATg                                                                                                  <  2:986035/239‑1 (MQ=255)
                                                                                                                                                    tgaATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAggg                                                                                                                              <  2:542400/199‑1 (MQ=255)
                                                                                                                                                       aTGCCTGGTGGCTGCTCCTTCCCGGCGCGCTGTTCTTGTTCGGGGCCGGGATGGTGTTCCCGCGGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtg                                                                                  <  2:447904/239‑1 (MQ=255)
                                                                                                                                                                        gTTCCCCCCGCGCTGTTCTTTTTGGGTGCCGGGATGCTGTTTCCGCGGGCGGCCGGCGGCGCGATGGGGCCGTTCCCCTTCCTGGGGGGCACGGCTGGCGGGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCt                                                                 <  1:646504/239‑1 (MQ=255)
                                                                                                                                                                           ccccctccgttgtgtttgtgcGGTGCCGGGTCGGGTTTTCCGGTGGCGACCGGCGGCGAGGGGGAGCCGTTCCCCGTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAAGGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGc                                                              <  1:1181745/219‑1 (MQ=255)
                                                                                                                                                                               ccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGGTTGTTGATCGTGCTTTGCTGGCGTCCCCTGGCGAc                                                           >  2:1136853/1‑238 (MQ=255)
                                                                                                                                                                               ccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAcg                                                          <  2:1028284/239‑1 (MQ=255)
                                                                                                                                                                                    cTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCCAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGGCGCGGAt                                                     >  2:1089393/1‑239 (MQ=255)
                                                                                                                                                                                      gTTCTTTTTCGGGGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATg                                                    <  2:287504/238‑1 (MQ=255)
                                                                                                                                                                                       ttCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCGGGCTGCCGCTGGCGACGCGGATGTc                                                  >  2:1194315/1‑239 (MQ=255)
                                                                                                                                                                                       ttCTTTTTCGGTGCCCGGGTGGTTTGTGCGCGGGGGACCGGGGGGGGGGTGGGGGCGTTCCCCTTCCGGGGGGGCCCGGTTGGCGCGCTGGGTGGGGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTc                                                  <  1:1136853/239‑1 (MQ=255)
                                                                                                                                                                                           tttctggTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCAt                                              <  2:505079/234‑1 (MQ=255)
                                                                                                                                                                                               cGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCCCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGGGCTGGCTGCCGCTGGCGACGCGGGTGGCGCTGCAgg                                          >  2:550033/1‑239 (MQ=255)
                                                                                                                                                                                                 gTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTCATGGGATTGTTGATCGTGCGGTGCTTGCGGCCGCTGGCGAccggggagtgggctccggg                                         >  2:251368/1‑220 (MQ=255)
                                                                                                                                                                                                    ccGGGATGGTGTTTCCGCTGGCGGCCAGCGGCGCGAGGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCa                                      <  2:1272561/238‑1 (MQ=255)
                                                                                                                                                                                                          gtctttttCCGCGGGCGGCCGTGGGCGGGCTGGGGCGGTTCCCCTGCCGGGCGGGCACGGTTGGCGCGGTTGTGCGCGGTCTTAAAACCGTGGGTTCCGGCGGGGTGGCGGCGCTCGCGTCGATTTTGCCGCAAACCGGTCAGGGCAGCCTGTGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTTGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCCGCCCGt                                <  1:1311042/233‑1 (MQ=255)
                                                                                                                                                                                                           ggtGTTTCCGCTGGCGACCAGGGGCGCGATGGGGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt                              <  2:301806/237‑1 (MQ=255)
                                                                                                                                                                                                                tcccGCGGGCGGCCGGGGGCGCGAGGGCGCCGTTCCCCTTCCTGGCGGGAACGGCTGGCGCGTTGGTCGGAGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgcg                         <  2:1022316/237‑1 (MQ=255)
                                                                                                                                                                                                                  ccGCGGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCa                       <  2:948545/239‑1 (MQ=255)
                                                                                                                                                                                                                     cggggggccgggggctggttGGGGCCGGGCCCGTTCCTGGCGGTCACGGTGGGCGGGCGGGGCGGGGGTCGGAAAAACAGTGGTTCCGTGGTGCGGTGGTCGCTCTATGCGAGTGTGCCGCAAGCCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGt                    <  1:850479/220‑1 (MQ=255)
                                                                                                                                                                                                                       ggCGACCAGCGGCGCGATGTAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGGTGGCCGGCGGTCAGCAAAAAATTTGGTCCGGCGGGCCGGCGGTCGCCCCTGCGATGGTGCCGCAAACCGGTCGAGGCAGCCCAGGGTGGTTTGCGACCTTTATGGGATTATATGACTTGGTGTGCTGGGTGCCGATGGCGCCCGGGAGGTGGCGTAAggggggacgccgatcacgacccgca                   >  1:184173/1‑217 (MQ=255)
                                                                                                                                                                                                                          gACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCg               <  2:1107582/239‑1 (MQ=255)
                                                                                                                                                                                                                            gccgggctccgttGGGGCTGTGCCCCTTCCGGGGGGGCCGGGCGGGCGCGCTGGGCGGCGGTCTGCAAACCCTTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGGCAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAAGGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCgca             <  1:740344/227‑1 (MQ=255)
                                                                                                                                                                                                                               gcggcgCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCCGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCgcagca          <  2:1114245/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gATGGAGCCGTTCCCCTCCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATc                                             <  2:562174/197‑1 (MQ=255)
                                                                                                                                                                                                                                      cttGGGGCCTTTCCCCTTCTTGGGGGGCACGGGTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCTTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGtcat   <  1:923726/237‑1 (MQ=255)
                                                                                                                                                                                                                                        tGGGGCCGTTCCCCTTCCTGGGGGGACCGGCTGGCGGGCTGGTCGGGGGTCGGCAAAACATTGGTTCCGGCGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGtcatc  <  2:1243312/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tGGGGCCGTTCCCCCTCCTGGCGGGCACGGCTGGCCGGCTGGTCGGTGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGtcatc  <  2:612017/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tGGAGCCGTTCCCCGTCCTGGCGGGCCCGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGtcatc  <  2:385669/238‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
TGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATC  >  NC_000913/3854663‑3855131

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