Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,353,842 A→T 14.6% C442S (TGC→AGC)  puuA ← glutamate‑‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,353,8420AT14.6% 562.0 / 75.0 199C442S (TGC→AGC) puuAglutamate‑‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base A (76/94);  new base T (17/12);  total (93/106)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 203 total aligned reads displayed.

AATGTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCAAGGCC  >  CP009273/1353696‑1353989
                                                                                                                                                  |                                                                                                                                                   
aaTGTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAg                                                                                                                                                    >  1:3141384/1‑148 (MQ=255)
 aTGTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGa                                                                                                                                                   <  2:3137712/148‑1 (MQ=255)
 aTGTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTATGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGa                                                                                                                                                   >  1:2742259/1‑148 (MQ=255)
   gTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAg                                                                                                                                                 <  1:955478/148‑1 (MQ=255)
     cAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCTAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCg                                                                                                                                               <  1:952856/148‑1 (MQ=255)
       ggCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGtt                                                                                                                                             <  2:2071152/148‑1 (MQ=255)
        gCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTc                                                                                                                                            >  1:1023629/1‑148 (MQ=255)
         ccTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCg                                                                                                                                           <  2:2232668/148‑1 (MQ=255)
          cTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGc                                                                                                                                          >  2:868963/1‑148 (MQ=255)
           tGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGcc                                                                                                                                         >  2:1854091/1‑148 (MQ=255)
               tCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTaaa                                                                                                                                     <  1:2199451/148‑1 (MQ=255)
               tCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTaaa                                                                                                                                     <  2:1228753/148‑1 (MQ=255)
                  gCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAg                                                                                                                                  <  1:2774547/148‑1 (MQ=255)
                   cTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAgc                                                                                                                                 >  2:1845765/1‑148 (MQ=255)
                      aGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAgcggc                                                                                                                              >  2:2888223/1‑148 (MQ=255)
                      aGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAgcggc                                                                                                                              >  2:2028816/1‑148 (MQ=255)
                            aTGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGt                                                                                                                        >  2:1672578/1‑148 (MQ=255)
                             tGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTg                                                                                                                       >  1:2841836/1‑148 (MQ=255)
                               aaaCAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGAt                                                                                                                     <  1:2176115/148‑1 (MQ=255)
                               aaaCAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGAt                                                                                                                     <  2:2828755/148‑1 (MQ=255)
                                 acaacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCCAAATAGCGGCGCAGGTGATCa                                                                                                                   >  2:1144917/1‑148 (MQ=255)
                                  caacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAt                                                                                                                  >  1:403731/1‑148 (MQ=255)
                                  caacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAt                                                                                                                  <  1:1923697/148‑1 (MQ=255)
                                  caacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAt                                                                                                                  >  1:405131/1‑148 (MQ=255)
                                  caacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAt                                                                                                                  >  2:2201888/1‑148 (MQ=255)
                                   aacCCCTCAAGGTGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAtt                                                                                                                 <  1:2940812/148‑1 (MQ=255)
                                   aacCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAtt                                                                                                                 <  1:3014237/148‑1 (MQ=255)
                                    acCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATtc                                                                                                                >  1:1842456/1‑148 (MQ=255)
                                        cGCAAGGGGTCTTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGa                                                                                                            >  1:2067663/1‑148 (MQ=255)
                                        cGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGa                                                                                                            <  2:242164/148‑1 (MQ=255)
                                        cGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGa                                                                                                            <  2:239565/148‑1 (MQ=255)
                                        cGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGa                                                                                                            <  1:1700815/148‑1 (MQ=255)
                                        cGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGa                                                                                                            >  1:696017/1‑148 (MQ=255)
                                            aGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATaaa                                                                                                        >  2:1224028/1‑148 (MQ=255)
                                            aGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATaaa                                                                                                        >  1:746111/1‑148 (MQ=255)
                                              gggTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACt                                                                                                      <  2:2670172/148‑1 (MQ=255)
                                              gggTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACt                                                                                                      >  2:1010858/1‑148 (MQ=255)
                                               ggTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTc                                                                                                     <  1:449770/148‑1 (MQ=255)
                                                  aTTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCAcc                                                                                                  <  1:70024/148‑1 (MQ=255)
                                                  aTTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCAcc                                                                                                  >  1:862063/1‑148 (MQ=255)
                                                    tACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCa                                                                                                <  2:2845911/148‑1 (MQ=255)
                                                    tACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCa                                                                                                <  2:2845785/148‑1 (MQ=255)
                                                        cgTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGc                                                                                            >  2:379960/1‑148 (MQ=255)
                                                          tttttCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCAt                                                                                          <  1:538397/148‑1 (MQ=255)
                                                          tttttCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCAt                                                                                          <  1:1558226/148‑1 (MQ=255)
                                                               cAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGctc                                                                                     <  2:3263716/148‑1 (MQ=255)
                                                                 aCATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTg                                                                                   >  1:2869068/1‑148 (MQ=255)
                                                                 aCATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTg                                                                                   <  2:649792/148‑1 (MQ=255)
                                                                 aCATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTg                                                                                   >  1:453271/1‑148 (MQ=255)
                                                                 aCATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTg                                                                                   >  1:2353299/1‑148 (MQ=255)
                                                                  cATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCAATCTCGATAAAATCACCCAGGGCATCGCTCTgc                                                                                  >  2:1757846/1‑147 (MQ=255)
                                                                    tCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCg                                                                                >  1:263484/1‑148 (MQ=255)
                                                                     ccACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGa                                                                               >  2:2980950/1‑148 (MQ=255)
                                                                     ccACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGa                                                                               >  2:2981616/1‑148 (MQ=255)
                                                                      cACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGaa                                                                              <  1:816700/148‑1 (MQ=255)
                                                                      cACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGaa                                                                              >  1:3098883/1‑148 (MQ=255)
                                                                        cTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATc                                                                            >  1:1930887/1‑148 (MQ=255)
                                                                         tCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCg                                                                           <  2:2377243/148‑1 (MQ=255)
                                                                         tCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCg                                                                           >  2:2482252/1‑148 (MQ=255)
                                                                         tCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCg                                                                           <  2:715871/148‑1 (MQ=255)
                                                                              ttctgtttctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAg                                                                      >  2:2515512/1‑148 (MQ=255)
                                                                              ttctgtttctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAg                                                                      >  2:1266247/1‑148 (MQ=255)
                                                                              ttctgtttctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAg                                                                      >  2:539193/1‑148 (MQ=255)
                                                                               tctgtttctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAgg                                                                     <  2:2139476/148‑1 (MQ=255)
                                                                               tctgtttctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAgg                                                                     <  2:476413/148‑1 (MQ=255)
                                                                                     tctgtGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGc                                                               <  2:1472754/148‑1 (MQ=255)
                                                                                        gtgATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCtt                                                            <  2:683982/148‑1 (MQ=255)
                                                                                         tgatgaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTc                                                           <  2:2734980/148‑1 (MQ=255)
                                                                                         tgatgaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTc                                                           <  1:25932/148‑1 (MQ=255)
                                                                                         tgatgaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTc                                                           <  2:2246669/148‑1 (MQ=255)
                                                                                          gatgaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTcc                                                          <  1:2630180/148‑1 (MQ=255)
                                                                                           atgaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCt                                                         >  2:48427/1‑148 (MQ=255)
                                                                                             gaggagCTAAAACTGCAACAGCTCATCTTTTTTTCAGGCGTGATAGACATGGCAGTATCGTTCGACTTAATAGCGGCGCAGGTTATCATTCTCGTTAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGt                                                       <  2:1346088/148‑1 (MQ=255)
                                                                                             gaGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGt                                                       <  2:746111/148‑1 (MQ=255)
                                                                                              aGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGtt                                                      <  1:1166205/148‑1 (MQ=255)
                                                                                               ggCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGACAGGGTAAGCCTTCCTGTTc                                                     >  1:2288385/1‑148 (MQ=255)
                                                                                                   cTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGc                                                 <  1:2790902/148‑1 (MQ=255)
                                                                                                   cTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGc                                                 <  1:2238946/148‑1 (MQ=255)
                                                                                                      aaccTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCg                                              <  2:1290610/145‑1 (MQ=255)
                                                                                                      aaaCTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCg                                              <  1:379960/148‑1 (MQ=255)
                                                                                                         cTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTg                                           >  1:2505549/1‑148 (MQ=255)
                                                                                                             aaCAGCTCATCATTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCtt                                       <  1:604315/148‑1 (MQ=255)
                                                                                                              aCAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTc                                      >  1:1062357/1‑148 (MQ=255)
                                                                                                               cAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCg                                     <  2:582194/148‑1 (MQ=255)
                                                                                                                    catcatTTTTACAGGCGTGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGActtc                                <  1:1375376/148‑1 (MQ=255)
                                                                                                                    catcatTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGActtc                                >  1:1880297/1‑148 (MQ=255)
                                                                                                                       catTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACttcttc                             <  2:3098883/148‑1 (MQ=255)
                                                                                                                       catTTTGACAGGCGTGAGAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACttcttc                             <  1:2301383/148‑1 (MQ=255)
                                                                                                                               cAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACg                     <  2:2462551/148‑1 (MQ=255)
                                                                                                                                  gCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCCACGGCa                  >  2:1572034/1‑148 (MQ=255)
                                                                                                                                     tGATAGACATGGCTGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCt               >  1:1438550/1‑148 (MQ=255)
                                                                                                                                         aGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGtt           <  1:2921344/148‑1 (MQ=255)
                                                                                                                                         aGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGtt           <  2:817157/148‑1 (MQ=255)
                                                                                                                                          gACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTa          >  2:2479388/1‑148 (MQ=255)
                                                                                                                                          gACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTa          <  2:2556008/148‑1 (MQ=255)
                                                                                                                                             aTGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCa       <  1:2486091/148‑1 (MQ=255)
                                                                                                                                              tGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCaa      <  2:3141384/148‑1 (MQ=255)
                                                                                                                                                 cAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCAAGGc   <  1:104919/148‑1 (MQ=255)
                                                                                                                                                  aGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATACACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGGTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCAAGGcc  <  2:2484863/148‑1 (MQ=255)
                                                                                                                                                  aGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCAAGGcc  <  2:2484928/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATGTCAGGCCTGGCTCCGCTCAGGCCGATGAAACAACCCCGCAAGGGGTATTACGCGTTTTTCAACATCCACTCAATTTCTGTTTCTGTGATGAGGCGCTCAAACTGCAACAGCTCATCATTTTTACAGGCGTGATAGACATGGCAGAAGCGTTCGCCTAAATAGCGGCGCAGGTGATCATTCTCGATAAACTCACCCAGGGCATCGCTCTGGCGAATCGGAAAGGGTAAGCCTTCCTGTTCCAGCCCGTTGCCTTCGACTTCTTCCTGCAACGGCAGCTCGTTATCAAGGCC  >  CP009273/1353696‑1353989

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