Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,412,205 C→A 100% D151Y (GAT→TAT)  yfcD ← putative NUDIX hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,412,2050CA97.0% 718.5 / 11.5 236D151Y (GAT→TAT) yfcDputative NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (142/86);  total (145/91)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 237 total aligned reads displayed.

AGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTG  >  CP009273/2412061‑2412349
                                                                                                                                                |                                                                                                                                                
aGTTGCTGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCa                                                                                                                                               <  2:1205137/148‑1 (MQ=255)
aGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCa                                                                                                                                               <  2:995556/148‑1 (MQ=255)
 gTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAg                                                                                                                                              >  1:2271868/1‑148 (MQ=255)
  ttGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATagcagc                                                                                                                                             >  2:914560/1‑148 (MQ=255)
     aTGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTg                                                                                                                                          >  2:1303484/1‑148 (MQ=255)
        ggCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCgg                                                                                                                                       <  1:1126398/148‑1 (MQ=255)
            aCCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGAt                                                                                                                                   >  2:1576833/1‑148 (MQ=255)
              cACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGAttt                                                                                                                                 >  2:3235785/1‑148 (MQ=255)
              cACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGAttt                                                                                                                                 >  2:3102881/1‑148 (MQ=255)
                cTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATttct                                                                                                                               >  1:2223871/1‑148 (MQ=255)
                 tAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATttctt                                                                                                                              >  2:1895761/1‑148 (MQ=255)
                 tAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATttctt                                                                                                                              >  2:3364540/1‑148 (MQ=255)
                  aaaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATttcttc                                                                                                                             >  1:1858708/1‑148 (MQ=255)
                   aaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCCTTTTTGGCATTTCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTcc                                                                                                                            >  1:2748893/1‑148 (MQ=255)
                     aaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTTGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCgg                                                                                                                          >  1:2104889/1‑148 (MQ=255)
                      aTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGc                                                                                                                         >  1:3173563/1‑148 (MQ=255)
                        gTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGt                                                                                                                       <  1:815971/148‑1 (MQ=255)
                        gTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGt                                                                                                                       <  1:2715201/148‑1 (MQ=255)
                            agagCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTcagc                                                                                                                   >  1:240987/1‑148 (MQ=255)
                            agagCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATACACAACGCAAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTcagc                                                                                                                   >  1:182475/1‑148 (MQ=255)
                               gCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTcagccag                                                                                                                <  1:1286258/148‑1 (MQ=255)
                                ccGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTcagccagc                                                                                                               >  2:2977922/1‑148 (MQ=255)
                                 cGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCa                                                                                                              >  2:1726610/1‑148 (MQ=255)
                                   cAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCaaa                                                                                                            >  1:2747994/1‑148 (MQ=255)
                                     aTTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTTGCATTTCGCGTCATCCACAAAGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATATCAGCGTGCGGTGATTTATTACGGCGTCAGCCAGCAAAca                                                                                                          >  1:3344953/1‑147 (MQ=255)
                                     aTTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACt                                                                                                          >  1:999342/1‑148 (MQ=255)
                                     aTTGCGGCTCTTTTTTCATTCTGCCGTTTAAGTCTCTACGGCTTCATTTTTTGCATTGCGCGTCATCCACAACCCGAGCGCATTAAGCGAGTCTGGAGTGAAAACATATCAGCGGGCCGTGATTTATTACGGCGTCAGCCACCAAAca                                                                                                          >  1:1605547/1‑147 (MQ=255)
                                      ttGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACtt                                                                                                         <  1:1329053/148‑1 (MQ=255)
                                        gCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCa                                                                                                       <  1:1763489/148‑1 (MQ=255)
                                        gCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCa                                                                                                       <  1:3500033/148‑1 (MQ=255)
                                         cGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCAc                                                                                                      >  2:3116089/1‑148 (MQ=255)
                                            ctctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGa                                                                                                   <  2:3454085/148‑1 (MQ=255)
                                            ctctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGa                                                                                                   <  2:2764919/148‑1 (MQ=255)
                                             tctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAc                                                                                                  >  2:1438041/1‑148 (MQ=255)
                                             tctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAc                                                                                                  >  2:1927879/1‑148 (MQ=255)
                                              ctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACt                                                                                                 >  2:2509145/1‑148 (MQ=255)
                                                          tGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGt                                                                                     >  2:2213226/1‑148 (MQ=255)
                                                           gCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTa                                                                                    >  1:1021001/1‑148 (MQ=255)
                                                                  cAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACACATAGCAGAGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTCACTTCATCTTCCTCTAGTGCGaa                                                                             >  1:2431377/1‑148 (MQ=255)
                                                                     tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          >  1:1522229/1‑148 (MQ=255)
                                                                     tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          <  2:248457/148‑1 (MQ=255)
                                                                     tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          >  1:1638452/1‑148 (MQ=255)
                                                                     tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          >  1:638755/1‑148 (MQ=255)
                                                                     tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          >  1:638740/1‑148 (MQ=255)
                                                                     tctcTACGGCTTCATATTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                          >  2:2761295/1‑148 (MQ=255)
                                                                      ctctACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGa                                                                         <  1:1725449/148‑1 (MQ=255)
                                                                         tACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCg                                                                      <  1:355122/148‑1 (MQ=255)
                                                                          aCGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGt                                                                     >  2:1655393/1‑148 (MQ=255)
                                                                              cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgaga                                                                 >  1:1970635/1‑148 (MQ=255)
                                                                              cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgaga                                                                 >  2:1731285/1‑148 (MQ=255)
                                                                                   tttttGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGc                                                            >  2:550603/1‑148 (MQ=255)
                                                                                    ttttGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCa                                                           >  2:1343353/1‑148 (MQ=255)
                                                                                          cATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGaa                                                     >  1:2866920/1‑148 (MQ=255)
                                                                                          cATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGaa                                                     >  1:2162961/1‑148 (MQ=255)
                                                                                            ttGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACa                                                   <  1:570443/148‑1 (MQ=255)
                                                                                              gcgcTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAAt                                                 <  2:639892/148‑1 (MQ=255)
                                                                                               cgcTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATg                                                <  1:1895761/148‑1 (MQ=255)
                                                                                                 cTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATgcg                                              <  2:2626199/148‑1 (MQ=255)
                                                                                                  ttCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATgcgc                                             <  1:2609312/148‑1 (MQ=255)
                                                                                                   tCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGcc                                            <  2:2300011/148‑1 (MQ=255)
                                                                                                   tCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGcc                                            <  1:2903437/148‑1 (MQ=255)
                                                                                                    cATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGccc                                           <  2:1255421/148‑1 (MQ=255)
                                                                                                       ccACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAg                                        >  1:730498/1‑148 (MQ=255)
                                                                                                       ccACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAg                                        >  2:2358465/1‑148 (MQ=255)
                                                                                                       ccACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAg                                        >  1:1277231/1‑148 (MQ=255)
                                                                                                       ccACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAg                                        >  2:2068032/1‑148 (MQ=255)
                                                                                                        cacaACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGa                                       <  2:2845172/148‑1 (MQ=255)
                                                                                                          caACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGaca                                     >  2:444575/1‑148 (MQ=255)
                                                                                                             cGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACgg                                  >  1:1247487/1‑148 (MQ=255)
                                                                                                             cGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACgg                                  >  1:1272534/1‑148 (MQ=255)
                                                                                                              gCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGc                                 <  1:2156538/148‑1 (MQ=255)
                                                                                                               cTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCa                                >  1:1087004/1‑148 (MQ=255)
                                                                                                               cTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCa                                <  2:1857603/148‑1 (MQ=255)
                                                                                                                tAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCaa                               >  2:580144/1‑148 (MQ=255)
                                                                                                                tAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCaa                               >  2:191005/1‑148 (MQ=255)
                                                                                                                tAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCaa                               <  2:2932425/148‑1 (MQ=255)
                                                                                                                tAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCaa                               >  2:1612936/1‑148 (MQ=255)
                                                                                                                  gcgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAtt                             >  2:2083753/1‑148 (MQ=255)
                                                                                                                   cgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                            <  2:2587438/148‑1 (MQ=255)
                                                                                                                   cgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                            >  1:1380190/1‑148 (MQ=255)
                                                                                                                   cgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                            >  1:2551242/1‑148 (MQ=255)
                                                                                                                   cgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                            <  2:828693/148‑1 (MQ=255)
                                                                                                                   cgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                            >  1:3056691/1‑148 (MQ=255)
                                                                                                                     cTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttttt                          <  2:182475/148‑1 (MQ=255)
                                                                                                                     cTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttttt                          <  2:182611/148‑1 (MQ=255)
                                                                                                                     cTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttttt                          <  2:709745/148‑1 (MQ=255)
                                                                                                                      tttCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAtttttc                         >  2:2663376/1‑147 (MQ=255)
                                                                                                                      tttCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTa                         <  1:2149301/148‑1 (MQ=255)
                                                                                                                      tttCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTa                         >  1:2981061/1‑148 (MQ=255)
                                                                                                                        tCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATc                       <  1:1505054/148‑1 (MQ=255)
                                                                                                                           gCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTc                    >  1:3480184/1‑148 (MQ=255)
                                                                                                                            cGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCg                   <  2:2271868/148‑1 (MQ=255)
                                                                                                                            cGAGTCTGGAGTGAACTCATAGCAGCGTGAGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAATTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCTCCCCAGACACGGCAATTTTTATCTTCg                   <  1:2793652/148‑1 (MQ=255)
                                                                                                                                tCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAAt               >  1:1542353/1‑148 (MQ=255)
                                                                                                                                    gAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCACACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGaa           <  2:2953137/148‑1 (MQ=255)
                                                                                                                                    gAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGaa           <  1:623091/148‑1 (MQ=255)
                                                                                                                                     aGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAAc          >  1:687093/1‑148 (MQ=255)
                                                                                                                                     aGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAAc          >  2:3024395/1‑148 (MQ=255)
                                                                                                                                     aGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAAc          <  1:893573/148‑1 (MQ=255)
                                                                                                                                       tGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTg        >  1:3358226/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTG  >  CP009273/2412061‑2412349

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: