Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,388,041 T→G T50P (ACG→CCG)  aaeB ← aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,0410TG97.6% 121.9 / NA 41T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (20/20);  total (20/21)
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.81e-01

CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACC  >  NC_000913/3387795‑3388278
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
caATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcg                                                                                                                                                                                                                                            <  2:370444/251‑1 (MQ=255)
     aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                       <  2:250780/251‑1 (MQ=255)
                      ggCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCGCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCGGGACCGGCGGCAACAATCGCCGCTGTCa                                                                                                                                                                                                                      <  2:6695/251‑1 (MQ=255)
                       gCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCATCGCCGCTGTc                                                                                                                                                                                                                       >  2:450475/1‑249 (MQ=255)
                                       aCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTgg                                                                                                                                                                                                      <  2:61525/250‑1 (MQ=255)
                                                cTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCCCACCCGAATCATCAATAATGGT‑GCACGGATCATCGCACTGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCAGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGCGCGTTTcccgc                                                                                                                                                                                            >  2:524176/1‑248 (MQ=255)
                                                  nttCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTGCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGAGAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAg                                                                                                                                                                                             <  2:467152/247‑1 (MQ=255)
                                                      gTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAg                                                                                                                                                                                      <  1:321148/251‑1 (MQ=255)
                                                          cAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCACAAAGCCACGATAGCGAATAGCGCCAGAAGACGGTTCGCCTCCCGCAGCAAAGGCCGGACCTGCGGCAACAATCGCCGCTGTCAGTACCGCACAGCGTGGCGTTTCCAGCTGGAAGTgga                                                                                                                                                                                  <  2:156262/251‑3 (MQ=255)
                                                              gAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                              >  1:96928/1‑251 (MQ=255)
                                                               aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCca                                                                                                                                                                             <  1:264076/251‑1 (MQ=255)
                                                                      ccAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCaacaaac                                                                                                                                                                       >  1:35550/1‑250 (MQ=255)
                                                                                   ccGGCCCAGATACCGCACACCATAATCATCA‑TAATGGTTGCGCGGATCATCGCAATGATGAACACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                         <  2:492372/251‑1 (MQ=255)
                                                                                       cccAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCg                                                                                                                                                     >  1:445233/1‑251 (MQ=255)
                                                                                       cccAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTACGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCAGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGCGGCGCTTCCAGCTGGACGTGAAAGCCAACACACAGCGCCAGGACAATGGCg                                                                                                                                                     >  2:277197/1‑251 (MQ=255)
                                                                                           gATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACTATAGCGAATAGCGCCAGAATACGGTGCGCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCTCTGTCAGTACCGCCCAGCGTGGCGTTTCCGGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCg                                                                                                                                                 <  1:217728/251‑1 (MQ=255)
                                                                                              aCAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGcc                                                                                                                                               <  1:247471/250‑1 (MQ=255)
                                                                                                aGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACACTGGCGGTCACCAGt                                                                                                                                            >  2:209118/1‑251 (MQ=255)
                                                                                                   acacCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGgttt                                                                                                                                         >  1:419936/1‑251 (MQ=255)
                                                                                                             catcaATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGc                                                                                                                               <  1:190292/251‑1 (MQ=255)
                                                                                                                 aataatGGT‑GCGCGGATCATCGCAATGATGATCACCACTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAat                                                                                                                           <  2:325899/251‑1 (MQ=255)
                                                                                                                     atGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGt                                                                                                                       >  1:169984/1‑251 (MQ=255)
                                                                                                                                            gatgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATAc                                                                                                >  1:246046/1‑251 (MQ=255)
                                                                                                                                                          cGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCACCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAAAAATCGCCGCTGTTAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACGCCAGTTTTAACGCAAAGAGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATAGTGAATAAc                                                                                  >  1:217607/1‑251 (MQ=255)
                                                                                                                                                              aaTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCaa                                                                               >  1:241787/1‑250 (MQ=255)
                                                                                                                                                                   aGCCAATAACTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTTCCGCAAAGCGAATATGGTGGGTAGACAA‑GGGGGAAATACCCATCGCGATTAACCAAACTccc                                                                         >  2:322404/1‑249 (MQ=255)
                                                                                                                                                                   aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCAc                                                                         <  2:7616/251‑1 (MQ=255)
                                                                                                                                                                    gCCAATAAATGTGCCGATGAAGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGTAAGTGAGAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGTATATGTGGGTTAG‑CAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACg                                                                        <  1:500291/251‑1 (MQ=255)
                                                                                                                                                                          aaaTGTGCCGATGATGCGCAAAAAGCCACGGTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCg                                                                  <  1:283924/251‑1 (MQ=255)
                                                                                                                                                                              cccgatctgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATg                                                               <  2:287223/243‑1 (MQ=255)
                                                                                                                                                                                      gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCAGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGGACAATGGCGGTCGCCAGGTTTAACGCAAAGCGAATATGGTGGTTAG‑CAAAGGAGAAAAAACACATCGTGATTAACCAAACTCACGCAGGCGATgagatagga                                                       >  1:287223/1‑242 (MQ=255)
                                                                                                                                                                                               aaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaagaa                                             <  2:241787/251‑1 (MQ=255)
                                                                                                                                                                                                aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATAACCATCGTGATTAACCAAACTCACGCAGGCGGTGGGCCATTTTACGgaagaa                                             >  1:216050/1‑250 (MQ=255)
                                                                                                                                                                                                aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaagaa                                             >  1:275951/1‑250 (MQ=255)
                                                                                                                                                                                                 aGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCTCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACg                                           <  2:7863/251‑1 (MQ=255)
                                                                                                                                                                                                        aTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGatccctg                                    <  2:486529/251‑7 (MQ=255)
                                                                                                                                                                                                             gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCttt                               <  1:322404/251‑1 (MQ=255)
                                                                                                                                                                                                                                ttCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTTGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGcc            >  1:229900/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATAGCCGCTGTCAGTACCGCACAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGGCAGTTTTACACCAAAGCGAATTTATTGGATAG‑CAAAGGCGAAAATACCCACCGTGATTAACCAAACTCACGCAGGAGAAGGGCCATTTTCCGGAAGAAAGAAACCAGGCCTTCCTCGCCATATTGttttcc            >  1:255447/1‑246 (MQ=255)
                                                                                                                                                                                                                                         cGCAGCAAAGGCCGGACCGGCGGAAACAATCGCCGCTGTCAGTACAGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGAACAATGGCGGCAGCCAGTTTTTCCGCCAAACGAATATGTTGGTTAG‑CAATTGAGAAAATACCCATCGTGATTAACCCAACCCCCGCCGGCCATGGGGCATTTTTCGGAAGAACGAAACCTGGATTTCGCCGCACTCCTGTTTGccccttaccac   >  2:297839/1‑242 (MQ=255)
                                                                                                                                                                                                                                          gcagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGGGAAAATACCCATCGTGATTAACCAAACTTACGCAGGCGATGGGGCATTTTACGGAAGAACGAATTCTTGCTTTCGTCGCGGTCTTGTTTGCCAGTGaccacc  >  1:35385/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
CAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGT‑GCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAG‑CAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACC  >  NC_000913/3387795‑3388278

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