Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,245 C→A 5.8% A494A (GCC→GCA 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,2450CA 5.8% 577.2 / 17.6 172A494A (GCC→GCArpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (89/73);  new base A (7/3);  total (96/76)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 186 total aligned reads displayed.

AAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGCCCT  >  NC_002947/543096‑543394
                                                                                                                                                     |                                                                                                                                                     
aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                                                       >  1:2341741/1‑150 (MQ=255)
aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                                                       >  2:2792034/1‑150 (MQ=255)
   gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                                                                    >  1:2571696/1‑150 (MQ=255)
   gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                                                                    >  2:2428392/1‑150 (MQ=255)
   gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                                                                    >  1:3254340/1‑150 (MQ=255)
     tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                                                                                                  >  1:2574624/1‑150 (MQ=255)
              gCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAat                                                                                                                                         <  2:4426338/150‑1 (MQ=255)
               cACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatc                                                                                                                                        <  2:3254340/150‑1 (MQ=255)
                aCCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatca                                                                                                                                       <  2:4899210/150‑1 (MQ=255)
                  ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                                                                     <  1:110285/150‑1 (MQ=255)
                  ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                                                                     >  1:4218743/1‑150 (MQ=255)
                  ccGCTCGTCTGTACAGAGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGatcatc                                                                                                                                     <  2:4624354/150‑1 (MQ=255)
                    gCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGt                                                                                                                                   >  2:1990173/1‑150 (MQ=255)
                        gTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCg                                                                                                                               >  1:2810719/1‑150 (MQ=255)
                           tgtgCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCg                                                                                                                            >  1:2024307/1‑150 (MQ=255)
                            gtgCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGc                                                                                                                           <  1:1445271/150‑1 (MQ=255)
                            gtacagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGc                                                                                                                           <  1:2706131/145‑1 (MQ=255)
                                tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                                                                                       >  1:1058247/1‑150 (MQ=255)
                                   gTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGt                                                                                                                    >  2:1280884/1‑150 (MQ=255)
                                   gTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGt                                                                                                                    >  2:1280762/1‑150 (MQ=255)
                                     tcAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAcatgg                                                                                                                  <  1:2299223/149‑5 (MQ=255)
                                     aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                                                                                  >  2:4439787/1‑150 (MQ=255)
                                      ccaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACACCCTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGt                                                                                                                 <  1:3401393/148‑1 (MQ=255)
                                       aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                                                                                                >  1:139435/1‑150 (MQ=255)
                                       aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                                                                                                >  1:4144696/1‑150 (MQ=255)
                                        aCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTAc                                                                                                               >  2:2024459/1‑150 (MQ=255)
                                          gCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTg                                                                                                             <  2:4374418/150‑1 (MQ=255)
                                          gCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTg                                                                                                             <  2:3019393/150‑1 (MQ=255)
                                           ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                                                                                            >  2:2626019/1‑150 (MQ=255)
                                            cGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTggg                                                                                                           >  1:3804275/1‑150 (MQ=255)
                                            cGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTggg                                                                                                           <  1:1627235/150‑1 (MQ=255)
                                             gACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGt                                                                                                          >  1:1603434/1‑150 (MQ=255)
                                                ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                                                                                       >  1:919176/1‑150 (MQ=255)
                                                  cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                                                                                     >  2:5113081/1‑150 (MQ=255)
                                                  cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                                                                                     >  2:3667521/1‑150 (MQ=255)
                                                   gACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtac                                                                                                    >  1:746822/1‑150 (MQ=255)
                                                     cGGTGACCAGATGGCCGTTCAGGTGCCGCTGAGCCTGGAAGCCCAGCTCGAAGGGCGTGCGCTGATGATGTCGACCAACTTCATCCTGTCGCCAGCCAACGGTGAGCCTATCATCGTTCCGTCGCAGGTCGTTGTACTGGGTCTGtacta                                                                                                  <  2:4435258/150‑1 (MQ=255)
                                                     cGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtacta                                                                                                  <  1:759732/150‑1 (MQ=255)
                                                       gTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACa                                                                                                >  2:1389247/1‑150 (MQ=255)
                                                              gATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCg                                                                                         >  1:670344/1‑150 (MQ=255)
                                                                tGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCggg                                                                                       >  2:337646/1‑150 (MQ=255)
                                                                  gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                     <  1:3510644/150‑1 (MQ=255)
                                                                  gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                     <  2:4817445/150‑1 (MQ=255)
                                                                  gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                     <  1:4030775/150‑1 (MQ=255)
                                                                   ccGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAg                                                                                    >  1:3892848/1‑150 (MQ=255)
                                                                    cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                   <  1:514686/150‑1 (MQ=255)
                                                                    cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                   >  1:3442844/1‑150 (MQ=255)
                                                                    cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                   >  2:480946/1‑150 (MQ=255)
                                                                       tCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCAt                                                                                <  2:540368/150‑1 (MQ=255)
                                                                         aCGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCa                                                                              >  2:2262037/1‑150 (MQ=255)
                                                                           gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt                                                                            <  1:2097508/150‑1 (MQ=255)
                                                                           gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt                                                                            <  2:2712415/150‑1 (MQ=255)
                                                                           gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt                                                                            <  2:4207182/150‑1 (MQ=255)
                                                                            tGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATg                                                                           <  1:3922006/150‑1 (MQ=255)
                                                                                gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa                                                                       <  2:3433725/150‑1 (MQ=255)
                                                                                gCTGACCCTGGAAGCCCAGCTCGAAGCGCGAGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAACCCATCAATGCCaa                                                                       <  2:2216256/150‑1 (MQ=255)
                                                                                    aCCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                   <  1:3494330/134‑1 (MQ=255)
                                                                                      ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                                                                 >  1:2253194/1‑150 (MQ=255)
                                                                                      ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                                                                 >  1:2170707/1‑150 (MQ=255)
                                                                                       cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                                <  2:3800808/150‑1 (MQ=255)
                                                                                       cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                                <  1:1962522/150‑1 (MQ=255)
                                                                                       cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                                <  1:888464/150‑1 (MQ=255)
                                                                                        tGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAg                                                               >  2:2232723/1‑150 (MQ=255)
                                                                                        tGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAg                                                               <  2:1645533/150‑1 (MQ=255)
                                                                                               ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                        >  2:4230234/1‑150 (MQ=255)
                                                                                               ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                        >  1:2161408/1‑150 (MQ=255)
                                                                                               ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                        >  1:461444/1‑150 (MQ=255)
                                                                                               ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                        <  1:4798558/150‑1 (MQ=255)
                                                                                               ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                        <  2:310624/150‑1 (MQ=255)
                                                                                                cAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCCTCAATGCCAAGGGCGAAGGTCGCgtg                                                       >  2:4689165/1‑150 (MQ=255)
                                                                                                cAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgtg                                                       >  2:4887576/1‑150 (MQ=255)
                                                                                                    tCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCg                                                   >  2:4976013/1‑150 (MQ=255)
                                                                                                    tCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCg                                                   >  2:4201431/1‑150 (MQ=255)
                                                                                                        aGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGa                                               <  2:4780774/150‑1 (MQ=255)
                                                                                                              tGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCa                                         <  1:5052324/150‑1 (MQ=255)
                                                                                                              tGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCa                                         <  1:5052315/150‑1 (MQ=255)
                                                                                                                  cTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGaa                                     <  2:2574624/150‑1 (MQ=255)
                                                                                                                  cTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGaa                                     <  1:1784759/150‑1 (MQ=255)
                                                                                                                  cTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGaa                                     <  1:3854907/150‑1 (MQ=255)
                                                                                                                   tgatgaTGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAg                                    >  1:3032025/1‑150 (MQ=255)
                                                                                                                     atgatgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCCATCATCGATCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGTAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCAGCAGGAAGTc                                  >  1:4199753/1‑150 (MQ=255)
                                                                                                                     atgatgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTc                                  <  2:2361266/150‑1 (MQ=255)
                                                                                                                         tgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGAcc                              >  1:2484382/1‑150 (MQ=255)
                                                                                                                          gTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCAACcg                             >  2:844941/1‑150 (MQ=255)
                                                                                                                              aCCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTa                         <  2:2508502/150‑1 (MQ=255)
                                                                                                                               ccAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAt                        <  2:2679137/150‑1 (MQ=255)
                                                                                                                                caacaaCATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAtt                       >  2:4480135/1‑150 (MQ=255)
                                                                                                                                caacaaCATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAtt                       >  1:441560/1‑150 (MQ=255)
                                                                                                                                  acaacaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTcc                     <  1:4279528/150‑1 (MQ=255)
                                                                                                                                  acaacaTCCTGTCGCCAGACAACGGTGAGCCAATCATCGTTCCGTCGCAGTACGTTGTACTGTGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGTCGAAGGTAGCGTGTTCGCCGACCTGCAGGAAGTAGACCGCGTATTcc                     <  1:2888204/150‑1 (MQ=255)
                                                                                                                                   caacaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcg                    >  2:3713302/1‑150 (MQ=255)
                                                                                                                                    aacaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgc                   <  2:3855154/150‑1 (MQ=255)
                                                                                                                                    aacaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgc                   >  1:1645645/1‑150 (MQ=255)
                                                                                                                                     acaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgcg                  <  2:2571770/150‑1 (MQ=255)
                                                                                                                                        tCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCg               <  2:3775166/150‑1 (MQ=255)
                                                                                                                                        tCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCg               <  2:1856804/150‑1 (MQ=255)
                                                                                                                                                    ccaACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGccc   >  2:4573377/1‑150 (MQ=255)
                                                                                                                                                     caACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGCCCt  <  1:1550789/150‑1 (MQ=255)
                                                                                                                                                     caACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGCCCt  <  2:2327557/150‑1 (MQ=255)
                                                                                                                                                     caACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGCCCt  >  1:12143/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGCCCT  >  NC_002947/543096‑543394

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: