Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,4320TG24.6% 601.3 / 48.9 367L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (130/143);  new base G (66/23);  total (200/167)
Fisher's exact test for biased strand distribution p-value = 1.49e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 377 total aligned reads displayed.

TGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTC  >  NC_002947/4160285‑4160580
                                                                                                                                                   |                                                                                                                                                    
tgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTg                                                                                                                                                     >  1:3330749/1‑149 (MQ=255)
 gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGggcggc                                                                                                                                                    >  1:4539250/1‑149 (MQ=255)
 gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGgc                                                                                                                                                    >  1:9016131/1‑149 (MQ=255)
  tCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATcggcgcggcg                                                                                                                                                   >  1:3441348/1‑149 (MQ=255)
   cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCggcgcggcgc                                                                                                                                                  >  1:5815374/1‑149 (MQ=255)
     gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGgcgccg                                                                                                                                                >  1:381562/1‑147 (MQ=255)
     gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAg                                                                                                                                                <  1:7468104/149‑1 (MQ=255)
          ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                           <  2:6839383/149‑1 (MQ=255)
            tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGACAGTGTTCAAGCTCGACCGCGCTACCGCCTTGCTGACCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCAGCGCAGCGccccc                                                                                                                                         >  2:6495037/1‑149 (MQ=255)
             aGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                        <  2:5092897/149‑1 (MQ=255)
               tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGccccccca                                                                                                                                      >  2:8558470/1‑147 (MQ=255)
                cACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATGGGCGCTGCGCAGCGCCCCCCATa                                                                                                                                     <  2:185729/149‑1 (MQ=255)
                cACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATa                                                                                                                                     <  2:8489523/149‑1 (MQ=255)
                   cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCagacaa                                                                                                                                  >  2:6013064/1‑144 (MQ=255)
                    aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAg                                                                                                                                 <  2:2388929/149‑1 (MQ=255)
                          ctgctgATCGGCCTGTGGTGTGGATTGAATTTTTTCACGCTCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAt                                                                                                                           <  1:798798/149‑1 (MQ=255)
                           ggctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGCTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCTTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATg                                                                                                                          <  2:5315373/148‑1 (MQ=255)
                                    gcttgtggggtggCTTGATTGTGTTCACGCTCGCCCGCGATCCCGGCTTGCTGACCGCCGCCGGCCGCGCCTTGGGCGGCGCCGCCGCGGTGCTGGATTTCGCTTGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGa                                                                                                                 <  1:5048038/146‑1 (MQ=255)
                                    gCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGa                                                                                                                 <  2:4230317/149‑1 (MQ=255)
                                    gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCCTGGCCGTCGGa                                                                                                                 >  2:2078675/1‑149 (MQ=255)
                                    gCCTGTGGTGGGGCTTGATTGTGTTCAAGCTCGACCGCGATACGGCTTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCGGCCGCGGTGGTGGCTTTCGATTCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGa                                                                                                                 <  1:4944970/149‑1 (MQ=255)
                                     ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCATGGCAGTCGGaa                                                                                                                >  2:10462518/1‑149 (MQ=255)
                                     ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCCCCCCCCCAAAAAGCGCCCGGGCAGTCGGac                                                                                                                >  2:8223027/1‑148 (MQ=255)
                                       tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCCAAAAAGCCCCAAGGCCGTCGGTAGc                                                                                                              >  1:3052709/1‑149 (MQ=255)
                                         tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAATCCCCAAGGCACTCGGAAGCGt                                                                                                            >  2:4383017/1‑149 (MQ=255)
                                         tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGAAACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAACGCCAAGGCAGTCGGAACCGt                                                                                                            >  2:7533062/1‑149 (MQ=255)
                                         tggtggggCATGAATGTGTTCACGCTCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCTAATCGGCGCTGCGCTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGt                                                                                                            <  1:6960756/149‑1 (MQ=255)
                                                 cATGAAAGTGTTCAAGCTCGACCGCGATCCCGCCTTGCTGCCCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCGGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGt                                                                                                    <  1:9551370/149‑1 (MQ=255)
                                                   tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCACAAAAACCCCCAGGCCCTCGGGAGCGGAGTGGTGttg                                                                                                  >  1:6606593/1‑148 (MQ=255)
                                                   tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAGCGCCCAGGCAGTCGGAAGCGGAGTGcgggtc                                                                                                  >  1:894168/1‑144 (MQ=255)
                                                    gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                 >  2:10058042/1‑149 (MQ=255)
                                                    gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAGCGCCAAGGCCCTCGGAAGCGTAGTGCTGttag                                                                                                 >  1:9350827/1‑147 (MQ=255)
                                                       atttttCAAGCTCTACCGCGATACCGCATTGCTGACCGCCGCCGGAAGCGACATCTGCGGCGCCGCCGCCGTTCTGGAATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                              <  2:1845921/145‑1 (MQ=255)
                                                       aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAATCGCCATGGCAGTCGGAAGCGGAGTGCTGTTAGGTa                                                                                              >  1:9616738/1‑149 (MQ=255)
                                                        gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                             <  2:9105804/149‑1 (MQ=255)
                                                        gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                             <  2:7702907/149‑1 (MQ=255)
                                                        gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAACGCCAAGGCACTCGGAAGCGGAGTGCTGGTGGGTaa                                                                                             >  1:7274557/1‑148 (MQ=255)
                                                                gCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGa                                                                                     <  2:3601320/149‑1 (MQ=255)
                                                                gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCGAAAAGGGCCAAGGCCGGCGGGAGCGTAGTGCTGTTCGGTGCCGTGTCgg                                                                                     >  1:2787316/1‑148 (MQ=255)
                                                                  aggcccGCTATACCGACTTGCTGACCGCCGCCGGCCGCGCAATCGGCGGCGCCGCCGCGGTGGTGGAATTCGACTCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATg                                                                                   <  1:1833123/145‑1 (MQ=255)
                                                                    gcccGCGATACGGCCTTGCTGCCCGCCGCCGGCAGCGCCATGGGCGGCGCCGCCGCCGTGGTGGCATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGtt                                                                                 <  1:2316774/147‑1 (MQ=255)
                                                                     aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCCGCAAAGCGCCAGGGCACGCGGAACCGGCGGGCGGGTCGGGACCGTGTCgcgggtc                                                                                >  2:8956496/1‑143 (MQ=255)
                                                                     aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCCCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                >  2:7853973/1‑149 (MQ=255)
                                                                      ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                               <  2:5003426/149‑1 (MQ=255)
                                                                      ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAACGCCCAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTGGATGttac                                                                               >  1:8560678/1‑147 (MQ=255)
                                                                         gtaccccGCGTTGCTCCCCGCCGCGGGCCGCGCCATGGGCGGCGCCGCCGCGGGGGTGGAATTCTAATCGGCGCTGCCCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                            <  1:2078675/144‑1 (MQ=255)
                                                                         cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAACCGCCCGGGCACGCGGGAGCGTAGTGCGGTTCGGTTCCGTGTCGATGTTCCttt                                                                            >  1:7615227/1‑149 (MQ=255)
                                                                          gATACCGCCTTGCTTACCGCCGCCGGCTGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCACCATAAAAGCGCCATGTCAGTCGGAAGCGTAGTGCTTTTCGGTACCTTGTCGATGTTCCTTTa                                                                           <  2:2164069/149‑1 (MQ=255)
                                                                            tACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGGTGTTCGGTACCGTGTGGATGTTGCGTTaac                                                                         >  1:780332/1‑147 (MQ=255)
                                                                             aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAccc                                                                        >  2:9521952/1‑149 (MQ=255)
                                                                             aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCGGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAccc                                                                        >  1:5086077/1‑149 (MQ=255)
                                                                               cGCCTTGCTGACCGCCGCGGGCCGCGCCATCTGGGGCGCCGCGGCGGTGCTGGATTTCGACTCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                      <  1:2310054/149‑1 (MQ=255)
                                                                               cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                      <  2:3508353/149‑1 (MQ=255)
                                                                                gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                     >  2:2822984/1‑149 (MQ=255)
                                                                                gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCTTGTCGTTGTTGCTTTAACCACt                                                                     >  2:4074485/1‑149 (MQ=255)
                                                                                 ccTTGCTGTCCGCCGCCGGCAGCGCACTCTGGGGCGCCGCGGCGGTGTTGTCATTCGTATCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAATCGTAGTGCTGTTCGGTACCTTGTCTATGTTCCTTTACCCACTg                                                                    <  1:10468048/149‑1 (MQ=255)
                                                                                   ttGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCAAAAAGCGCCCGGGCAGGCGGAAGCGGAGTGCTGGTGGGGACCATGTCGATGTTGCTTTAACCACAGGc                                                                  >  1:1441953/1‑149 (MQ=255)
                                                                                   ttGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAGCGCCCTGGCAGTCGGAAGCGTCGTGCTGTTCGGTACCGTGTCGATGTTACATTACCCACTGGc                                                                  >  2:3036190/1‑149 (MQ=255)
                                                                                   ttGCTGACCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                                  <  2:6166172/149‑1 (MQ=255)
                                                                                    tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCg                                                                 <  2:5533626/149‑1 (MQ=255)
                                                                                     gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGa                                                                <  2:9329353/149‑1 (MQ=255)
                                                                                      cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCACAAAAGCGCCAGGGCAGTCGGAAGCGGAGTGCTGGTCGGTGCCTTGTCGATGTTCCGTTACCCACGGGCgcg                                                               >  1:4272459/1‑147 (MQ=255)
                                                                                      cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAACCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                               <  2:10288051/149‑1 (MQ=255)
                                                                                        gACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCa                                                             <  2:3125550/149‑1 (MQ=255)
                                                                                        gACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCGGGCATTCGCATCGGCGCTGCGCAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGGGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGATCa                                                             <  2:9455951/149‑1 (MQ=255)
                                                                                          ccgcccccgCCAGCCCCATCGGCGGCCCCGCCGCCGTGCTGGAATTCCCACCGGCGCCGCCCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                           <  1:6773374/149‑1 (MQ=255)
                                                                                           cgccgccgGCCGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGATTCGGCGATGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAACTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                          <  2:6536092/149‑1 (MQ=255)
                                                                                            gccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCGACAAGCGCCAAGGCGGTCGGGAGCGGAGTGCTGGGCGGGACCAAGGCGATGTGCCGGTACCCCCAGGCCGGCAATGc                                                         >  1:8659007/1‑149 (MQ=255)
                                                                                              cgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCTTTTGGATGTTGCTTTAACCACAGGCGATAAATGCCg                                                       >  1:10250878/1‑149 (MQ=255)
                                                                                               gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCCCATCGCCCCCCATAAAAGCCCCATGGCACTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGTTCCTTTACCCAAAGGCGATAAAATCCgg                                                      >  2:3789270/1‑149 (MQ=255)
                                                                                                 cgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCTAAAAGCGCCCAGGCCGTCGGAAGCGTAGTGCTGTTCGGTAACTTGTCGATGTTCCATTACCCACAGGCGAAAAAAGCCggct                                                    >  1:9329131/1‑149 (MQ=255)
                                                                                                  ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATGGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctg                                                   <  1:527545/149‑1 (MQ=255)
                                                                                                     aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCggctggct                                                >  1:6088499/1‑149 (MQ=255)
                                                                                                       cgcCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGc                                              >  2:5831287/1‑149 (MQ=255)
                                                                                                          cctCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATc                                           <  2:7297257/147‑1 (MQ=255)
                                                                                                          cATCTGGGGGGCCGCCGCCGGGCTGGCATGCGACGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATc                                           <  1:10132579/149‑1 (MQ=255)
                                                                                                           aTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                          >  2:5783326/1‑149 (MQ=255)
                                                                                                               gcggcgCCGCCGCGGTGGTGGAATTCTACGCGGGGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGAt                                      <  1:976011/149‑1 (MQ=255)
                                                                                                                 ggcggcGCCGGGGTGCTGGCTTTCGACTGGGGGCGGCGCCGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAc                                    <  1:8239221/149‑1 (MQ=255)
                                                                                                                    gccgccgccGGGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGGCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                                 <  2:161313/149‑1 (MQ=255)
                                                                                                                     ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATCAATGCCGGCGGGCTGCATCTGGATAACaag                                >  1:4115230/1‑147 (MQ=255)
                                                                                                                      cgccgccgTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCGTGTCGATGTTCCGTTAACCACAGGCGATCAATGCCCGCTGGCTGCATCTGGATACCATgg                               >  1:4243209/1‑149 (MQ=255)
                                                                                                                         cgcggTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgcg                            <  1:5560810/149‑1 (MQ=255)
                                                                                                                         cgccgTGCTGGCATCCGAAGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgcg                            <  2:9016131/149‑1 (MQ=255)
                                                                                                                             gTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                        >  2:7415767/1‑149 (MQ=255)
                                                                                                                               gCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgc                      <  2:9773212/149‑1 (MQ=255)
                                                                                                                                  ggatttCGAATCGGCGCTGCGCTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                   <  2:7184674/145‑1 (MQ=255)
                                                                                                                                  ggCATCCGAATCGGCGCTGCGCAGCGCCCCCCATAACCGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                   <  2:7028629/149‑1 (MQ=255)
                                                                                                                                    cATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATCAATTCCGGCTGGCTGCATCTGGATAACAAGGGCGCCCGgctgctg                 >  2:8758726/1‑149 (MQ=255)
                                                                                                                                    cATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAAGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctg                 >  2:8072509/1‑149 (MQ=255)
                                                                                                                                    cATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTAGATACCATGGGCGCCGGgctgctg                 >  2:2642748/1‑149 (MQ=255)
                                                                                                                                      ttCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct               <  2:2493812/149‑1 (MQ=255)
                                                                                                                                      tgcgacgcggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct               <  1:6142115/142‑1 (MQ=255)
                                                                                                                                       tctaaTGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCTGCTGGCTGCATCTGGATACCATGGGCGCCGTGctgctgctg              <  2:1525268/146‑1 (MQ=255)
                                                                                                                                           cgcGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg          <  1:3152316/147‑1 (MQ=255)
                                                                                                                                             cGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAGGGGCGCCGGGCTGCTGCTGggcggc        >  1:1569061/1‑149 (MQ=255)
                                                                                                                                             cGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcggc        >  1:1568638/1‑149 (MQ=255)
                                                                                                                                                cgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt     >  2:3366398/1‑149 (MQ=255)
                                                                                                                                                cgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt     >  1:6403015/1‑149 (MQ=255)
                                                                                                                                                  cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt   <  2:9858548/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTC  >  NC_002947/4160285‑4160580

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