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,8160GT95.6% 1143.2 / 35.0 384A27A (GCG→GCT) nrfCformate‑dependent nitrite reductase, 4Fe4S subunit
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base T (191/176);  total (197/187)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 384 total aligned reads displayed.

CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGCC  >  CP009273/4279671‑4279962
                                                                                                                                                 |                                                                                                                                                  
cacaGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTaaaa                                                                                                                                                >  2:3005476/1‑150 (MQ=255)
cacaGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGaaaa                                                                                                                                                <  2:2601484/150‑1 (MQ=255)
   aGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGt                                                                                                                                             >  1:954191/1‑150 (MQ=255)
         cAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAata                                                                                                                                       <  2:469433/150‑1 (MQ=255)
            cgcACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCa                                                                                                                                    <  1:681327/150‑1 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAAt                                                                                                                                  <  1:100614/150‑1 (MQ=255)
              cACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAAt                                                                                                                                  >  2:525432/1‑150 (MQ=255)
               aCCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATg                                                                                                                                 <  1:2748997/150‑1 (MQ=255)
                  aaTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATgggg                                                                                                                              >  2:2382448/1‑150 (MQ=255)
                  aaTCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATgggg                                                                                                                              >  2:5087612/1‑150 (MQ=255)
                        aaCTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGtt                                                                                                                        <  1:2576274/150‑1 (MQ=255)
                               aCCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCAt                                                                                                                 >  1:379993/1‑150 (MQ=255)
                                cccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATCTCAATGGGGTGCGTTACGGCagg                                                                                                                >  2:831675/1‑148 (MQ=255)
                                cccGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATg                                                                                                                <  1:3977301/150‑1 (MQ=255)
                                  ctgCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGt                                                                                                              >  1:859099/3‑150 (MQ=255)
                                     cGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAGTATCAATGGGGTGCGTTACGGCATGGTGCa                                                                                                           >  2:5033733/1‑150 (MQ=255)
                                       tCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATg                                                                                                         <  1:1825925/150‑1 (MQ=255)
                                        cggTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGa                                                                                                        <  2:4385772/148‑1 (MQ=255)
                                        cTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGa                                                                                                        >  1:2221272/1‑150 (MQ=255)
                                        cTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGa                                                                                                        <  2:1812963/150‑1 (MQ=255)
                                         tGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGAc                                                                                                       >  2:4591936/1‑150 (MQ=255)
                                          gTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGAca                                                                                                      <  2:4472598/150‑2 (MQ=255)
                                          gTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACg                                                                                                      >  2:3264260/1‑150 (MQ=255)
                                           ttCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGa                                                                                                     >  2:3038276/1‑150 (MQ=255)
                                             ccGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGt                                                                                                   >  1:442675/1‑150 (MQ=255)
                                              cGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTc                                                                                                  <  2:3856729/150‑1 (MQ=255)
                                               gTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCg                                                                                                 >  2:428847/1‑150 (MQ=255)
                                                 tGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGtt                                                                                               <  1:4611172/150‑1 (MQ=255)
                                                   cTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTAt                                                                                             <  1:4676571/150‑1 (MQ=255)
                                                    ttAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATg                                                                                            >  1:1727251/1‑150 (MQ=255)
                                                       aGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCAt                                                                                         <  2:1515004/150‑1 (MQ=255)
                                                         gAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCGTGACGAGTCGTTATGCATCg                                                                                       <  1:4047290/150‑1 (MQ=255)
                                                         gAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCg                                                                                       <  1:538720/150‑1 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCt                                                                                    <  1:483897/150‑1 (MQ=255)
                                                            cagcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCt                                                                                    >  2:1986349/1‑150 (MQ=255)
                                                             agcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTg                                                                                   >  1:3874884/1‑150 (MQ=255)
                                                             agcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTg                                                                                   <  2:3371438/150‑1 (MQ=255)
                                                              gcCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGt                                                                                  <  2:4800428/150‑1 (MQ=255)
                                                               ccATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                                 <  1:1313087/150‑1 (MQ=255)
                                                               ccATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTa                                                                                 >  2:454983/1‑150 (MQ=255)
                                                                cATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTAc                                                                                <  2:398273/150‑1 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTAcc                                                                               >  1:1872262/1‑150 (MQ=255)
                                                                 aTGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTAcc                                                                               >  2:1133672/1‑150 (MQ=255)
                                                                  tGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCg                                                                              >  2:1906972/1‑150 (MQ=255)
                                                                     ccTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                           <  1:287273/150‑1 (MQ=255)
                                                                     ccTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCt                                                                           <  1:1596432/150‑1 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTg                                                                          >  2:4920906/1‑150 (MQ=255)
                                                                      cTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTg                                                                          <  1:4015352/150‑1 (MQ=255)
                                                                       tGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGt                                                                         >  1:3376844/1‑150 (MQ=255)
                                                                        ggTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                        >  1:1966676/1‑150 (MQ=255)
                                                                        ggTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTa                                                                        >  2:1677765/1‑150 (MQ=255)
                                                                         gTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTAt                                                                       >  1:4984941/1‑150 (MQ=255)
                                                                             cgtcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGAt                                                                   >  1:1447198/1‑150 (MQ=255)
                                                                              gtcgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATg                                                                  <  2:1271127/150‑1 (MQ=255)
                                                                                cgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCt                                                                >  1:2947306/1‑150 (MQ=255)
                                                                                cgCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCt                                                                >  2:406789/1‑150 (MQ=255)
                                                                                agcaaGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTAAGGCATGGTGCATGACGAGTCGTTATGCATCTGCTGTAACGCCTGTCTGGATGCt                                                                >  1:1773113/5‑150 (MQ=255)
                                                                                  ccAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCtta                                                              <  2:1428564/150‑2 (MQ=255)
                                                                                    aGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTc                                                            <  1:3614970/150‑1 (MQ=255)
                                                                                    aGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTc                                                            >  1:3211718/1‑150 (MQ=255)
                                                                                     gTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCg                                                           >  2:3797656/1‑150 (MQ=255)
                                                                                       ttCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCggg                                                         <  2:4183645/150‑1 (MQ=255)
                                                                                        tctcACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGa                                                        >  2:3606312/1‑150 (MQ=255)
                                                                                           cACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGt                                                     <  1:317396/150‑1 (MQ=255)
                                                                                            aCCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTg                                                    <  2:1176784/150‑1 (MQ=255)
                                                                                            aCCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTg                                                    <  2:37135/150‑1 (MQ=255)
                                                                                             ccGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGa                                                   <  1:2161174/150‑1 (MQ=255)
                                                                                             ccGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGa                                                   <  2:1615001/150‑1 (MQ=255)
                                                                                              cggcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGaa                                                  >  1:2298950/1‑150 (MQ=255)
                                                                                              cggcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGaa                                                  <  2:250188/150‑1 (MQ=255)
                                                                                              cggcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGaa                                                  >  1:4893172/1‑150 (MQ=255)
                                                                                              cggcgtcggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGaa                                                  <  2:4991214/150‑1 (MQ=255)
                                                                                                    cggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGt                                            >  1:4510705/1‑150 (MQ=255)
                                                                                                    cggcgtGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGt                                            >  1:76123/1‑150 (MQ=255)
                                                                                                           cTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGaa                                     >  1:1256466/1‑150 (MQ=255)
                                                                                                           cTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGaa                                     <  2:369114/150‑1 (MQ=255)
                                                                                                           cTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGaa                                     <  1:3134406/150‑1 (MQ=255)
                                                                                                            tGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAg                                    >  1:3563727/1‑150 (MQ=255)
                                                                                                            tGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAg                                    >  1:1348582/1‑150 (MQ=255)
                                                                                                              gCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGc                                  >  1:3295080/1‑150 (MQ=255)
                                                                                                                   cGTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGtctc                             <  2:2979819/150‑1 (MQ=255)
                                                                                                                    gTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGGCTGTCTGGTTGCTTGTGGGGAAGTGAAGAGGGTGCCGGTCGGCGtctgg                            >  2:2647069/1‑148 (MQ=255)
                                                                                                                    gTCAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTcg                            <  1:748569/150‑1 (MQ=255)
                                                                                                                      cAGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTcgcg                          >  1:1900523/1‑150 (MQ=255)
                                                                                                                       aGCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCgcgc                         >  1:1122931/1‑150 (MQ=255)
                                                                                                                        gCGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGcc                        <  2:636677/150‑1 (MQ=255)
                                                                                                                         cGGCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCt                       >  2:909979/1‑150 (MQ=255)
                                                                                                                           gCACCGCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGa                     >  1:4581531/1‑150 (MQ=255)
                                                                                                                                gCAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGaat                >  2:681179/1‑148 (MQ=255)
                                                                                                                                 cAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTCTGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta               <  2:2899285/150‑1 (MQ=255)
                                                                                                                                 cAGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGatta               >  1:4172522/1‑150 (MQ=255)
                                                                                                                                  aGGGCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGattat              <  1:2589920/150‑1 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTc            <  2:4494968/150‑1 (MQ=255)
                                                                                                                                    ggCGCGTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTc            <  2:320006/150‑1 (MQ=255)
                                                                                                                                      cgcgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt          >  2:2132447/1‑150 (MQ=255)
                                                                                                                                      cgcgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGt          <  2:1014776/150‑1 (MQ=255)
                                                                                                                                        cgTGGTGGCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAg        <  2:2235384/150‑1 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGcc  <  1:5087612/150‑1 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGcc  <  2:586041/150‑1 (MQ=255)
                                                                                                                                              ggCTAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGcc  <  2:785529/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACAGCGATCAGCGCACCAATCCGAACTTTAACCCGGCGTCTGTTCCGTTGCTTAAGGAGCAGCCATGACCTGGTCTCGTCGCCAGTTTCTCACCGGCGTCGGCGTGCTGGCAGCCGTCAGCGGCACCGCAGGGCGCGTGGTGGCGAAAACGTTGAATATCAATGGGGTGCGTTACGGCATGGTGCATGACGAGTCGTTATGCATCGGCTGTACCGCCTGTATGGATGCTTGTCGGGAAGTGAACAAGGTGCCGGAAGGCGTCTCGCGCCTGACGATTATTCGTAGCGAGCC  >  CP009273/4279671‑4279962

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: ‑