Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,577,060 T→A 5.3% E396D (GAA→GAT ttgC ← efflux pump outer membrane protein TtgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,577,0600TA 5.3% 235.6 / 2.5 76E396D (GAA→GATttgCefflux pump outer membrane protein TtgC
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base A (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

GCGCTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTGAAGATCGG  >  NC_002947/1576918‑1577203
                                                                                                                                               |                                                                                                                                                
gcgcTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCCGCGCTGGGCCTGCAACTGCTCTTCGaa                                                                                                                                            <  1:17283/149‑1 (MQ=255)
 cgcTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAgg                                                                                                                                          >  4:115351/1‑149 (MQ=255)
 cgcTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAgg                                                                                                                                          >  1:98551/1‑149 (MQ=255)
  gcTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCTGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGTTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAgg                                                                                                                                          >  2:311363/1‑148 (MQ=255)
   cTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTc                                                                                                                                        <  3:64743/149‑1 (MQ=255)
    ttGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCt                                                                                                                                       >  3:126756/1‑149 (MQ=255)
    ttGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCt                                                                                                                                       >  6:238561/1‑149 (MQ=255)
      gCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCtt                                                                                                                                      <  3:251358/148‑1 (MQ=255)
      gCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCtt                                                                                                                                      <  1:88701/148‑1 (MQ=255)
        tGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCg                                                                                                                                   >  2:264412/1‑149 (MQ=255)
            cACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGAAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTgcc                                                                                                                               >  3:82730/1‑149 (MQ=255)
              cTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTgccgc                                                                                                                             <  3:115351/149‑1 (MQ=255)
               tGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTgccgc                                                                                                                             >  2:89844/1‑148 (MQ=255)
               tGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTCCGTgccgc                                                                                                                             >  7:46129/1‑148 (MQ=255)
                 tagaGGTTGCGTTGGGCATCGAGGATGGTCAGGTTGCTGTCAATCCCGATGCGGTTACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGACGGTCTTGCGTGCCGCCAg                                                                                                                          <  7:71533/146‑1 (MQ=255)
                 aaCAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCa                                                                                                                           <  4:334983/148‑1 (MQ=255)
                 aaCAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAg                                                                                                                          <  4:159499/149‑1 (MQ=255)
                  aCAGGTTGCGTTGGGCATCTAGGAAGGGCTGGTTGCTGTCAATCCCGATGCGTTAACGGCGTTCGGCCAAGCGTTAGTAGTCTTAGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAg                                                                                                                          <  1:311363/148‑1 (MQ=255)
                     ggTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCa                                                                                                                      >  6:117013/1‑149 (MQ=255)
                          cGTGGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCGCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGga                                                                                                                 <  8:355238/149‑1 (MQ=255)
                            ttGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGgaga                                                                                                               <  1:89844/149‑1 (MQ=255)
                                 cATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCATCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTcc                                                                                                          >  3:150324/1‑149 (MQ=255)
                                 cATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCATCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTcc                                                                                                          >  3:156092/1‑149 (MQ=255)
                                    cGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTcc                                                                                                          <  1:264412/147‑1 (MQ=255)
                                     gAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTgg                                                                                                       <  6:260088/148‑1 (MQ=255)
                                     gAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTgg                                                                                                       >  3:258616/1‑148 (MQ=255)
                                       ggAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTGCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAg                                                                                                    <  7:108818/149‑1 (MQ=255)
                                          aGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGcc                                                                                                 <  1:389184/149‑1 (MQ=255)
                                               aGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCt                                                                                             <  5:134335/148‑1 (MQ=255)
                                                 gTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGa                                                                                          >  4:210195/1‑149 (MQ=255)
                                                    gCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTc                                                                                                                                        <  3:239288/100‑1 (MQ=255)
                                                    gCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTc                                                                                                                                        >  4:239288/1‑100 (MQ=255)
                                                          aaTCCCGATGCGGTCACGGCGTTCGGCCAAGCGGTAGTAGTCCGGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGGCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGtttatt                                                                                 <  6:82470/149‑1 (MQ=255)
                                                          aaTCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGtttttt                                                                                 <  4:126756/149‑1 (MQ=255)
                                                           aTCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGc                                                                                                                                    <  5:168334/97‑1 (MQ=255)
                                                           aTCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGc                                                                                                                                    >  6:168334/1‑97 (MQ=255)
                                                                  tGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTg                                                                         <  1:105172/149‑1 (MQ=255)
                                                                     ggTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGc                                                                       >  2:192122/1‑148 (MQ=255)
                                                                       tAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCa                                                                                                                           >  1:10376/1‑94 (MQ=255)
                                                                       tAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCa                                                                                                                           <  2:10376/94‑1 (MQ=255)
                                                                        aaCGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACg                                                                   >  3:131464/1‑149 (MQ=255)
                                                                         aCGGCGGTCTGCCAATCGGTAGTAGTCCTGGTTCTCCTGCACCATGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTGTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCGGTCTGGATGGTTTTTTCGTACTTGGCGACGt                                                                  <  8:46129/149‑1 (MQ=255)
                                                                           ggCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCATCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTg                                                                <  4:156092/149‑1 (MQ=255)
                                                                           ggCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCATCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTg                                                                <  4:150324/149‑1 (MQ=255)
                                                                              gTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATg                                                             >  5:229085/1‑149 (MQ=255)
                                                                                  ggCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCt                                                         <  4:258616/149‑1 (MQ=255)
                                                                                   gCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGgaga                                                                                                               >  7:275610/1‑94 (MQ=255)
                                                                                   gCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGgaga                                                                                                               <  8:275610/94‑1 (MQ=255)
                                                                                   gCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCtt                                                        >  7:11981/1‑149 (MQ=255)
                                                                                    ccAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTc                                                       >  5:191381/1‑149 (MQ=255)
                                                                                        gCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGa                                                   <  8:35865/149‑1 (MQ=255)
                                                                                          ggTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGAtt                                                 <  2:98551/149‑1 (MQ=255)
                                                                                           gtagtaGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGTTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGAttt                                                >  3:170965/1‑149 (MQ=255)
                                                                                            tagtagTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGAtttt                                               >  1:377464/1‑149 (MQ=255)
                                                                                             agtagtCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTc                                              <  4:82730/149‑1 (MQ=255)
                                                                                              gtagtCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTc                                              <  5:238561/148‑1 (MQ=255)
                                                                                               tagtCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTc                                              >  4:175259/1‑147 (MQ=255)
                                                                                                 gtCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGag                                           <  3:210195/148‑1 (MQ=255)
                                                                                                    cTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAGGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGagtagt                                       >  7:203123/1‑149 (MQ=255)
                                                                                                       gTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTc                                      >  7:311335/1‑147 (MQ=255)
                                                                                                          cGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAgg                                  <  3:175259/148‑1 (MQ=255)
                                                                                                            ccTGCACCAGGTCGCGCTGGTCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAggctg                               <  2:105433/149‑1 (MQ=255)
                                                                                                            ccTGCACCAGCTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAggctg                               >  3:324753/1‑149 (MQ=255)
                                                                                                               gCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGa                                            <  8:229715/133‑1 (MQ=255)
                                                                                                               gCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGa                                            >  7:229715/1‑133 (MQ=255)
                                                                                                                cACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAggctggct                            <  5:117013/148‑1 (MQ=255)
                                                                                                                cACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCtt                           >  4:364209/1‑149 (MQ=255)
                                                                                                                 aCCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCttt                          >  6:345636/1‑149 (MQ=255)
                                                                                                                    aGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAGGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAg                       <  8:203123/149‑1 (MQ=255)
                                                                                                                     ggTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAgg                      <  3:364209/149‑1 (MQ=255)
                                                                                                                           cgcTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGACATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGcc                 >  3:6552/1‑149 (MQ=255)
                                                                                                                               gggCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCg             <  5:345636/149‑1 (MQ=255)
                                                                                                                                  ccTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGtt           >  6:356204/1‑148 (MQ=255)
                                                                                                                                   cTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTg          >  8:38798/1‑148 (MQ=255)
                                                                                                                                       aaCTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTTCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTGAAGAt     <  1:192122/149‑1 (MQ=255)
                                                                                                                                          tGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTGAAGATCgg  >  1:147132/1‑149 (MQ=255)
                                                                                                                                            ctctTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTg          >  3:155371/1‑139 (MQ=255)
                                                                                                                                            ctctTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTg          >  1:176538/1‑139 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCGCTTGCTGGGCACTGAACAGGTTGCGTTGGGCATCGAGGAAGGTCAGGTTGCTGTCAATCCCGATGCGGTAACGGCGTTCGGCCAAGCGGTAGTAGTCCTGGTTCGCCTGCACCAGGTCGCGCTGGGCCTGCAACTGCTCTTCGAAGGTCTTGCGTGCCGCCAGGCCATCGGAGACTTCCTGGAAGGCCGTCTGGATGGTTTTTTCGTACTTGGCGACGTTGATGTCCTTCTGGATTTTCGAGTAGTCCAGGCTGGCTTTGAGGCTGCCGGCGTTGAAGATCGG  >  NC_002947/1576918‑1577203

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