Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,279,816 G→T 73.8% 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,8160GT73.8% 270.4 / 182.1 214A27A (GCG→GCTnrfCformate‑dependent nitrite reductase, 4Fe4S subunit
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base T (64/94);  total (89/125)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

Only 100 of 214 total aligned reads displayed.

CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAG  >  CP009273/4279671‑4279960
                                                                                                                                                 |                                                                                                                                                
cacaGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTaa                                                                                                                                                <  1:1285435/148‑1 (MQ=255)
 acaGCGATCAGCGCACCACTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGaaa                                                                                                                                               <  2:1780558/148‑1 (MQ=255)
  caGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTaaaa                                                                                                                                              >  1:1021952/1‑148 (MQ=255)
    gCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACg                                                                                                                                            <  2:2326910/148‑1 (MQ=255)
    gCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACg                                                                                                                                            <  1:1390417/148‑1 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                           <  2:830054/148‑1 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                           >  1:1140795/1‑148 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                           <  2:3430425/148‑1 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                           >  1:2089232/1‑148 (MQ=255)
       aTCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTg                                                                                                                                         <  1:1498727/148‑1 (MQ=255)
        tCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGa                                                                                                                                        >  1:873217/1‑148 (MQ=255)
        tCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGa                                                                                                                                        >  2:763138/1‑148 (MQ=255)
        tCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGa                                                                                                                                        <  2:2073735/148‑1 (MQ=255)
         cAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGaa                                                                                                                                       <  2:2980045/148‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAat                                                                                                                                      <  1:1210111/148‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAat                                                                                                                                      <  2:2654154/148‑1 (MQ=255)
           gcgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAata                                                                                                                                     >  1:464781/1‑148 (MQ=255)
            cgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGGGGTGGCTAAAACGTTGAatat                                                                                                                                    >  1:1660826/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                  >  1:226177/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                  >  1:914953/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                  >  2:1741735/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                  >  2:2128819/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCa                                                                                                                                  <  2:3277183/148‑1 (MQ=255)
                ccAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAAt                                                                                                                                >  1:174741/1‑148 (MQ=255)
                 cAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATg                                                                                                                               <  2:2881104/148‑1 (MQ=255)
                   aTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATggg                                                                                                                             >  1:1233346/1‑148 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATgggg                                                                                                                            >  1:575352/1‑148 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATgggg                                                                                                                            <  2:2198672/148‑1 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCAGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCTTCAGCGGCACCTCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATgggg                                                                                                                            <  2:436955/148‑1 (MQ=255)
                      cGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTg                                                                                                                          <  1:1874426/148‑1 (MQ=255)
                       gAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGc                                                                                                                         <  2:3324262/148‑1 (MQ=255)
                       gAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGc                                                                                                                         <  2:3163057/148‑1 (MQ=255)
                         aCTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGt                                                                                                                       <  2:114763/148‑1 (MQ=255)
                          cTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGtt                                                                                                                      >  2:1520239/1‑148 (MQ=255)
                          cTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGtt                                                                                                                      >  2:2896412/1‑148 (MQ=255)
                           tttAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTa                                                                                                                     <  2:3202895/148‑1 (MQ=255)
                                cccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCa                                                                                                                <  1:2007947/148‑1 (MQ=255)
                                 ccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCCTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCAt                                                                                                               <  2:1902813/148‑1 (MQ=255)
                                  cGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATg                                                                                                              >  2:2580922/1‑148 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                             <  1:2534832/148‑1 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                             <  1:2533580/148‑1 (MQ=255)
                                    gCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGt                                                                                                            >  2:756385/1‑148 (MQ=255)
                                     cGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTg                                                                                                           <  1:763138/148‑1 (MQ=255)
                                     cGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTg                                                                                                           <  1:869218/148‑1 (MQ=255)
                                     cGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTg                                                                                                           <  1:1455797/148‑1 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                    <  1:2349247/148‑1 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                    >  1:809565/1‑148 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                    <  2:524689/148‑1 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGAGTTACGGCATGGTGCATGACg                                                                                                    >  1:425602/1‑148 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                 <  2:2625199/148‑1 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                 <  1:1490739/148‑1 (MQ=255)
                                                ttGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTc                                                                                                <  1:1074123/148‑1 (MQ=255)
                                                ttGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTc                                                                                                <  1:2983237/148‑1 (MQ=255)
                                                    ttAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTa                                                                                            <  2:3341038/148‑1 (MQ=255)
                                                         gAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCAt                                                                                       >  2:1996717/1‑148 (MQ=255)
                                                          agcagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATc                                                                                      >  2:1875024/1‑148 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCgg                                                                                    >  1:2790617/1‑148 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCgg                                                                                    <  1:971932/148‑1 (MQ=255)
                                                               ccATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTg                                                                                 <  2:1282819/148‑1 (MQ=255)
                                                                cATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGt                                                                                <  1:226541/148‑1 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGTCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                               >  1:3126959/1‑148 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                               <  2:174741/148‑1 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTCCCGCCTGt                                                                       >  2:2484882/1‑148 (MQ=255)
                                                                          tctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                      <  1:3312622/148‑1 (MQ=255)
                                                                          tctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                      <  1:3493164/148‑1 (MQ=255)
                                                                              gtcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGa                                                                  >  2:3033127/1‑148 (MQ=255)
                                                                              gtcgCCAGGTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGa                                                                  <  1:389535/148‑1 (MQ=255)
                                                                               tcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGAt                                                                 <  2:226177/148‑1 (MQ=255)
                                                                                   cAGTATCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCAAGACGAGTCGTTATGCATATGCAGAACCGCCTGTATGGAAGCtt                                                             >  2:1314295/1‑148 (MQ=255)
                                                                                     gTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGt                                                           <  1:1539828/148‑1 (MQ=255)
                                                                                      tttCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTc                                                          >  1:2873491/1‑148 (MQ=255)
                                                                                      tttCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTc                                                          >  2:3081629/1‑148 (MQ=255)
                                                                                         ctcACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCggg                                                       <  2:702420/148‑1 (MQ=255)
                                                                                          tcACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCGTGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGa                                                      <  1:2860224/148‑1 (MQ=255)
                                                                                           cACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGaa                                                     <  1:924964/148‑1 (MQ=255)
                                                                                            aCCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAg                                                    <  1:2008003/148‑1 (MQ=255)
                                                                                             ccGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGt                                                   >  2:1796318/1‑148 (MQ=255)
                                                                                                gcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGaa                                                <  1:3437373/148‑1 (MQ=255)
                                                                                                  gtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGAATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACa                                              <  1:605737/148‑1 (MQ=255)
                                                                                                      gcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGt                                          >  2:601486/1‑148 (MQ=255)
                                                                                                          gCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCg                                      >  1:3202453/1‑148 (MQ=255)
                                                                                                              gCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAg                                  >  2:1399370/1‑148 (MQ=255)
                                                                                                               cAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAgg                                 <  2:2402277/148‑1 (MQ=255)
                                                                                                                 gCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCg                               <  2:1959857/148‑1 (MQ=255)
                                                                                                                  ccGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGt                              <  1:2128819/148‑1 (MQ=255)
                                                                                                                  ccGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGt                              >  2:1276178/1‑148 (MQ=255)
                                                                                                                      cAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTcg                          >  1:2312366/1‑148 (MQ=255)
                                                                                                                           gCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCt                     <  2:1728693/148‑1 (MQ=255)
                                                                                                                            cACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTg                    >  1:110940/1‑148 (MQ=255)
                                                                                                                             aCCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGa                   <  1:2026876/148‑1 (MQ=255)
                                                                                                                                gCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGa                <  1:2729918/148‑1 (MQ=255)
                                                                                                                                gCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGa                <  2:1698140/148‑1 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat            <  2:873217/148‑1 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat            >  1:3140950/1‑148 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat            <  1:553026/148‑1 (MQ=255)
                                                                                                                                         gtggtgGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTa       <  1:3274827/148‑1 (MQ=255)
                                                                                                                                          tggtggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAg      <  1:981623/148‑1 (MQ=255)
                                                                                                                                           ggtggCGAAAACGTTGAATATCAATGGGGTGCGTTACGACATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGc     >  2:3284000/1‑148 (MQ=255)
                                                                                                                                           ggtggCGAAAACGTTGAATATCAATGGGGTGCGTTACGACATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCGGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGc     >  2:1626579/1‑148 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAg  <  2:731145/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAG  >  CP009273/4279671‑4279960

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: