Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2680TA40.4% 65.9 / 30.2 57L24F (TTA→TTT) emrKmultidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (20/14);  new base A (0/23);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGC  >  NC_000913/2483118‑2483417
                                                                                                                                                       |                                                                                                                                                      
gTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCaataaat                                                                                                                                                        <  1:714849/151‑1 (MQ=255)
  tCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCaataaata                                                                                                                                                       >  1:851330/1‑150 (MQ=255)
          ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTa                                                                                                                                                <  2:321475/148‑1 (MQ=255)
          ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCg                                                                                                                                             >  1:201665/1‑151 (MQ=255)
            aTGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGcc                                                                                                                                           >  2:972635/1‑151 (MQ=255)
              gATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCa                                                                                                                                          <  2:1148235/150‑1 (MQ=255)
                    cGACAGTGACCCTACCTGAGTCTTGTGCAAAAATTGGATCTGCATCCCCCGGGACATAGGCGTCATCTGTACTAATAATTTCTTCTAATCCCATGGCCA‑AAAAGGCATAGGCCCCTGAAAACGAAAAAAAAAAAA‑CTACCGCCAATAAAGa                                                                                                                                   <  1:1166049/151‑1 (MQ=255)
                     gACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAg                                                                                                                                    >  1:399598/1‑149 (MQ=255)
                     gACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGa                                                                                                                                   >  2:491879/1‑150 (MQ=255)
                      aCAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGaa                                                                                                                                  >  2:9523/1‑150 (MQ=255)
                           gACACTACCGAATGCTTGTGAAGAAATTGGATCTGCTTCCCCGGTGAAAGGGGTGTAATCTGTACTAATCATGTCTTCTATTCCCATTGACC‑AAGAGGAATAGGCACTTGAAAACAAAATAAATAAAT‑CTACCGCCAATAAAg                                                                                                                                    <  1:530478/143‑1 (MQ=255)
                                tACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTttct                                                                                                                       >  2:2021/1‑151 (MQ=255)
                                 aCCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTttctt                                                                                                                      >  2:549951/1‑151 (MQ=255)
                                  ccTGAGATTTGTGCAAAAATTGGATCTCCTTTCCCCGGGACATAGGCTTCATCTGTACTAATATTGTCTTCTAATCCAATTGACC‑AATGGGAAAAGCCCCCTAAAAACGCAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTttct                                                                                                                       <  2:1143867/149‑1 (MQ=255)
                                     gagaCTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGt                                                                                                                  >  2:1145396/1‑151 (MQ=255)
                                         tttGTGGAAAAATTGTGCCTGCATCCCCCGTGACATAGGCTAAATCTGACATAAACATGTTTTAAAATTCCATTCCCA‑AATGGGCAGAGGCACCAAAAACCACAAAAAAAAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                                                                                               <  2:546956/149‑1 (MQ=255)
                                          ttGTGAAAAAATGTGATTTCCTTTCCCCGGGACATAGGGGTCATTTGTACTAATAATGTTTTCAAATCCCA‑TGACCAAAAAGAAATAGGCACCT‑AAAACCCAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                              <  1:809603/150‑1 (MQ=255)
                                            gtgaaaaaaTTGGATCTCCATTCCCCGTGAAATAGGCTTCATCTTTACTAATAATTTCTTCTAATCCCATTGACC‑AATAGGCATAGCCCCCTAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAAt                                                                                                            <  1:1063327/144‑1 (MQ=255)
                                                   atttGTATCTGCTTCCCCCGTGACATAGGCGTCATCTGTACTAATCATTTCTTCAAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTa                                                                                                    <  2:737471/149‑1 (MQ=255)
                                                    aTTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                                                                                                   >  1:1023294/1‑151 (MQ=255)
                                                          tCTGCTTCCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATACCATTGACA‑AATAGGAAAAGCCCCTTAAAAACGAAAAAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                <  1:271919/148‑1 (MQ=255)
                                                              cATTCCCAGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGGTAGAATGTTTTTTATTTGAATTa                                                                                           >  1:558819/1‑149 (MQ=255)
                                                                 cccccGACACAAAGGCTTCATCTGTATAAACAATGTCTTTAAATCCCATTACCC‑ATTAGAAAAAGGCACTTAAAACCTCAAAAAAAAATT‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                                                                                                  <  1:241002/138‑1 (MQ=255)
                                                                    ccGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc                                                                                   >  2:404705/1‑151 (MQ=255)
                                                                        aaaaTAGGCCTCATTTTTATAAATTATGTTTTCAATTACCATTAACA‑AAAGGGCATAGGCCTTTAAAATCCCAATAAAAAAAT‑GTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa                                                                                 <  1:1081351/146‑1 (MQ=255)
                                                                         aCATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt                                                                               >  2:167225/1‑150 (MQ=255)
                                                                         aCATAGGCGTCATCTATACAAAACATGTCTTCTAATCCCATGCACA‑AAAAGGAATAGACACTAAAAAACGCAAAAAATAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa                                                                                 <  1:294555/148‑1 (MQ=255)
                                                                                 gTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctct                                                                      >  2:92642/1‑151 (MQ=255)
                                                                                       tatattaataaTTTCTTTAAATCCCATTACCC‑AATGGACATAGTCCCCAGAAACCGAAAAAAAAAAAA‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctctc                                                                     <  2:1095843/141‑1 (MQ=255)
                                                                                       tGTACTAATCATGTTTTCAAATCCCATTGACA‑AAAGGGAATAGGCCCCTAAAACCGAAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAttt                                                                 <  2:1014290/150‑1 (MQ=255)
                                                                                       tGTACAAACAATGTCTTCAAATTCCATTCACC‑AAAAGGAATAGGCTCTTAAAAACAAAAAAAAAAAAT‑CTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCa                                                                    <  1:9521/147‑1 (MQ=255)
                                                                                          aCTAATCTTGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCCCCTAAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtct                                                               <  2:694482/149‑1 (MQ=255)
                                                                                             aatgagtttttctAATTCAATGGACA‑AAAAGGCATAGCCACCTGAAAACCAAATAAAAAAAC‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATa                                                           <  2:19983/141‑1 (MQ=255)
                                                                                             aaTCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACTTGAAAACGAAAAAAAAAAAT‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                                                <  1:273317/145‑1 (MQ=255)
                                                                                                  tgtttacaaatccaaTTGCCA‑AATAGGAAAAGGCACCTGAAACCAAAAAAAATAAAT‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGAt                                                        <  1:708562/134‑1 (MQ=255)
                                                                                                        ctAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAt                                               >  1:1167756/1‑151 (MQ=255)
                                                                                                        ctAATTCCATTGACC‑AAAAGGCATAGGCCCCTGAAACCAAAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGaa                                                     <  2:386349/145‑1 (MQ=255)
                                                                                                         tAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtg                                              >  2:643250/1‑151 (MQ=255)
                                                                                                          aaTTCCATTGCCC‑AAAAGGCATAGCCACCAAAAACCGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGa                                                      <  2:905596/142‑1 (MQ=255)
                                                                                                          aaTTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtg                                              >  2:1001640/1‑150 (MQ=255)
                                                                                                          aaTTCCATTGACA‑AATAGGCATAGGCCCCAGAAAACGAAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtg                                              <  2:183756/150‑1 (MQ=255)
                                                                                                              ccATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                         >  1:970295/1‑151 (MQ=255)
                                                                                                               cATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                         >  2:1170449/1‑150 (MQ=255)
                                                                                                                  tGACA‑AATAGGCATAGGCCCCTAAAAACACAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattgaat                                      <  2:863725/150‑1 (MQ=255)
                                                                                                                    accaaaTAGGCATAGCCACCAAAAAACAAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatg                                     <  2:505238/146‑1 (MQ=255)
                                                                                                                      ccaaaagtaataggcaCTTAAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                         <  2:1121130/135‑1 (MQ=255)
                                                                                                                        aataggaataggCCCCTAAAAACCAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaaa                                 <  2:479696/150‑1 (MQ=255)
                                                                                                                        aataggaataggCCCCTAAAAACCAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                                  <  1:197229/149‑1 (MQ=255)
                                                                                                                        aaTAGGCATAGGCACCTAAAAACGCAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAag                                <  2:866175/151‑1 (MQ=255)
                                                                                                                        aaTAGGCATAGCCACCTAAAAACGAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAag                                <  1:739851/151‑1 (MQ=255)
                                                                                                                         ataggaatagacACCTAAAAACAAAAAAAAAAAAT‑TTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAaga                               <  1:420926/151‑1 (MQ=255)
                                                                                                                           aggcataggcaCCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGc                             >  2:339358/1‑151 (MQ=255)
                                                                                                                            ggcataggccCCTAAAAACGCAAAAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGca                            <  2:808396/151‑1 (MQ=255)
                                                                                                                             gcataggcaCCTAAAAACACAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcat                           <  2:330309/151‑1 (MQ=255)
                                                                                                                             gaatgggcaCCTGAAAACAAAAAAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAagag                              <  2:1190563/143‑1 (MQ=255)
                                                                                                                                    caCCAGAAAACGAAAAAAATAAAC‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattcatt                      <  2:374482/149‑1 (MQ=255)
                                                                                                                                      ccTAAAAACGCAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttct                    <  2:320687/149‑1 (MQ=255)
                                                                                                                                        tGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat                >  1:167996/1‑151 (MQ=255)
                                                                                                                                          aaaaCAAAATAAAAAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat                <  1:1071171/149‑1 (MQ=255)
                                                                                                                                                 aataaataAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAacac       >  2:1109605/1‑151 (MQ=255)
                                                                                                                                                    aaataaaA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAg     >  1:425119/1‑150 (MQ=255)
                                                                                                                                                       taaaA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGc  >  1:242913/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACC‑AATAGGCATAGGCACCTGAAAACGCAATAAATAAAA‑CTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGC  >  NC_000913/2483118‑2483417

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