Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,176 G→T 16.6% P471P (CCG→CCT rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,1760GT16.6% 523.6 / 78.5 196P471P (CCG→CCTrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (101/62);  new base T (19/14);  total (120/76)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

TGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAA  >  NC_002947/543028‑543325
                                                                                                                                                    |                                                                                                                                                     
tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                                                                                                      >  1:1645533/1‑150 (MQ=255)
tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                                                                                                      >  2:3958509/1‑150 (MQ=255)
  cTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGACTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTg                                                                                                                                                    >  1:5017721/1‑150 (MQ=255)
   tgaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgt                                                                                                                                                   >  2:1155098/3‑149 (MQ=255)
   tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGATCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                                                                                                                                   <  2:4216121/150‑1 (MQ=255)
        cGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTg                                                                                                                                              <  2:452896/150‑1 (MQ=255)
        cGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTg                                                                                                                                              <  1:2227588/150‑1 (MQ=255)
          tGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGa                                                                                                                                            >  1:2848313/1‑150 (MQ=255)
           gCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGaa                                                                                                                                           <  2:2209115/150‑1 (MQ=255)
           gCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGaa                                                                                                                                           <  1:3878836/150‑1 (MQ=255)
             aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                                                                                                         >  2:1784759/1‑150 (MQ=255)
              ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTGCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                                                                                                        >  1:2712415/1‑150 (MQ=255)
               cGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                                                                                                                                       >  1:4899210/1‑150 (MQ=255)
                 aCCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAg                                                                                                                                     <  2:2878118/150‑1 (MQ=255)
                  cgctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                                                                                                    >  2:3340049/6‑150 (MQ=255)
                   gctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCt                                                                                                                                   <  1:4708438/146‑1 (MQ=255)
                   gctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCt                                                                                                                                   >  2:2949635/5‑150 (MQ=255)
                   ccTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCt                                                                                                                                   >  1:3433725/1‑150 (MQ=255)
                    cTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTc                                                                                                                                  >  2:1390206/1‑150 (MQ=255)
                     ttCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAATGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCg                                                                                                                                 >  2:4342544/1‑150 (MQ=255)
                       cACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGaa                                                                                                                               <  1:4858093/150‑1 (MQ=255)
                        aCCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGAGTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCTGCTTGAAg                                                                                                                              >  2:1764194/1‑150 (MQ=255)
                         ccGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgc                                                                                                                             <  2:2116320/150‑1 (MQ=255)
                           gTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAgcac                                                                                                                           <  1:283836/150‑3 (MQ=255)
                             cttggTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGcacgc                                                                                                                         <  1:588985/147‑2 (MQ=255)
                             cTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGt                                                                                                                         <  1:715388/150‑1 (MQ=255)
                              tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                        >  2:2313911/1‑150 (MQ=255)
                              tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                        >  2:4279528/1‑150 (MQ=255)
                              tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                        >  2:3077956/1‑150 (MQ=255)
                               tggTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGcacgcgt                                                                                                                       >  2:848516/2‑147 (MQ=255)
                               gggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                                                                                                       >  2:741244/1‑150 (MQ=255)
                                 gTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcgc                                                                                                                     <  1:925609/150‑1 (MQ=255)
                                     ccAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgat                                                                                                                 >  1:4207182/1‑150 (MQ=255)
                                       aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga                                                                                                               >  2:5052315/1‑150 (MQ=255)
                                          cGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGt                                                                                                            <  1:585077/150‑1 (MQ=255)
                                           gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                                                                                           <  2:712477/150‑1 (MQ=255)
                                           gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                                                                                           <  1:2836770/150‑1 (MQ=255)
                                            tttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCg                                                                                                          >  1:2361266/1‑150 (MQ=255)
                                              tGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGAc                                                                                                        <  2:3213298/150‑1 (MQ=255)
                                               gAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa                                                                                                         <  1:1265795/148‑1 (MQ=255)
                                                   ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                                                                                                   >  1:2679137/1‑150 (MQ=255)
                                                   ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                                                                                                   >  1:2063334/1‑150 (MQ=255)
                                                           aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTg                                                                                           >  1:2327557/1‑150 (MQ=255)
                                                           aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTg                                                                                           >  2:4912077/1‑150 (MQ=255)
                                                            tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt                                                                                          >  2:607743/1‑150 (MQ=255)
                                                            tCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTgc                                                                                          <  1:3516083/150‑2 (MQ=255)
                                                             cGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTc                                                                                         <  2:1288153/150‑1 (MQ=255)
                                                               aaggcaaggcCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTgcacc                                                                                       >  2:2143742/1‑146 (MQ=255)
                                                                  gTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccag                                                                                    >  1:3775166/1‑150 (MQ=255)
                                                                  gTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccag                                                                                    >  1:3174095/1‑150 (MQ=255)
                                                                    aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                  >  1:1856804/1‑150 (MQ=255)
                                                                       gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                               >  1:2571696/1‑150 (MQ=255)
                                                                       gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                               >  2:2428392/1‑150 (MQ=255)
                                                                         tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                             >  1:2574624/1‑150 (MQ=255)
                                                                                  gCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAat                                                                    <  2:4426338/150‑1 (MQ=255)
                                                                                      ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                >  1:4218784/1‑150 (MQ=255)
                                                                                        gCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGt                                                              >  2:266346/1‑150 (MQ=255)
                                                                                               tgtgCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCg                                                       >  1:2024307/1‑150 (MQ=255)
                                                                                                gtgCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGc                                                      <  1:1445271/150‑1 (MQ=255)
                                                                                                   cagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAgg                                                   >  2:516437/3‑150 (MQ=255)
                                                                                                    tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                  >  1:4466464/1‑150 (MQ=255)
                                                                                                    tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                  >  1:1058247/1‑150 (MQ=255)
                                                                                                     gcGTACAACGCCGACTTCGACGGTGACAAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGAc                                                 >  2:4206142/1‑150 (MQ=255)
                                                                                                      cGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAca                                                >  2:1089992/1‑149 (MQ=255)
                                                                                                       gTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGt                                               >  2:1280762/1‑150 (MQ=255)
                                                                                                        tACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGtt                                              <  1:3609758/150‑1 (MQ=255)
                                                                                                         aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                             <  1:3809817/150‑1 (MQ=255)
                                                                                                          cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                                                                                                       >  2:3258469/1‑75 (MQ=255)
                                                                                                          cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                                                                                                       <  1:3258469/75‑1 (MQ=255)
                                                                                                          cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGt                                            <  1:3401271/150‑1 (MQ=255)
                                                                                                           aaCGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAcatggtt                                           >  2:2246800/1‑144 (MQ=255)
                                                                                                           aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                           >  1:139435/1‑150 (MQ=255)
                                                                                                           aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                           >  2:1863320/1‑150 (MQ=255)
                                                                                                            aCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTAc                                          >  2:2024459/1‑150 (MQ=255)
                                                                                                             cGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACt                                         >  2:1550789/1‑150 (MQ=255)
                                                                                                              gCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTg                                        <  2:3019393/150‑1 (MQ=255)
                                                                                                               ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                       >  1:4214784/1‑150 (MQ=255)
                                                                                                               ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                       <  1:551267/150‑1 (MQ=255)
                                                                                                               ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                       >  2:2626019/1‑150 (MQ=255)
                                                                                                                cGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTggg                                      >  1:3804275/1‑150 (MQ=255)
                                                                                                                 gACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGt                                     >  1:1603434/1‑150 (MQ=255)
                                                                                                                    ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGGTGTACTGGGTCTg                                  >  1:919092/1‑150 (MQ=255)
                                                                                                                      cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                >  2:3667521/1‑150 (MQ=255)
                                                                                                                      cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                <  2:373442/150‑1 (MQ=255)
                                                                                                                       gACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtac                               >  1:746914/1‑150 (MQ=255)
                                                                                                                         cGGTGACCAGATGGCCGTTCAGGTGCCGCTGAGCCTGGAAGCCCAGCTCGAAGGGCGTGCGCTGATGATGTCGACCAACTTCATCCTGTCGCCAGCCAACGGTGAGCCTATCATCGTTCCGTCGCAGGTCGTTGTACTGGGTCTGtacta                             <  2:4435258/150‑1 (MQ=255)
                                                                                                                           gTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACa                           <  1:4139624/150‑1 (MQ=255)
                                                                                                                                    tGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCggg                  >  2:337646/1‑150 (MQ=255)
                                                                                                                                      gCCTTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTTAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                <  2:4817673/150‑1 (MQ=255)
                                                                                                                                      gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                <  1:4030775/150‑1 (MQ=255)
                                                                                                                                       ccGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAg               >  1:3892848/1‑150 (MQ=255)
                                                                                                                                        cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc              >  1:3442844/1‑150 (MQ=255)
                                                                                                                                        cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc              >  2:480946/1‑150 (MQ=255)
                                                                                                                                        cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc              >  2:5074572/1‑150 (MQ=255)
                                                                                                                                             aCGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCa         >  2:2262037/1‑150 (MQ=255)
                                                                                                                                               gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt       <  2:2712415/150‑1 (MQ=255)
                                                                                                                                               gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt       <  1:2097508/150‑1 (MQ=255)
                                                                                                                                                tGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATg      <  1:3922006/150‑1 (MQ=255)
                                                                                                                                                 gCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGc     <  2:4714279/150‑1 (MQ=255)
                                                                                                                                                    gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa  <  2:3433725/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAA  >  NC_002947/543028‑543325

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