Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 545,513 G→T 6.1% K1251N (AAG→AAT)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947545,5130GT 6.1% 234.2 / 2.7 76K1251N (AAG→AAT) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base T (2/3);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTAAACCGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCATCTC  >  NC_002947/545367‑545659
                                                                                                                                                   |                                                                                                                                                  
cAGGTAAACCGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGca                                                                                                                                                  <  3:101426/149‑1 (MQ=255)
   gTAAACCGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACAt                                                                                                                                               >  4:207809/1‑149 (MQ=255)
      aaCCGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATc                                                                                                                                              >  6:19667/1‑147 (MQ=255)
        ccGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCgaga                                                                                                                                          >  6:334731/1‑149 (MQ=255)
         cgcgGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTTCGTCTGCTGGGTGTGTTTGCGCTTGCGAAGTACATCGTCAACGAGATTCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCgaga                                                                                                                                          <  1:299361/149‑1 (MQ=255)
            ggTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCTCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGTCGTGAAGATCAACGACAAGCACATCGAGACCAt                                                                                                                                      >  4:92788/1‑149 (MQ=255)
               gAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCt                                                                                                                                   >  3:127510/1‑149 (MQ=255)
                 aGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGc                                                                                                                                 <  1:31033/149‑1 (MQ=255)
                 aGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGc                                                                                                                                 <  5:19667/149‑1 (MQ=255)
                   ttATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGt                                                                                                                               >  4:281823/1‑149 (MQ=255)
                     ctctcCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTc                                                                                                                              <  6:118819/147‑1 (MQ=255)
                     aTCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGt                                                                                                                               >  8:236678/1‑147 (MQ=255)
                     aTCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGCCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGt                                                                                                                               >  3:148788/1‑147 (MQ=255)
                       ctcCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGt                                                                                                                               >  7:219599/1‑145 (MQ=255)
                        tcCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAt                                                                                                                          >  5:333576/1‑149 (MQ=255)
                            aCGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAtgctg                                                                                                                      <  2:257116/149‑1 (MQ=255)
                             cGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAtgctgc                                                                                                                     <  4:284631/149‑1 (MQ=255)
                             cGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAtgctgc                                                                                                                     <  8:217611/149‑1 (MQ=255)
                                cccGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTgcgc                                                                                                                   >  7:181218/1‑148 (MQ=255)
                                 ccGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTgcgc                                                                                                                   >  4:63482/1‑147 (MQ=255)
                                   gAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAgg                                                                                                               <  3:208598/149‑1 (MQ=255)
                                     gCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGtt                                                                                                             >  7:194345/1‑149 (MQ=255)
                                              aTGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTcc                                                                                                    <  3:207809/149‑1 (MQ=255)
                                                gACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAATCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTcc                                                                                                    >  6:159721/1‑147 (MQ=255)
                                                gACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAATCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTcc                                                                                                    >  6:165715/1‑147 (MQ=255)
                                                    tCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACACCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCg                                                                                              >  2:220238/1‑149 (MQ=255)
                                                      cTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGTCCGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCggg                                                                                            <  5:34491/149‑1 (MQ=255)
                                                      cTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCggg                                                                                            <  8:219599/149‑1 (MQ=255)
                                                               ctgGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTc                                                                                      <  8:194345/146‑1 (MQ=255)
                                                                 gggTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGTTTCCAGCt                                                                                 >  2:280330/1‑149 (MQ=255)
                                                                  ggTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCt                                                                                 <  4:127510/148‑1 (MQ=255)
                                                                         gTGCGCTGGCGAAGTACATCGTAAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCaa                                                                          <  5:334731/148‑1 (MQ=255)
                                                                                 gCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGAc                                                                   >  1:28849/1‑147 (MQ=255)
                                                                                    aaGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGAt                                                              >  3:246692/1‑149 (MQ=255)
                                                                                      gTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAtg                                                                                                                         >  1:161341/1‑88 (MQ=255)
                                                                                      gTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGAtg                                                                                                                         <  2:161341/88‑1 (MQ=255)
                                                                                      gTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATg                                                             >  7:131318/1‑148 (MQ=255)
                                                                                        aCATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCCACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGGTGAGAGCTCCGAGTCGGGCGCTTCCAGCTTCATCAAGGGCGACCAGATGGaa                                                          >  5:117007/1‑149 (MQ=255)
                                                                                         cATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGaa                                                          <  8:131318/148‑1 (MQ=255)
                                                                                          aTCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACt                                                        >  3:351688/1‑149 (MQ=255)
                                                                                                aaCGAGA‑TCCAGGCCGTTTACCGTCTGCAAGGCGTGAAGATCTACGAGAAGCACATCGAGACCATTCTGAGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTAGGGCGATTCCAGCTTCCTCAAGGGCGCCCAGATGGAGCTGAcc                                                    >  2:298917/1‑147 (MQ=255)
                                                                                                  cGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCCgg                                                >  4:166191/1‑149 (MQ=255)
                                                                                                     ga‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGc                                             >  3:232971/1‑149 (MQ=255)
                                                                                                          cAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGc                                             <  6:315224/145‑1 (MQ=255)
                                                                                                          cAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGc                                             >  5:315224/1‑145 (MQ=255)
                                                                                                          cAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACACCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTgg                                          <  1:220238/148‑1 (MQ=255)
                                                                                                           aGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCggg                                                                                            <  4:283515/97‑1 (MQ=255)
                                                                                                               cGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAATCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAgaga                                    <  5:159721/149‑1 (MQ=255)
                                                                                                               cGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAATCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAgaga                                    <  5:165715/149‑1 (MQ=255)
                                                                                                                 tttACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAAc                                  <  3:281823/149‑1 (MQ=255)
                                                                                                                     ccGTCTGCAAGGCGTGAAGATCAACGACAATCACATGGAGACGTTTCTGCGTCAGATGCTGCGCAAGTTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGc                              <  2:28849/149‑1 (MQ=255)
                                                                                                                     ccGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGc                              >  1:350060/1‑149 (MQ=255)
                                                                                                                      cGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGAGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCg                             >  5:196709/1‑149 (MQ=255)
                                                                                                                        tCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGtc                           >  8:220820/1‑149 (MQ=255)
                                                                                                                        tCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGAGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGtc                           >  3:146507/1‑149 (MQ=255)
                                                                                                                           gCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCg                        >  1:337802/1‑149 (MQ=255)
                                                                                                                            cAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGtctc                         <  4:148788/147‑1 (MQ=255)
                                                                                                                            cAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGtct                          <  4:246692/146‑1 (MQ=255)
                                                                                                                             aaGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGc                       >  4:175873/1‑148 (MQ=255)
                                                                                                                              aGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTa                     <  2:337802/149‑1 (MQ=255)
                                                                                                                                gCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGAGCAGATGGAACTGACCCAGGTCCTGGTAGAGAACGAGCGTCtagctagc                   >  4:61739/1‑149 (MQ=255)
                                                                                                                                gCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGGCGGGCGATTCCCGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACCAGCGTCGCGCTAGc                   <  8:181218/149‑1 (MQ=255)
                                                                                                                                  gTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAAAGAGCGTCTCGCTAGc                   >  3:343168/1‑147 (MQ=255)
                                                                                                                                  cagacgATCAACGCCAAGCACCTCCCGACCCTTCTGCGGCACCTGCTGCGCAACGGTGACATCTCCGCCTCCGGCGAATCCCCCTTCATCAAGGGCGACCAGATGGAACGGACCCCGGGGCGGGGAGAGAACGAGCGGCTCGCTAGCGa                 <  8:297317/144‑1 (MQ=255)
                                                                                                                                    gAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAg                >  3:276588/1‑148 (MQ=255)
                                                                                                                                       gATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGc                              >  8:12494/1‑131 (MQ=255)
                                                                                                                                       gATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGc                              <  7:12494/131‑1 (MQ=255)
                                                                                                                                       gATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGAc             >  4:123380/1‑148 (MQ=255)
                                                                                                                                         tCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAg          >  3:178044/1‑149 (MQ=255)
                                                                                                                                            aCGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGtt        >  3:24026/1‑148 (MQ=255)
                                                                                                                                             cGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCa      >  3:231886/1‑149 (MQ=255)
                                                                                                                                             cGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCa      >  3:226934/1‑149 (MQ=255)
                                                                                                                                              gACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCAt     <  3:63482/149‑1 (MQ=255)
                                                                                                                                              gACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCAt     <  4:351688/149‑1 (MQ=255)
                                                                                                                                               aCAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCAtc    <  3:175873/149‑1 (MQ=255)
                                                                                                                                                 aaGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCAtctc  >  7:87873/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGTAAACCGCGGTGAAGTTATCTCCGACGGCCCGAGCGATCCGCATGACATCCTGCGTCTGCTGGGTGTGAGTGCGCTGGCGAAGTACATCGTCAACGAGA‑TCCAGGACGTTTACCGTCTGCAAGGCGTGAAGATCAACGACAAGCACATCGAGACCATTCTGCGTCAGATGCTGCGCAAGGTTGAGATCTCCGAGTCGGGCGATTCCAGCTTCATCAAGGGCGACCAGATGGAACTGACCCAGGTGCTGGTAGAGAACGAGCGTCTCGCTAGCGAGGACAAGTTCATCTC  >  NC_002947/545367‑545659

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