Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,567,982 A→G 5.7% S180P (TCA→CCA)  mdtM ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,567,9820AG 5.7% 219.9 / 3.3 70S180P (TCA→CCA) mdtMmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base G (2/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

ATATAGCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCTTTTGTCTGTCC  >  NC_000913/4567834‑4568129
                                                                                                                                                    |                                                                                                                                                   
atatAGCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATga                                                                                                                                                     <  4:216244/149‑1 (MQ=255)
atatAGCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATga                                                                                                                                                     <  3:28466/149‑1 (MQ=255)
    aGCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATa                                                                                                                                                 >  4:235645/1‑149 (MQ=255)
     gCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATaa                                                                                                                                                >  3:82166/1‑149 (MQ=255)
       ttAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATaaaa                                                                                                                                              >  4:580993/1‑149 (MQ=255)
            gagaTGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCAt                                                                                                                                         >  2:284831/1‑149 (MQ=255)
             agaTGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATa                                                                                                                                        >  3:428262/1‑149 (MQ=255)
             agaTGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATa                                                                                                                                        >  1:565710/1‑149 (MQ=255)
                    tgctgcGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCaa                                                                                                                                 <  4:195988/149‑1 (MQ=255)
                     gctgcGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAat                                                                                                                                <  3:295615/149‑1 (MQ=255)
                     gctgcGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAat                                                                                                                                <  3:425669/149‑1 (MQ=255)
                        gcgcCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAatgat                                                                                                                             <  3:40408/149‑1 (MQ=255)
                         cgcCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAatgatg                                                                                                                            <  2:153751/149‑1 (MQ=255)
                                gagGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAat                                                                                                                                >  3:223524/1‑138 (MQ=255)
                                gagGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAat                                                                                                                                <  4:223524/138‑1 (MQ=255)
                                gagGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCaaaaa                                                                                                                     <  3:576800/149‑1 (MQ=255)
                                  ggAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAgg                                                                                                                   >  1:461420/1‑149 (MQ=255)
                                   gAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGa                                                                                                                  <  2:565710/149‑1 (MQ=255)
                                           gATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAg                                                                                                          >  4:382143/1‑149 (MQ=255)
                                            aTTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCGCCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGt                                                                                                         <  4:561419/149‑1 (MQ=255)
                                            aTTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGt                                                                                                         >  1:59471/1‑149 (MQ=255)
                                            aTTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGt                                                                                                         <  3:469363/149‑1 (MQ=255)
                                             ttGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTg                                                                                                        <  4:141833/149‑1 (MQ=255)
                                             ttGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTg                                                                                                        >  1:52999/1‑149 (MQ=255)
                                                cAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCAt                                                                                                     <  1:284831/149‑1 (MQ=255)
                                                 aaaaGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATa                                                                                                    >  3:341385/1‑149 (MQ=255)
                                                  aaaGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATaa                                                                                                   >  4:239777/1‑149 (MQ=255)
                                                     gACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGGGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGt                                                                                                >  3:501716/1‑149 (MQ=255)
                                                         ttACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCAt                                                                                            <  2:300549/149‑1 (MQ=255)
                                                             gAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCaga                                                                                        <  1:367180/149‑1 (MQ=255)
                                                               aaaTCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGc                                                                                      <  1:3556/149‑1 (MQ=255)
                                                                aaTCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCt                                                                                     >  3:445028/1‑149 (MQ=255)
                                                                     cgcAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGcc                                                                                >  1:575344/1‑149 (MQ=255)
                                                                      gcAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCg                                                                               >  4:326101/1‑149 (MQ=255)
                                                                      gcAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCg                                                                               >  1:537369/1‑149 (MQ=255)
                                                                        aaGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAatgat                                                                                                                             >  3:546001/1‑101 (MQ=255)
                                                                        aaGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAatgat                                                                                                                             <  4:546001/101‑1 (MQ=255)
                                                                         aGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGaa                                                                            >  2:479487/1‑149 (MQ=255)
                                                                          gACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGGGATAAAACCCATAACCGCAatga                                                                                                                              <  2:595862/98‑1 (MQ=255)
                                                                          gACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGGGATAAAACCCATAACCGCAatga                                                                                                                              >  1:595862/1‑98 (MQ=255)
                                                                               ttttGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCggg                                                                      >  3:443142/1‑149 (MQ=255)
                                                                                        tAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATc                                                             <  2:291186/149‑1 (MQ=255)
                                                                                        tAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATc                                                             <  3:382143/149‑1 (MQ=255)
                                                                                           aCGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACATATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          >  2:145838/1‑149 (MQ=255)
                                                                                           aCGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGGGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          <  4:501716/149‑1 (MQ=255)
                                                                                           aCGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          <  4:235948/149‑1 (MQ=255)
                                                                                           aCGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          >  4:91664/1‑149 (MQ=255)
                                                                                            cGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgcg                                                         >  4:180963/1‑149 (MQ=255)
                                                                                             ggAACCGCGCCGCGCTTCACCGTCTATGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCgg                                                                       >  1:380871/1‑134 (MQ=255)
                                                                                             ggAACCGCGCCGCGCTTCACCGTCTATGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCgg                                                                       <  2:380871/134‑1 (MQ=255)
                                                                                                  cgcgccgcgcTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCa                                                   >  3:441008/1‑149 (MQ=255)
                                                                                                     gccgcgcTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTa                                                <  2:23020/149‑1 (MQ=255)
                                                                                                      ccgcgcTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTAc                                               <  2:566416/149‑1 (MQ=255)
                                                                                                      ccgcgcTTCACCGTCTCTGGCATCGCCAACAGTAAGCAAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTAc                                               <  3:495799/149‑1 (MQ=255)
                                                                                                            ttCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGa                                         <  1:62809/149‑1 (MQ=255)
                                                                                                             tCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAg                                        >  2:178821/1‑149 (MQ=255)
                                                                                                               aCCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGt                                      <  4:190272/149‑1 (MQ=255)
                                                                                                                 cGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          >  3:286636/1‑127 (MQ=255)
                                                                                                                 cGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGgc                                                          <  4:286636/127‑1 (MQ=255)
                                                                                                                      ctGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATc                               <  3:144232/149‑1 (MQ=255)
                                                                                                                        ggCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGc                             <  2:575344/149‑1 (MQ=255)
                                                                                                                        ggCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGc                             <  1:389433/149‑1 (MQ=255)
                                                                                                                         gCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGCCTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGcc                            >  1:33847/1‑149 (MQ=255)
                                                                                                                               ccAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAAc                      <  2:537369/149‑1 (MQ=255)
                                                                                                                                 aaCAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACtt                    <  3:326101/149‑1 (MQ=255)
                                                                                                                                    aGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGAt                 <  3:235645/149‑1 (MQ=255)
                                                                                                                                         gCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCtt            >  2:372065/1‑149 (MQ=255)
                                                                                                                                         gCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCtt            >  2:333768/1‑149 (MQ=255)
                                                                                                                                          ccAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGCTCCCttt           <  2:320085/149‑1 (MQ=255)
                                                                                                                                           caacaaATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCtttt          <  1:106161/149‑1 (MQ=255)
                                                                                                                                             acaaATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTa                                                <  4:551616/109‑1 (MQ=255)
                                                                                                                                             acaaATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTa                                                >  3:551616/1‑109 (MQ=255)
                                                                                                                                                   gagaTAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCTTTTGTCTGTcc  <  4:428262/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATATAGCTTAAAGAGATGGTTGCTGCGCCAAAGAGGAACAGCCGATTGCAAAAGACATTACGAAAATCGCGCAAGACGCTTTTGGCGCTAAACGGAACCGCGCCGCGCTTCACCGTCTCTGGCATCGCCAACAGTAAGCCAACAAATGAGATAAAACCCATAACCGCAATGATGGCAAAAAGGACTTTCCAGTGCATAAAGTGCATCAGAGCTGCGCCGGAAAGCGGGCCGATAATCGGCGCAATCAGTACGATGGAGGTGATAATCGCCATCAACTTGATCCCTTTTGTCTGTCC  >  NC_000913/4567834‑4568129

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