Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,708,534 T→G 9.3% E49A (GAG→GCG)  yfiE ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,708,5340TG 9.3% 464.3 / 43.0 183E49A (GAG→GCG) yfiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (74/92);  new base G (8/9);  total (82/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 187 total aligned reads displayed.

TCGCCCGAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA  >  CP009273/2708387‑2708681
                                                                                                                                                   |                                                                                                                                                   
tCGCCCGAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTtct                                                                                                                                                     >  1:1312038/1‑148 (MQ=255)
   cccGAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCaa                                                                                                                                                  >  1:1707647/1‑148 (MQ=255)
     cGAAACAACGCGTAGGTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAAt                                                                                                                                                <  1:2421903/148‑1 (MQ=255)
      gAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATa                                                                                                                                               <  1:3284857/148‑1 (MQ=255)
      gAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATa                                                                                                                                               <  1:367876/148‑1 (MQ=255)
           aacGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTgg                                                                                                                                          <  2:893024/148‑1 (MQ=255)
           aacGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTgg                                                                                                                                          <  1:2655118/148‑1 (MQ=255)
           aacGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTgg                                                                                                                                          >  2:1476079/1‑148 (MQ=255)
            acGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGa                                                                                                                                         >  1:2680341/1‑148 (MQ=255)
            acGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGa                                                                                                                                         <  2:1787561/148‑1 (MQ=255)
            acGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGa                                                                                                                                         <  2:1787393/148‑1 (MQ=255)
            acGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGa                                                                                                                                         <  2:165612/148‑1 (MQ=255)
              gcgTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACt                                                                                                                                       <  1:2033812/148‑1 (MQ=255)
               cgTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCCGCCAATTTTCTCAAATAACTGGACTg                                                                                                                                      >  1:163547/1‑148 (MQ=255)
               cgTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGAGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTg                                                                                                                                      >  2:1008796/1‑148 (MQ=255)
                  aGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaaa                                                                                                                                   <  2:1169314/148‑1 (MQ=255)
                   gTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaaaa                                                                                                                                  >  1:2707034/1‑148 (MQ=255)
                   gTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaaaa                                                                                                                                  >  2:1381283/1‑148 (MQ=255)
                      ccccGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTc                                                                                                                               <  2:2986617/148‑1 (MQ=255)
                      ccccGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTc                                                                                                                               >  2:1485173/1‑148 (MQ=255)
                      ccccGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTc                                                                                                                               >  2:2815249/1‑148 (MQ=255)
                        ccGTCCGGATCTGACTCCTTTGTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCt                                                                                                                             <  2:3202425/148‑1 (MQ=255)
                           tCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCt                                                                                                                          >  2:1382270/1‑148 (MQ=255)
                           tCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCt                                                                                                                          >  1:1007291/1‑148 (MQ=255)
                               gATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAg                                                                                                                      <  1:2419333/148‑1 (MQ=255)
                               gATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAg                                                                                                                      <  1:1226374/148‑1 (MQ=255)
                                  cTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAgctg                                                                                                                   <  2:2596082/148‑1 (MQ=255)
                                      cTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGa                                                                                                               <  1:501944/148‑1 (MQ=255)
                                      cTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGa                                                                                                               >  1:2717400/1‑148 (MQ=255)
                                        acTTTTTGGCATCCTCACGAAGCGTATCCATCACTCTGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAat                                                                                                             <  2:2220802/147‑1 (MQ=255)
                                            tttGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATACCTGGACTGAAAACTCCTGCTCAAGCTTCTGAATATGa                                                                                                         >  2:153502/1‑148 (MQ=255)
                                                 cAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGt                                                                                                    >  1:2207928/1‑148 (MQ=255)
                                                  aGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                                                                                                   <  2:1143124/148‑1 (MQ=255)
                                                   gCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCa                                                                                                  >  1:1911487/1‑148 (MQ=255)
                                                          gAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCNGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa                                                                                           <  1:1680557/148‑1 (MQ=255)
                                                            aGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAtt                                                                                         <  1:1968700/148‑1 (MQ=255)
                                                            aGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAtt                                                                                         <  1:1968634/148‑1 (MQ=255)
                                                               gTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg                                                                                      <  1:64022/148‑1 (MQ=255)
                                                                tATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTgtgt                                                                                     <  2:2546793/148‑1 (MQ=255)
                                                                 aTCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTa                                                                                    >  1:2625614/1‑148 (MQ=255)
                                                                  tCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTAc                                                                                   <  1:2278475/148‑1 (MQ=255)
                                                                  tCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTAc                                                                                   <  1:2278478/148‑1 (MQ=255)
                                                                   ccATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACa                                                                                  >  2:2786443/1‑148 (MQ=255)
                                                                   ccATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACa                                                                                  <  2:1414780/148‑1 (MQ=255)
                                                                     aTCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGc                                                                                >  1:2897397/1‑148 (MQ=255)
                                                                     aTCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGc                                                                                >  1:2897286/1‑148 (MQ=255)
                                                                     aTCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGc                                                                                >  1:2806040/1‑148 (MQ=255)
                                                                          tCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAAt                                                                           <  2:17572/148‑1 (MQ=255)
                                                                             ggTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAAtttt                                                                        <  1:1866010/148‑1 (MQ=255)
                                                                               tAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTg                                                                      >  2:2154130/1‑148 (MQ=255)
                                                                               tAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTg                                                                      >  2:1309214/1‑148 (MQ=255)
                                                                                 aGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTTTGTTCAGCACAATTTTTTCg                                                                    >  1:2900110/1‑148 (MQ=255)
                                                                                  gTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGa                                                                   <  1:3252594/148‑1 (MQ=255)
                                                                                   ttCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGaa                                                                  >  2:2318909/1‑148 (MQ=255)
                                                                                   ttCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGaa                                                                  >  1:2188076/1‑148 (MQ=255)
                                                                                   ttCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGaa                                                                  >  1:95501/1‑148 (MQ=255)
                                                                                     cATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGc                                                                <  2:738410/148‑1 (MQ=255)
                                                                                        aaaTGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCg                                                             >  1:3144723/1‑148 (MQ=255)
                                                                                        aaaTGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCg                                                             >  1:3142696/1‑148 (MQ=255)
                                                                                             tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCaaaa                                                        <  2:2707034/148‑1 (MQ=255)
                                                                                                ggCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAgg                                                     <  1:3211989/148‑1 (MQ=255)
                                                                                                 gcagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGa                                                    <  1:2336909/148‑1 (MQ=255)
                                                                                                 gcagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGa                                                    <  1:677823/148‑1 (MQ=255)
                                                                                                   agcaACTTTTTTCCTTCGCGGTTAATGCACATCAGTCGTCCAATTTTCTCAAATAACTTGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAc                                                  <  2:1965286/148‑1 (MQ=255)
                                                                                                     caACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAAccc                                                <  2:437396/148‑1 (MQ=255)
                                                                                                       aCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctc                                              >  2:2109/1‑148 (MQ=255)
                                                                                                       aCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctc                                              <  1:2417367/148‑1 (MQ=255)
                                                                                                        cTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAATGAACCctct                                             >  1:375010/1‑148 (MQ=255)
                                                                                                          tttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCACGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  2:2762241/148‑1 (MQ=255)
                                                                                                          tttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  1:1388823/148‑1 (MQ=255)
                                                                                                          tttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  2:248012/148‑1 (MQ=255)
                                                                                                          tttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  2:1259437/148‑1 (MQ=255)
                                                                                                          tttgtCCTGCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  2:2762245/148‑1 (MQ=255)
                                                                                                          attttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                           <  2:3085299/147‑1 (MQ=255)
                                                                                                                    gcgGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGtt                                 <  2:1312038/148‑1 (MQ=255)
                                                                                                                     cgGGTAAGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAAATTTTTCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGttt                                >  1:203885/1‑148 (MQ=255)
                                                                                                                         tAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAg                            <  2:2417244/148‑1 (MQ=255)
                                                                                                                          aaGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCCAAAAGGAACCCTCTTCCACCACGGTTTTAAGc                           >  1:252394/1‑148 (MQ=255)
                                                                                                                           aGGCACATCCGTCGGCCAATTTTCGCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCg                          >  1:2065801/1‑148 (MQ=255)
                                                                                                                                caTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaata                     >  2:2723022/1‑148 (MQ=255)
                                                                                                                                caTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaata                     >  1:1516749/1‑148 (MQ=255)
                                                                                                                                caTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAGTGTGTACAGCACAAATTTTGCGAAGCTCGCAAAAAGTAACCCTATTCCACCACGGTTTTAAGCGtaata                     >  1:1517157/1‑148 (MQ=255)
                                                                                                                                 aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGtaataa                    <  1:1149163/148‑1 (MQ=255)
                                                                                                                                  tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaa                   <  1:1332813/148‑1 (MQ=255)
                                                                                                                                  tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaa                   >  2:3201868/1‑148 (MQ=255)
                                                                                                                                  tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaa                   <  1:1332868/148‑1 (MQ=255)
                                                                                                                                   ccGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                  <  1:838712/148‑1 (MQ=255)
                                                                                                                                   ccGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAAATTTTGAGAAGCACGCAAAAAGGAACCCGCATCCACCACCGTTTTAAGCGTAaaaaaa                  >  1:2644645/1‑148 (MQ=255)
                                                                                                                                    cGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAAc                 <  2:2200787/148‑1 (MQ=255)
                                                                                                                                      tCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACgg               >  2:3111243/1‑148 (MQ=255)
                                                                                                                                           cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTTTGTACAGCACAATTTTTGCGAAGCTCGCCAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAg          >  1:2093124/1‑148 (MQ=255)
                                                                                                                                            aaTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGACGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGa         >  2:2960727/1‑148 (MQ=255)
                                                                                                                                            aaTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCCAAAAGGAACCCACTTCCCCCACGGTTTTAAGCGTAATAAAACGGCGCAGa         >  2:2258772/1‑148 (MQ=255)
                                                                                                                                            aaTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGa         <  2:170731/148‑1 (MQ=255)
                                                                                                                                             aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTTCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGGTTTAAGCGTAATAAAACGGCGCAGAt        >  1:2492004/1‑148 (MQ=255)
                                                                                                                                             aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAATGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAAAGGCGCAGAt        >  1:473475/1‑148 (MQ=255)
                                                                                                                                             aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGAt        >  1:1738407/1‑148 (MQ=255)
                                                                                                                                               tttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATcc      >  2:3210060/1‑148 (MQ=255)
                                                                                                                                                  ctcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCCAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATa   >  2:471518/1‑148 (MQ=255)
                                                                                                                                                  ctcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATa   <  2:388602/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCCCGAAACAACGCGTAGTTCCCCGTCCGGATCTGACTCCTTTTTGGCAGCCTCACGAAGCGTATCCATCACTCGGGTAAGTTCATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAA  >  CP009273/2708387‑2708681

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: