Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,334,393 T→A 5.3% L10Q (CTG→CAG)  PP_4690 → Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,334,3930TA 5.3% 1360.4 / 44.6 415L10Q (CTG→CAG) PP_4690Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base T (181/210);  new base A (9/13);  total (190/225)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 419 total aligned reads displayed.

GCCTTGTGGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGC  >  NC_002947/5334244‑5334541
                                                                                                                                                     |                                                                                                                                                    
      tgGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGa                                                                                                                                                <  1:6758932/150‑1 (MQ=255)
        gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                              <  2:7112516/150‑1 (MQ=255)
        gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                              >  1:4773917/1‑150 (MQ=255)
        gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                              <  2:6053592/150‑1 (MQ=255)
        gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                              <  1:723923/150‑1 (MQ=255)
          gAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGc                                                                                                                                            >  2:1558295/1‑150 (MQ=255)
          gAGCGGCGGAACTGCGCACCAGAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGGTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAggga                                                                                                                                            <  2:3752842/150‑2 (MQ=255)
           aGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGcc                                                                                                                                           >  2:1443358/1‑150 (MQ=255)
           aGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGcc                                                                                                                                           >  1:2278434/1‑150 (MQ=255)
           aGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                           <  2:1208165/150‑3 (MQ=255)
              ggcggAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCaaa                                                                                                                                        <  2:4543132/150‑1 (MQ=255)
              ggcggAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCaaa                                                                                                                                        <  1:4381678/150‑1 (MQ=255)
               gcggAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAg                                                                                                                                       >  2:1841078/1‑150 (MQ=255)
                cggAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGc                                                                                                                                      >  1:6046748/1‑150 (MQ=255)
                cggAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGc                                                                                                                                      <  2:1963933/150‑1 (MQ=255)
                   aaCTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCcgc                                                                                                                                   <  2:2571342/150‑1 (MQ=255)
                     cTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCgcgc                                                                                                                                 <  1:7399970/150‑1 (MQ=255)
                     cTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCgcgc                                                                                                                                 >  1:1855217/1‑150 (MQ=255)
                     cTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTGTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCgcgc                                                                                                                                 >  2:5465626/1‑150 (MQ=255)
                      tGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGcc                                                                                                                                >  1:5747604/1‑150 (MQ=255)
                         gcACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTc                                                                                                                             <  1:2413925/150‑1 (MQ=255)
                           aCCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAg                                                                                                                           <  2:2004033/150‑1 (MQ=255)
                               tAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTa                                                                                                                       <  2:5363833/150‑1 (MQ=255)
                               tAGGCAGCCGGGCGCAGAGAATGAGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTa                                                                                                                       <  2:7134299/150‑1 (MQ=255)
                                  gCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGAc                                                                                                                    <  1:3799009/150‑1 (MQ=255)
                                   cAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACg                                                                                                                   <  2:220520/150‑1 (MQ=255)
                                       cGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCt                                                                                                               <  1:1187862/150‑1 (MQ=255)
                                          gcgcAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGa                                                                                                            >  1:1522981/1‑150 (MQ=255)
                                            gcAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAAc                                                                                                          >  2:3540874/1‑150 (MQ=255)
                                              agagAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTg                                                                                                        >  1:6903668/1‑150 (MQ=255)
                                               gagaATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGt                                                                                                       <  1:1154894/150‑1 (MQ=255)
                                                agaATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGtt                                                                                                      >  2:2137803/1‑150 (MQ=255)
                                                agaATGCGATCAATTTGATAACTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGtt                                                                                                      <  2:5979548/150‑1 (MQ=255)
                                                 gaAGGCGATCAATTTGATAATGAGTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCGTTGGACATCCCGCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTc                                                                                                     <  2:4664389/150‑1 (MQ=255)
                                                  aaTGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCg                                                                                                    <  2:909101/150‑1 (MQ=255)
                                                     gCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcg                                                                                                 <  1:4336147/150‑1 (MQ=255)
                                                      cGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgt                                                                                                <  1:5466757/150‑1 (MQ=255)
                                                        aTCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGtt                                                                                                      >  2:2364997/1‑142 (MQ=255)
                                                        aTCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgc                                                                                              <  1:4158797/150‑1 (MQ=255)
                                                        aTCAATTGGATAATTATGCGCAGTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgc                                                                                              <  2:4827948/150‑1 (MQ=255)
                                                              ttGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGcc                                                                                        >  1:6290239/1‑150 (MQ=255)
                                                                gATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCaa                                                                                      >  1:25198/1‑150 (MQ=255)
                                                                 aTAATTATTCGCATTGTAGCGTCAAGCCATCACGTTTTTCACGGGTAATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAg                                                                                     <  1:288830/150‑1 (MQ=255)
                                                                   aaTTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCa                                                                                                                            >  2:2842607/1‑109 (MQ=255)
                                                                   aaTTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGc                                                                                   >  1:6550310/1‑150 (MQ=255)
                                                                   aaTTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGc                                                                                   >  1:1279821/1‑150 (MQ=255)
                                                                     ttattCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCa                                                                                 <  2:4055703/150‑1 (MQ=255)
                                                                     ttattCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCa                                                                                 >  2:2991153/1‑150 (MQ=255)
                                                                        ttCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGgt                                                                              <  2:6379465/150‑1 (MQ=255)
                                                                         tCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGgtg                                                                             >  2:3854132/1‑150 (MQ=255)
                                                                         tCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGgtg                                                                             >  1:1000644/1‑150 (MQ=255)
                                                                            cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCGTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                          <  1:3878850/150‑1 (MQ=255)
                                                                            cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                          <  2:3965486/150‑1 (MQ=255)
                                                                            cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                          <  1:5937171/150‑1 (MQ=255)
                                                                            cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                          <  1:2577381/150‑1 (MQ=255)
                                                                             aTTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacc                                                                         >  1:23100/1‑150 (MQ=255)
                                                                             aTTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacc                                                                         >  2:651054/1‑150 (MQ=255)
                                                                              ttGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacct                                                                        >  1:5240526/1‑150 (MQ=255)
                                                                               tGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctg                                                                       >  2:7041311/1‑150 (MQ=255)
                                                                                ggAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCGACGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgt                                                                      >  2:2747348/1‑150 (MQ=255)
                                                                                ggAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgt                                                                      >  1:7143680/1‑150 (MQ=255)
                                                                                 gAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAg                                                                                >  1:6220139/1‑139 (MQ=255)
                                                                                   gCGTCAAGCCATCTCGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgtacc                                                                   <  1:775018/150‑1 (MQ=255)
                                                                                     gTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                           <  2:5634784/76‑3 (MQ=255)
                                                                                       cAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCaa                                                               <  2:5527142/150‑1 (MQ=255)
                                                                                         aGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAAcc                                                             >  2:6690501/1‑150 (MQ=255)
                                                                                            cATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCt                                                          <  1:6108841/150‑1 (MQ=255)
                                                                                              tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGc                                                        <  1:487424/150‑1 (MQ=255)
                                                                                                    gtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCAc                                                  >  2:7164063/1‑150 (MQ=255)
                                                                                                    gtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCAc                                                  >  2:5213264/1‑150 (MQ=255)
                                                                                                        cACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCggg                                              <  2:1646694/150‑1 (MQ=255)
                                                                                                         aCGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGa                                             <  2:6765029/150‑1 (MQ=255)
                                                                                                            ggTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACAt                                          >  1:5612694/1‑150 (MQ=255)
                                                                                                                aTCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCg                                      <  2:4899662/150‑1 (MQ=255)
                                                                                                                aTCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCg                                      >  2:1522080/1‑150 (MQ=255)
                                                                                                                aTCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCg                                      <  1:3248267/150‑1 (MQ=255)
                                                                                                                 tCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGc                                     >  1:5677846/1‑150 (MQ=255)
                                                                                                                   cTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTg                                   >  1:158617/1‑150 (MQ=255)
                                                                                                                         aTGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTgg                             >  1:6637181/1‑150 (MQ=255)
                                                                                                                          tGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTggg                            <  1:2448041/150‑1 (MQ=255)
                                                                                                                           gCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGc                                                                                                                          >  2:3270744/1‑55 (MQ=255)
                                                                                                                           gCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGc                           >  2:6841987/1‑150 (MQ=255)
                                                                                                                            cACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTggg                            >  2:4237216/1‑148 (MQ=255)
                                                                                                                              cTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCaga                        <  1:4559535/150‑1 (MQ=255)
                                                                                                                               ttCCTTTGTACATCCCTCGGCCTGGCCGGGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCagag                       >  1:917528/1‑150 (MQ=255)
                                                                                                                               ttCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCagag                       <  2:592800/150‑1 (MQ=255)
                                                                                                                                tCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAgg                      >  1:6375258/1‑150 (MQ=255)
                                                                                                                                  cTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAc                    >  1:3538450/1‑150 (MQ=255)
                                                                                                                                   tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAcg                   >  1:1985639/1‑150 (MQ=255)
                                                                                                                                    ttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAcgc                  <  1:3127035/150‑1 (MQ=255)
                                                                                                                                         cATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTc             >  1:2643284/1‑150 (MQ=255)
                                                                                                                                          aTCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCt            <  1:5133534/150‑1 (MQ=255)
                                                                                                                                           tCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTg           <  1:2795071/150‑1 (MQ=255)
                                                                                                                                            cccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTgg          <  2:6209193/150‑1 (MQ=255)
                                                                                                                                            cccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTgg          <  2:7406472/150‑1 (MQ=255)
                                                                                                                                             ccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGa         >  1:6437736/1‑150 (MQ=255)
                                                                                                                                             ccTCGGCCAGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGa         >  2:512502/1‑150 (MQ=255)
                                                                                                                                                  gCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACa    <  1:2151660/150‑1 (MQ=255)
                                                                                                                                                  gCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACa    <  1:5403215/150‑1 (MQ=255)
                                                                                                                                                   ccTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAg   <  1:4039409/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCTTGTGGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGATGACAGC  >  NC_002947/5334244‑5334541

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: