Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,279,816 G→T 33.2% A27A (GCG→GCT nrfC → formate‑dependent nitrite reductase, 4Fe4S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,279,8160GT33.2% 374.8 / 171.3 217A27A (GCG→GCTnrfCformate‑dependent nitrite reductase, 4Fe4S subunit
Reads supporting (aligned to +/- strand):  ref base G (73/71);  new base T (28/44);  total (101/116)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGC  >  CP009273/4279671‑4279961
                                                                                                                                                 |                                                                                                                                                 
cacaGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGaa                                                                                                                                                 <  1:1995638/148‑1 (MQ=255)
cacaGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGaa                                                                                                                                                 >  1:2259515/1‑148 (MQ=255)
  caGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGaaaa                                                                                                                                               <  1:3088907/148‑1 (MQ=255)
   aGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGTGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAAc                                                                                                                                              >  1:2040374/1‑148 (MQ=255)
     cgctCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                            <  2:1979531/145‑1 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGt                                                                                                                                            <  1:3309133/148‑1 (MQ=255)
       aTCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTg                                                                                                                                          <  1:56347/148‑1 (MQ=255)
       aTCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTg                                                                                                                                          >  2:220157/1‑148 (MQ=255)
        tCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGa                                                                                                                                         <  1:1521553/148‑1 (MQ=255)
         cAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGaa                                                                                                                                        <  1:2039611/148‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAat                                                                                                                                       <  2:1174788/148‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAat                                                                                                                                       <  2:1668813/148‑1 (MQ=255)
           gcgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAata                                                                                                                                      >  1:537362/1‑148 (MQ=255)
            cgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAatat                                                                                                                                     >  2:751113/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCa                                                                                                                                   <  2:86021/148‑1 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCa                                                                                                                                   <  2:1718680/148‑1 (MQ=255)
               aCCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCaa                                                                                                                                  <  1:2670485/148‑1 (MQ=255)
                ccAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAAt                                                                                                                                 <  2:501410/148‑1 (MQ=255)
                 cAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAAAGTTGAATATCAATg                                                                                                                                >  1:616506/1‑148 (MQ=255)
                 cAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATg                                                                                                                                <  2:2112069/148‑1 (MQ=255)
                 cAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCCTCAGCGGCACCGCAGGGCGCCAGGAGGCGAAATCGTTGAATATCaaag                                                                                                                                >  1:2470468/1‑146 (MQ=255)
                  aaTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATgg                                                                                                                               <  1:1426123/148‑1 (MQ=255)
                   aTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATggg                                                                                                                              <  2:2162175/148‑1 (MQ=255)
                   aTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATggg                                                                                                                              <  2:1413133/148‑1 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATgggg                                                                                                                             >  2:2436558/1‑148 (MQ=255)
                        aaCTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCg                                                                                                                         <  1:397717/148‑1 (MQ=255)
                           tttAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTa                                                                                                                      >  2:2379390/1‑148 (MQ=255)
                              aaCCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACgg                                                                                                                   >  1:1935746/1‑148 (MQ=255)
                               aCCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGc                                                                                                                  >  2:2108197/1‑148 (MQ=255)
                                cccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCa                                                                                                                 >  1:363116/1‑148 (MQ=255)
                                 ccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCAt                                                                                                                >  1:2916551/1‑148 (MQ=255)
                                  cGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATg                                                                                                               >  1:373678/1‑148 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                              >  2:1847865/1‑148 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                              <  2:1233989/148‑1 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                              >  1:2989046/1‑148 (MQ=255)
                                    gCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGt                                                                                                             >  2:195090/1‑148 (MQ=255)
                                       tCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCa                                                                                                          <  1:983847/148‑1 (MQ=255)
                                          gTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGa                                                                                                       <  1:2927878/148‑1 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                     >  2:1285873/1‑148 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                     >  2:44140/1‑148 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                     >  2:2999741/1‑148 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                  <  1:748767/148‑1 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                  <  1:751113/148‑1 (MQ=255)
                                                  gCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGt                                                                                               <  1:3178074/148‑1 (MQ=255)
                                                   cTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGtt                                                                                              >  2:1772037/1‑148 (MQ=255)
                                                       aGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGc                                                                                          >  2:1072784/1‑148 (MQ=255)
                                                           gcagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCg                                                                                      >  2:3468860/1‑148 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCgg                                                                                     >  2:721648/1‑148 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCgg                                                                                     >  2:1504400/1‑148 (MQ=255)
                                                               ccATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTg                                                                                  <  1:1897146/148‑1 (MQ=255)
                                                                     ccTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGc                                                                            <  2:2824673/148‑1 (MQ=255)
                                                                     ccTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGc                                                                            <  2:2528460/148‑1 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGcc                                                                           >  2:2527500/1‑148 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGcc                                                                           >  2:3523539/1‑148 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGcc                                                                           >  2:3261217/1‑148 (MQ=255)
                                                                       tGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                          >  2:1163508/1‑148 (MQ=255)
                                                                       tGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                          >  1:1625821/1‑148 (MQ=255)
                                                                        ggTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGTCTAAAACGTTGAATATCAATGGGGTGCGTAACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTg                                                                         >  2:849639/1‑148 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGt                                                                        <  1:2911960/148‑1 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGt                                                                        >  2:363294/1‑148 (MQ=255)
                                                                          tctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                       >  2:671618/1‑148 (MQ=255)
                                                                          tctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                       <  2:710239/148‑1 (MQ=255)
                                                                           ctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTAt                                                                      >  1:2550259/1‑148 (MQ=255)
                                                                             cgtcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATgg                                                                    >  1:2709044/1‑148 (MQ=255)
                                                                                cgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATg                                                                 >  2:1391701/1‑148 (MQ=255)
                                                                                    aGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCGGTATGGATGCTTg                                                             >  2:1226324/1‑148 (MQ=255)
                                                                                     gTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGt                                                            >  1:2314814/1‑148 (MQ=255)
                                                                                     gTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGt                                                            <  2:625827/148‑1 (MQ=255)
                                                                                           cACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGaa                                                      <  1:1096410/148‑1 (MQ=255)
                                                                                            aCCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAg                                                     <  1:2770838/148‑1 (MQ=255)
                                                                                                     ggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAAAGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAgg                                            <  1:130731/148‑1 (MQ=255)
                                                                                                     ggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAgg                                            <  1:3044047/148‑1 (MQ=255)
                                                                                                       cgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTg                                          <  2:12184/148‑1 (MQ=255)
                                                                                                          gCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCg                                       >  1:1781709/1‑148 (MQ=255)
                                                                                                           cTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCgg                                      <  2:168664/148‑1 (MQ=255)
                                                                                                             ggCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGTCGGaa                                    <  2:2727930/148‑1 (MQ=255)
                                                                                                                aGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGc                                 <  2:2693817/148‑1 (MQ=255)
                                                                                                                aGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGc                                 <  2:2693963/148‑1 (MQ=255)
                                                                                                                 gCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCg                                >  2:1066528/1‑148 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAACACGTTGACTATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                              <  2:3542928/148‑1 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                              <  1:3469094/148‑1 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                              <  1:3468536/148‑1 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                              <  2:750815/148‑1 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                              <  1:2655395/148‑1 (MQ=255)
                                                                                                                    gTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtct                             <  1:3002179/148‑1 (MQ=255)
                                                                                                                         cGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCgcgc                        >  2:3023909/1‑148 (MQ=255)
                                                                                                                          ggCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGcc                       <  2:2326358/148‑1 (MQ=255)
                                                                                                                          ggCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGcc                       <  1:3553024/148‑1 (MQ=255)
                                                                                                                            cACCGCAGGGCGCGTGGTGGCTCAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTAGTTATGCATCGGCTGTACCGCCTGTATGGATGATTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTg                     >  2:2170774/1‑148 (MQ=255)
                                                                                                                               cGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACg                  >  1:1912513/1‑148 (MQ=255)
                                                                                                                                 cAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGat                >  2:750152/1‑148 (MQ=255)
                                                                                                                                   gggCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAAAAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta              <  1:2379390/148‑1 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat             >  2:1298136/1‑148 (MQ=255)
                                                                                                                                     gcgcGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattatt            <  2:3157394/148‑1 (MQ=255)
                                                                                                                                        cgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt         <  1:841977/148‑1 (MQ=255)
                                                                                                                                        cgTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt         <  2:764096/148‑1 (MQ=255)
                                                                                                                                          tggtggCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAg       >  1:2805803/1‑148 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAg   <  2:2700090/148‑1 (MQ=255)
                                                                                                                                               gCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGc  >  2:37360/1‑148 (MQ=255)
                                                                                                                                               gCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGc  >  2:74777/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGC  >  CP009273/4279671‑4279961

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: