Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,279,816 G→T 100% 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,8160GT99.5% 644.7 / ‑4.3 200A27A (GCG→GCTnrfCformate‑dependent nitrite reductase, 4Fe4S subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (86/113);  total (86/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAG  >  CP009273/4279671‑4279960
                                                                                                                                                 |                                                                                                                                                
  caGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTaaaa                                                                                                                                              >  2:1373221/1‑148 (MQ=255)
     cGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                           <  2:2378622/148‑1 (MQ=255)
      gATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGtt                                                                                                                                          <  2:528237/148‑1 (MQ=255)
       aTCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTg                                                                                                                                         <  1:951518/148‑1 (MQ=255)
        tCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGa                                                                                                                                        <  2:1346348/148‑1 (MQ=255)
         cAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGaa                                                                                                                                       >  1:3339233/1‑148 (MQ=255)
         atgcgcACGAATGCGAACTTTAACTCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCAGCGGCGTCGGCGTGGTGGCTGCCGTCAGCGGCTCCGCAGGGCGCGTGGTGGCTAAATGGTTGaa                                                                                                                                       <  2:2021803/146‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAat                                                                                                                                      <  2:1008294/148‑1 (MQ=255)
          aGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAat                                                                                                                                      <  1:2302021/148‑1 (MQ=255)
            cgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGaggtg                                                                                                                                    >  2:3207673/1‑144 (MQ=255)
            cgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAatat                                                                                                                                    >  2:3252935/1‑148 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                  <  1:3081110/148‑1 (MQ=255)
               aCCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCaa                                                                                                                                 <  1:794737/148‑1 (MQ=255)
               aCCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCaa                                                                                                                                 <  1:1358609/148‑1 (MQ=255)
                 cAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATg                                                                                                                               >  2:925547/1‑148 (MQ=255)
                  aaTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATgg                                                                                                                              <  2:449010/148‑1 (MQ=255)
                   aTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCANCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATggg                                                                                                                             <  1:1783921/148‑1 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTCAAACGTTGAATATCAATgggg                                                                                                                            >  2:2997802/1‑148 (MQ=255)
                    tCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATgggg                                                                                                                            >  2:3158605/1‑148 (MQ=255)
                         aCTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGt                                                                                                                       >  1:1214764/1‑148 (MQ=255)
                         aCTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGt                                                                                                                       >  1:1215104/1‑148 (MQ=255)
                           tttAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTa                                                                                                                     >  2:3343255/1‑148 (MQ=255)
                                cccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCa                                                                                                                <  2:1506448/148‑1 (MQ=255)
                                  cGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATg                                                                                                              <  1:1615627/148‑1 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                             <  2:1909255/148‑1 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                             <  2:964357/148‑1 (MQ=255)
                                   ggCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATgg                                                                                                             <  2:2842609/148‑1 (MQ=255)
                                     cGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTg                                                                                                           <  1:2111715/148‑1 (MQ=255)
                                       tCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAACACGTTGAAGATCAATGGGGTGCGTTACGGCATGGTGCa                                                                                                         <  2:2297956/148‑1 (MQ=255)
                                       tCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCa                                                                                                         <  2:2298033/148‑1 (MQ=255)
                                         tGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATg                                                                                                       >  1:267998/1‑148 (MQ=255)
                                         tGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATg                                                                                                       >  2:1905840/1‑148 (MQ=255)
                                          gTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGa                                                                                                      <  2:1275897/148‑1 (MQ=255)
                                           ttCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGAc                                                                                                     >  1:1526991/1‑148 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                    >  2:1075832/1‑148 (MQ=255)
                                            tCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                    >  2:1075670/1‑148 (MQ=255)
                                             ccGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGa                                                                                                   <  1:2527846/148‑1 (MQ=255)
                                             ccGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGa                                                                                                   <  2:2261845/148‑1 (MQ=255)
                                             ccGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGa                                                                                                   <  2:1612706/148‑1 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                 >  1:762320/1‑148 (MQ=255)
                                                   cTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGtt                                                                                             <  2:584961/148‑1 (MQ=255)
                                                   cTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGtt                                                                                             >  1:1632861/1‑148 (MQ=255)
                                                    ttAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGTTGCGTTACGGCATGGTGCATGACGAGTCGTTa                                                                                            <  2:3236659/148‑1 (MQ=255)
                                                    ttAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTa                                                                                            <  1:249397/148‑1 (MQ=255)
                                                     tAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTAt                                                                                           >  2:2945337/1‑148 (MQ=255)
                                                       aGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGc                                                                                         >  1:546871/1‑148 (MQ=255)
                                                        ggAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCa                                                                                        >  1:1160140/1‑148 (MQ=255)
                                                        ggAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCa                                                                                        >  2:3137924/1‑148 (MQ=255)
                                                        ggAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCa                                                                                        >  1:3489436/1‑148 (MQ=255)
                                                         gAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCAt                                                                                       >  1:284342/1‑148 (MQ=255)
                                                         gAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCAt                                                                                       <  2:3303886/148‑1 (MQ=255)
                                                          agcagcCATGACCTGGTCTCGTCGCCAGTTTCTCACAGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATc                                                                                      <  1:952346/148‑1 (MQ=255)
                                                          agcagcCATGACCTGGTCTCGTCGCCAGTTTCTCACAGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATc                                                                                      <  1:952387/148‑1 (MQ=255)
                                                               ccATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTg                                                                                 <  2:2236681/148‑1 (MQ=255)
                                                                cATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGt                                                                                <  1:2410095/148‑1 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTTGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                               <  1:1480996/148‑1 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                               <  1:1929344/148‑1 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                               <  1:2530456/148‑1 (MQ=255)
                                                                     ccTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGc                                                                           <  1:157001/148‑1 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGcc                                                                          <  1:2570525/148‑1 (MQ=255)
                                                                       tGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                         >  1:2799817/1‑148 (MQ=255)
                                                                       tGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                         >  1:1048609/1‑148 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGt                                                                       <  2:790177/148‑1 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGt                                                                       <  1:454289/148‑1 (MQ=255)
                                                                           ctcGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTAt                                                                     <  1:2028845/148‑1 (MQ=255)
                                                                             cgtcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTGGCATATGAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATgg                                                                   <  2:578731/148‑1 (MQ=255)
                                                                                cgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATg                                                                <  2:342963/148‑1 (MQ=255)
                                                                                 gCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGc                                                               <  1:2402583/148‑1 (MQ=255)
                                                                                    aGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTg                                                            <  1:3058640/148‑1 (MQ=255)
                                                                                     gTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGt                                                           <  2:3419122/148‑1 (MQ=255)
                                                                                           cACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGaa                                                     >  2:1737092/1‑148 (MQ=255)
                                                                                           cACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGaa                                                     <  2:3273233/148‑1 (MQ=255)
                                                                                              cggcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTg                                                  >  1:387992/1‑148 (MQ=255)
                                                                                                 cgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAAc                                               >  2:2586930/1‑148 (MQ=255)
                                                                                                   tcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACaa                                             <  2:1048609/148‑1 (MQ=255)
                                                                                                   tcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACaa                                             <  2:1523352/148‑1 (MQ=255)
                                                                                                    cggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGACGTGAACAAg                                            >  2:1396708/1‑148 (MQ=255)
                                                                                                      gcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGt                                          <  2:632672/148‑1 (MQ=255)
                                                                                                            tGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGa                                    <  2:2522587/148‑1 (MQ=255)
                                                                                                             ggCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGaa                                   >  2:3244857/1‑148 (MQ=255)
                                                                                                                aGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGc                                >  2:2984704/1‑148 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtc                             <  2:1586249/148‑1 (MQ=255)
                                                                                                                         cGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCgcgc                       <  2:2905052/148‑1 (MQ=255)
                                                                                                                         cGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCgcgc                       >  1:3062640/1‑148 (MQ=255)
                                                                                                                         cGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCAGGAAGGCGTCTCgcgc                       >  2:220473/1‑148 (MQ=255)
                                                                                                                            cACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTg                    <  1:925547/148‑1 (MQ=255)
                                                                                                                              ccGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGAc                  >  1:963822/1‑148 (MQ=255)
                                                                                                                                 cAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGat               >  1:1615200/1‑148 (MQ=255)
                                                                                                                                   gggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta             <  2:353525/148‑1 (MQ=255)
                                                                                                                                   gggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta             <  2:354035/148‑1 (MQ=255)
                                                                                                                                   gggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta             >  2:2862849/1‑148 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat            <  2:2360658/148‑1 (MQ=255)
                                                                                                                                     gcgcGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattatt           <  2:1208831/148‑1 (MQ=255)
                                                                                                                                        cgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt        >  1:3368046/1‑148 (MQ=255)
                                                                                                                                        cgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt        >  2:2740766/1‑148 (MQ=255)
                                                                                                                                        cgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt        <  1:2914648/148‑1 (MQ=255)
                                                                                                                                            gtggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCg    <  1:1342911/148‑1 (MQ=255)
                                                                                                                                            gtggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCg    >  1:3467182/1‑148 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAg  >  2:478725/1‑148 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAg  <  1:3241955/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: