Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,153,033 C→A 9.4% L491L (CTG→CTT ygiQ ← Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,153,0330CA 9.4% 758.9 / 32.7 216L491L (CTG→CTTygiQRadical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (93/102);  new base A (9/12);  total (102/114)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

CATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGC  >  CP009273/3152886‑3153180
                                                                                                                                                   |                                                                                                                                                   
cATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATc                                                                                                                                                     <  2:2414757/148‑1 (MQ=255)
cATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATc                                                                                                                                                     >  1:1683668/1‑148 (MQ=255)
  tGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAg                                                                                                                                                   >  1:2409873/1‑148 (MQ=255)
     ccGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGAt                                                                                                                                                <  2:1317740/148‑1 (MQ=255)
     ccGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGAt                                                                                                                                                <  1:2657241/148‑1 (MQ=255)
      cGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATc                                                                                                                                               <  2:2667007/148‑1 (MQ=255)
       ggCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCt                                                                                                                                              >  2:2910306/1‑148 (MQ=255)
       ggCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCt                                                                                                                                              <  2:1829224/148‑1 (MQ=255)
        gCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCtt                                                                                                                                             >  2:2262858/1‑148 (MQ=255)
         cTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCttt                                                                                                                                            <  1:2257662/148‑1 (MQ=255)
           tcatcaTCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCttttt                                                                                                                                          <  2:684768/148‑1 (MQ=255)
           tcatcaTCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCttttt                                                                                                                                          >  2:3157111/1‑148 (MQ=255)
            catcatCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTa                                                                                                                                         >  2:1190537/1‑148 (MQ=255)
             atcatcTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTaa                                                                                                                                        <  2:1802492/148‑1 (MQ=255)
              tcatcTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAAt                                                                                                                                       >  1:1445617/1‑148 (MQ=255)
                    tCGATACCGGCCCTTCTTCGGTCTGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAAGGCCttt                                                                                                                                 <  2:2907030/148‑1 (MQ=255)
                    tCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCttt                                                                                                                                 <  2:1343311/148‑1 (MQ=255)
                     cGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTc                                                                                                                                <  1:2342608/148‑1 (MQ=255)
                     cGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTc                                                                                                                                <  2:2967153/148‑1 (MQ=255)
                       aTAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTc                                                                                                                                <  2:2227399/146‑1 (MQ=255)
                       aTAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAg                                                                                                                              <  2:2849206/148‑1 (MQ=255)
                                 ttcttcGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCa                                                                                                                    >  1:1576054/1‑148 (MQ=255)
                                 ttcttcGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCa                                                                                                                    <  1:1182384/148‑1 (MQ=255)
                                    ttcGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACgg                                                                                                                 <  1:3138406/148‑1 (MQ=255)
                                       ggTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGa                                                                                                              <  2:499411/148‑1 (MQ=255)
                                          aTGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATag                                                                                                           <  1:219289/148‑1 (MQ=255)
                                          aTGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATag                                                                                                           <  1:379872/148‑1 (MQ=255)
                                            gTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATagag                                                                                                         <  1:1360097/148‑1 (MQ=255)
                                            gTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATagag                                                                                                         <  1:444999/148‑1 (MQ=255)
                                             ttCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTCATGCCTTTCAGATCACGCGCACGGCGATAGAgg                                                                                                        <  1:3157111/148‑1 (MQ=255)
                                              tCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGt                                                                                                       >  1:102016/1‑148 (MQ=255)
                                                  gggCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGAt                                                                                                   >  1:1837995/1‑148 (MQ=255)
                                                   ggCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATc                                                                                                  >  2:2725612/1‑148 (MQ=255)
                                                   ggCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATc                                                                                                  >  1:458717/1‑148 (MQ=255)
                                                   ggCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATc                                                                                                  >  1:72779/1‑148 (MQ=255)
                                                    gCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCg                                                                                                 <  1:2109326/148‑1 (MQ=255)
                                                     cAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGt                                                                                                >  1:2200810/1‑148 (MQ=255)
                                                      aaTCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGtt                                                                                               <  2:2695126/148‑1 (MQ=255)
                                                        tCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtgg                                                                                             >  1:903115/1‑148 (MQ=255)
                                                         cTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggt                                                                                            >  2:3258473/1‑148 (MQ=255)
                                                         cTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggt                                                                                            <  2:2552432/148‑1 (MQ=255)
                                                          ttCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGTATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTAAGATCACGCGCACGGCGAGAGAGGTTGATCGTtggtt                                                                                           >  2:3109596/1‑148 (MQ=255)
                                                          ttCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggtt                                                                                           <  2:877006/148‑1 (MQ=255)
                                                            cAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcg                                                                                         <  1:1916831/148‑1 (MQ=255)
                                                              gATAACCGCCGACGTGATGGGTCGCCAGTTATTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtg                                                                                       <  1:2299255/148‑1 (MQ=255)
                                                               aTAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtgg                                                                                      >  1:2674056/1‑148 (MQ=255)
                                                                tAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggt                                                                                     <  1:831232/148‑1 (MQ=255)
                                                                tAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggt                                                                                     <  2:52833/148‑1 (MQ=255)
                                                                 aaCCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggtt                                                                                    <  2:2514331/148‑1 (MQ=255)
                                                                  aCCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggttc                                                                                   >  1:1507385/1‑148 (MQ=255)
                                                                  aCCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggttc                                                                                   >  1:321911/1‑148 (MQ=255)
                                                                  aCCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTtggttcgtggttc                                                                                   <  2:343623/148‑1 (MQ=255)
                                                                        gACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTcc                                                                             <  2:1644361/148‑1 (MQ=255)
                                                                          cGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCAt                                                                           >  1:1409058/1‑148 (MQ=255)
                                                                           gTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCAtg                                                                          <  2:691453/148‑1 (MQ=255)
                                                                             gATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCAtgtg                                                                        >  2:2459603/1‑148 (MQ=255)
                                                                              aTGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGc                                                                       >  1:2788221/1‑148 (MQ=255)
                                                                               tGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCg                                                                      >  2:2881138/1‑148 (MQ=255)
                                                                                 ggTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGa                                                                    >  1:1192313/1‑148 (MQ=255)
                                                                                      ccAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAat                                                               <  1:2977131/148‑1 (MQ=255)
                                                                                      ccAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAat                                                               <  2:3221943/148‑1 (MQ=255)
                                                                                       cAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAata                                                              >  2:1427150/1‑148 (MQ=255)
                                                                                        aGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAatat                                                             >  2:2593892/1‑148 (MQ=255)
                                                                                         gTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATc                                                            >  1:112668/1‑148 (MQ=255)
                                                                                         gTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATc                                                            <  1:1000122/148‑1 (MQ=255)
                                                                                              ttGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGAt                                                       <  2:173852/148‑1 (MQ=255)
                                                                                                gATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATaa                                                     <  1:2060933/148‑1 (MQ=255)
                                                                                                gATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATaa                                                     <  2:262402/148‑1 (MQ=255)
                                                                                                gATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATaa                                                     <  2:3165282/148‑1 (MQ=255)
                                                                                                  tataGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAAc                                                   >  1:2875582/1‑148 (MQ=255)
                                                                                                        gcgGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGa                                             <  1:1049836/148‑1 (MQ=255)
                                                                                                         cgGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGAc                                            >  1:791212/1‑148 (MQ=255)
                                                                                                           gATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACaa                                          <  1:3069659/148‑1 (MQ=255)
                                                                                                           gATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGACAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACaa                                          <  2:401291/148‑1 (MQ=255)
                                                                                                            aTCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACaaa                                         <  2:1556477/148‑1 (MQ=255)
                                                                                                             tcttctACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAAc                                        >  2:559235/1‑148 (MQ=255)
                                                                                                               ttctACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACgg                                      <  1:976315/148‑1 (MQ=255)
                                                                                                                      gCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACaaa                               <  2:1449925/148‑2 (MQ=255)
                                                                                                                      gCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAg                               <  1:2606356/148‑1 (MQ=255)
                                                                                                                      gCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAg                               >  1:2680949/1‑148 (MQ=255)
                                                                                                                      gCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAg                               <  2:1683668/148‑1 (MQ=255)
                                                                                                                       cTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAACCGCACGACAAACGGCGACAAGt                              <  1:872799/148‑1 (MQ=255)
                                                                                                                         aTGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATAGCCGGATAAACGCACGACAAACGGCGACAAGttt                            >  1:1960226/1‑148 (MQ=255)
                                                                                                                         aTGTCATAACGCACACCAGAGGCAATAAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGttt                            >  2:1040820/1‑148 (MQ=255)
                                                                                                                          tGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTg                           <  2:1021865/148‑1 (MQ=255)
                                                                                                                          tGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTg                           <  2:3022205/148‑1 (MQ=255)
                                                                                                                               tAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAg                      >  2:213746/1‑148 (MQ=255)
                                                                                                                               tAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAg                      >  1:2720707/1‑148 (MQ=255)
                                                                                                                                   gCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGcgcg                  <  2:2066376/148‑1 (MQ=255)
                                                                                                                                      cacCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTgg               >  1:645801/1‑148 (MQ=255)
                                                                                                                                        ccAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCg             <  2:1165414/148‑1 (MQ=255)
                                                                                                                                         cAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGa            <  1:1308617/148‑1 (MQ=255)
                                                                                                                                          agagGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGAt           <  1:1064906/148‑1 (MQ=255)
                                                                                                                                          agagGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGAt           <  2:2276727/148‑1 (MQ=255)
                                                                                                                                          agagGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGAt           >  1:605676/1‑148 (MQ=255)
                                                                                                                                          agagGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGAt           <  2:2277479/148‑1 (MQ=255)
                                                                                                                                          agagGCAATCAGGATCTTTTTAAAGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGAt           >  1:1672374/1‑148 (MQ=255)
                                                                                                                                             ggCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTg        >  2:2074194/1‑148 (MQ=255)
                                                                                                                                             ggCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTg        >  2:2283124/1‑148 (MQ=255)
                                                                                                                                              gCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGc       >  2:2027691/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATAGCGCGGATCTTCTACGGCTATGTCATAACGCACACCAGAGGCAATCAGGATCTTTTTAATGCCTTTCAGATCACGCGCACGGCGATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGC  >  CP009273/3152886‑3153180

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: