Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,242 A→G 5.8% P493P (CCA→CCG 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,2420AG 5.8% 548.9 / 14.1 172P493P (CCA→CCGrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (89/73);  new base G (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.44e-01

Only 100 of 182 total aligned reads displayed.

GTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGC  >  NC_002947/543094‑543391
                                                                                                                                                    |                                                                                                                                                     
gTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccag                                                                                                                                                      >  1:3775166/1‑150 (MQ=255)
  aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                                                    >  1:2341741/1‑150 (MQ=255)
     gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                                                                 >  1:3254340/1‑150 (MQ=255)
     gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                                                                 >  1:2571770/1‑150 (MQ=255)
       tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                                                                                               >  1:3903502/1‑150 (MQ=255)
       tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                                                                                               >  2:3037304/1‑150 (MQ=255)
        aTTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGt                                                                                                                                              >  1:3887497/1‑150 (MQ=255)
            aGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGc                                                                                                                                          >  1:3019393/1‑150 (MQ=255)
                gCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAat                                                                                                                                      <  2:4426338/150‑1 (MQ=255)
                  aCCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatca                                                                                                                                    <  2:4899210/150‑1 (MQ=255)
                    ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                                                                  <  1:110285/150‑1 (MQ=255)
                    ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                                                                  >  1:4218784/1‑150 (MQ=255)
                    ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                                                                  >  1:2508502/1‑150 (MQ=255)
                      gCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGt                                                                                                                                >  2:266346/1‑150 (MQ=255)
                      gCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGt                                                                                                                                >  2:1990173/1‑150 (MQ=255)
                          gTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCg                                                                                                                            >  1:2810719/1‑150 (MQ=255)
                             tgtgCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCg                                                                                                                         >  1:2024307/1‑150 (MQ=255)
                             tgtacagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCg                                                                                                                         >  2:2706131/7‑150 (MQ=255)
                               tgctgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCa                                                                                                                       <  2:3393157/150‑1 (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:1280884/1‑150 (MQ=255)
                                      tACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGtt                                                                                                                <  1:1440/150‑1 (MQ=255)
                                       tcAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAcatgg                                                                                                               <  1:2299223/149‑5 (MQ=255)
                                       aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                                                                               <  1:3809817/150‑1 (MQ=255)
                                        ccaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACACCCTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGt                                                                                                              <  1:3401393/148‑1 (MQ=255)
                                         aaCGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAcatggtt                                                                                                             >  2:2246800/1‑144 (MQ=255)
                                         aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                                                                                             >  2:1863320/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)
                                             ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                                                                                         >  1:4214784/1‑150 (MQ=255)
                                               gACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGt                                                                                                       <  2:3687502/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)
                                                     gACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtac                                                                                                 >  2:3349426/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)
                                                                    gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                  <  1:4030775/150‑1 (MQ=255)
                                                                    gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                  <  1:3510644/150‑1 (MQ=255)
                                                                    gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                                                                  <  2:4817445/150‑1 (MQ=255)
                                                                      cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                <  1:514686/150‑1 (MQ=255)
                                                                      cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                >  1:3442844/1‑150 (MQ=255)
                                                                         tCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCAt                                                                             <  2:540368/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:4098232/150‑1 (MQ=255)
                                                                                  gCTGACCCTGGAAGCCCAGCTCGAAGCGCGAGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAACCCATCAATGCCaa                                                                    <  2:2216256/150‑1 (MQ=255)
                                                                                     gACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAggg                                                                 <  2:421113/150‑1 (MQ=255)
                                                                                      aCCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                                                                >  2:3494330/1‑134 (MQ=255)
                                                                                        ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                                                              >  2:4798558/1‑150 (MQ=255)
                                                                                        ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                                                              >  1:2253194/1‑150 (MQ=255)
                                                                                         cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                             <  2:3800808/150‑1 (MQ=255)
                                                                                         cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                             <  1:1962522/150‑1 (MQ=255)
                                                                                         cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                             <  1:2167824/150‑1 (MQ=255)
                                                                                         cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                                                             <  1:888464/150‑1 (MQ=255)
                                                                                          tGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAg                                                            <  2:1645533/150‑1 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     >  2:4230234/1‑150 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     >  1:461444/1‑150 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     >  1:1625353/1‑150 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     <  1:4798558/150‑1 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     >  1:1526783/1‑150 (MQ=255)
                                                                                                 ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                                                     <  2:310624/150‑1 (MQ=255)
                                                                                                  cAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgtg                                                    >  2:4887576/1‑150 (MQ=255)
                                                                                                    gCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGtt                                                  >  1:2097693/1‑150 (MQ=255)
                                                                                                      tCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCg                                                >  2:4976013/1‑150 (MQ=255)
                                                                                                        gAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGcc                                              >  1:2680463/1‑150 (MQ=255)
                                                                                                                tGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCa                                      <  1:5052315/150‑1 (MQ=255)
                                                                                                                    cTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGaa                                  <  2:2574624/150‑1 (MQ=255)
                                                                                                                    cTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGaa                                  <  1:1784759/150‑1 (MQ=255)
                                                                                                                     tgatgaTGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAg                                 >  1:2451902/1‑150 (MQ=255)
                                                                                                                         gatgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGa                             >  1:1278023/1‑150 (MQ=255)
                                                                                                                           tgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGAcc                           <  1:3077956/150‑1 (MQ=255)
                                                                                                                           tgTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGAcc                           <  2:2848313/150‑1 (MQ=255)
                                                                                                                            gTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCAACcg                          >  2:844941/1‑150 (MQ=255)
                                                                                                                                aCCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTa                      <  2:232593/150‑1 (MQ=255)
                                                                                                                                aCCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTa                      >  2:2749680/1‑150 (MQ=255)
                                                                                                                                 ccAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAt                     >  1:5064595/1‑150 (MQ=255)
                                                                                                                                  caacaaCATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAtt                    >  2:4647012/1‑150 (MQ=255)
                                                                                                                                  caacaaCATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAtt                    >  1:441560/1‑150 (MQ=255)
                                                                                                                                  caacaaCATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTAtt                    >  2:4480135/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)
                                                                                                                                      aacaTCCTGTCGCCAGCCAACGGTGAGCCATTCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgc                <  1:3490036/150‑1 (MQ=255)
                                                                                                                                      aacaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgc                >  1:1645645/1‑150 (MQ=255)
                                                                                                                                       acaTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCcgcg               <  2:2571696/150‑1 (MQ=255)
                                                                                                                                         aTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGcc             <  2:2484382/150‑1 (MQ=255)
                                                                                                                                          tCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCg            <  2:1856804/150‑1 (MQ=255)
                                                                                                                                          tCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCg            <  2:3775166/150‑1 (MQ=255)
                                                                                                                                           ccTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCCAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACAGCGTATTCCGCGCCgg           >  2:2354283/1‑150 (MQ=255)
                                                                                                                                             tGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCg         >  2:2373948/1‑150 (MQ=255)
                                                                                                                                                   cagccaACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAgccg   >  2:2605438/1‑150 (MQ=255)
                                                                                                                                                    agccaACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAgccgc  >  2:445539/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCAGGAAGTCGACCGCGTATTCCGCGCCGGCGAAGCCGC  >  NC_002947/543094‑543391

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: