Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,318,829 G→T 100% E159* (GAA→TAA)  yciO → putative RNA binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,318,8290GT96.8% 715.2 / 11.3 234E159* (GAA→TAA) yciOputative RNA binding protein
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (112/114);  total (115/119)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

TCTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTAT  >  CP009273/1318684‑1318973
                                                                                                                                                 |                                                                                                                                                
tCTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGgaataa                                                                                                                                                <  1:2995670/148‑1 (MQ=255)
tCTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGgaagaa                                                                                                                                                >  1:116129/1‑148 (MQ=255)
tCTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGATGCCAGGCAGCGAATTTACCGAATCGGACCCGgaataa                                                                                                                                                >  2:2979751/1‑148 (MQ=255)
 cTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATaaa                                                                                                                                               >  1:2216151/1‑148 (MQ=255)
 cTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATaaa                                                                                                                                               <  2:277841/148‑1 (MQ=255)
 cTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATaaa                                                                                                                                               <  1:909553/148‑1 (MQ=255)
 cTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATaaa                                                                                                                                               <  2:79154/148‑1 (MQ=255)
       cAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAg                                                                                                                                         <  2:1789934/148‑1 (MQ=255)
         ggAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGAt                                                                                                                                       <  1:2615041/148‑1 (MQ=255)
         ggAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGAt                                                                                                                                       >  2:2497870/1‑148 (MQ=255)
          gAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATc                                                                                                                                      >  2:206293/1‑148 (MQ=255)
            aaaaCGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGt                                                                                                                                    <  2:3086090/148‑1 (MQ=255)
             aaaCGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTc                                                                                                                                   >  1:2928484/1‑148 (MQ=255)
              aaCGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCt                                                                                                                                  >  2:829015/1‑148 (MQ=255)
              aaCGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCt                                                                                                                                  <  2:1474175/148‑1 (MQ=255)
               aCGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTg                                                                                                                                 >  1:3252856/1‑148 (MQ=255)
                     aaCCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGaaaaa                                                                                                                           >  1:320022/1‑148 (MQ=255)
                       ccATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGACTAAATTAAAGATCGTCTGGAacaaca                                                                                                                         >  2:2090371/1‑148 (MQ=255)
                       ccATCGGTATGCGTGTGCCTTCTAAACCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACa                                                                                                                         <  1:2066869/148‑1 (MQ=255)
                         aTCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGaaaaa                                                                                                                           <  2:2086462/144‑1 (MQ=255)
                         aTCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGaaaaa                                                                                                                           >  1:2086764/1‑144 (MQ=255)
                         aTCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGAACGTCTGGAAAAACAgg                                                                                                                       >  2:178486/1‑148 (MQ=255)
                             gTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGa                                                                                                                   >  1:3480361/1‑148 (MQ=255)
                             gTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGa                                                                                                                   >  1:1491953/1‑148 (MQ=255)
                              tATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGAt                                                                                                                  >  2:2766248/1‑148 (MQ=255)
                                tGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGAttt                                                                                                                >  1:2243609/1‑148 (MQ=255)
                                tGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGAttt                                                                                                                <  1:1746533/148‑1 (MQ=255)
                                 gCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTg                                                                                                               >  2:2752020/1‑148 (MQ=255)
                                 gCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTg                                                                                                               >  1:1440752/1‑148 (MQ=255)
                                  cGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGa                                                                                                              <  1:2766248/148‑1 (MQ=255)
                                      tgCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGattat                                                                                                          <  2:2050895/148‑1 (MQ=255)
                                       gCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGattatt                                                                                                         <  2:1569922/148‑1 (MQ=255)
                                         cTTCTAACCCTATATCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCAGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCa                                                                                                       <  2:1374085/148‑1 (MQ=255)
                                          ttCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAc                                                                                                      >  1:3210162/1‑148 (MQ=255)
                                           tCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcg                                                                                                     >  2:3301588/1‑148 (MQ=255)
                                            cTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcgg                                                                                                    <  2:2274779/148‑1 (MQ=255)
                                            cTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcgg                                                                                                    <  2:513882/148‑1 (MQ=255)
                                               aCCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcggcgg                                                                                                 <  2:1645105/148‑1 (MQ=255)
                                               aCCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcggcgg                                                                                                 <  1:1154073/148‑1 (MQ=255)
                                               aCCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCAcggcgg                                                                                                 <  1:311670/148‑1 (MQ=255)
                                                   tATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTAt                                                                                             >  2:3193730/1‑148 (MQ=255)
                                                   tATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTAt                                                                                             <  2:822208/148‑1 (MQ=255)
                                                    aTCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTAtc                                                                                            <  1:2350512/148‑1 (MQ=255)
                                                      cGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTAtctc                                                                                          >  1:2472404/1‑148 (MQ=255)
                                                         tCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCggg                                                                                       <  1:1339059/148‑1 (MQ=255)
                                                          cAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGc                                                                                      <  2:2141294/148‑1 (MQ=255)
                                                          cAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGc                                                                                      <  1:708233/148‑1 (MQ=255)
                                                           aGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCa                                                                                     >  2:663295/1‑148 (MQ=255)
                                                             gCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGa                                                                                   >  1:2220215/1‑148 (MQ=255)
                                                                ttACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAAc                                                                                >  2:1936617/1‑148 (MQ=255)
                                                                    ttGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGAc                                                                            <  1:1359513/148‑1 (MQ=255)
                                                                          cATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGt                                                                      >  1:577673/1‑148 (MQ=255)
                                                                                  gAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATc                                                              <  2:3322901/148‑1 (MQ=255)
                                                                                    aCCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCtt                                                            <  1:350274/148‑1 (MQ=255)
                                                                                     ccGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTa                                                           >  2:564304/1‑148 (MQ=255)
                                                                                     ccGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTa                                                           >  1:2854267/1‑148 (MQ=255)
                                                                                     ccGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTa                                                           >  1:2012135/1‑148 (MQ=255)
                                                                                       gATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTAcc                                                         >  2:353983/1‑148 (MQ=255)
                                                                                            tttCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                                    >  2:1073449/1‑148 (MQ=255)
                                                                                            tttCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                                    >  1:121302/1‑148 (MQ=255)
                                                                                            tttCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCgatga                                                    >  1:121191/1‑148 (MQ=255)
                                                                                              tCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGATTAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGTGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATa                                                  <  2:2730174/148‑1 (MQ=255)
                                                                                                cACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGGTATTGATCTTACCGATGATACg                                                >  1:2140217/1‑148 (MQ=255)
                                                                                                 aCTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGc                                               <  1:537586/148‑1 (MQ=255)
                                                                                                   ttCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCg                                             >  2:705077/1‑148 (MQ=255)
                                                                                                   ttCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCg                                             >  1:1032086/1‑148 (MQ=255)
                                                                                                     cGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGt                                           <  1:1419810/148‑1 (MQ=255)
                                                                                                       cTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCg                                         >  1:2700837/1‑148 (MQ=255)
                                                                                                        tAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgt                                        <  2:2339626/148‑1 (MQ=255)
                                                                                                         aaTGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtg                                       <  2:785192/148‑1 (MQ=255)
                                                                                                            gctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtggtg                                    >  1:988811/1‑148 (MQ=255)
                                                                                                            gctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtggtg                                    <  2:1199864/148‑1 (MQ=255)
                                                                                                            gctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtggtg                                    >  1:1573686/1‑148 (MQ=255)
                                                                                                            gctgcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCCATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCgtggtg                                    >  1:1573425/1‑148 (MQ=255)
                                                                                                               gcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGt                                 >  2:232220/1‑148 (MQ=255)
                                                                                                               gcCAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGt                                 >  2:1413795/1‑148 (MQ=255)
                                                                                                                 cAGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGa                               <  2:714387/148‑1 (MQ=255)
                                                                                                                  aGGCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGaa                              >  1:2749395/1‑148 (MQ=255)
                                                                                                                   ggCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAg                             <  2:2366622/148‑1 (MQ=255)
                                                                                                                   ggCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAg                             >  2:3226656/1‑148 (MQ=255)
                                                                                                                    gCAGCGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAgg                            >  2:3301781/1‑148 (MQ=255)
                                                                                                                        cGAATTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTAGTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAGGGCGTa                        >  1:1861259/1‑148 (MQ=255)
                                                                                                                          aaTTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAgg                      <  1:612385/148‑1 (MQ=255)
                                                                                                                           aTTTACCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGt                     >  2:1706123/1‑148 (MQ=255)
                                                                                                                               aCCGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtg                 <  2:2971354/148‑1 (MQ=255)
                                                                                                                                 cGAATCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGAtgtg               <  2:1854902/148‑1 (MQ=255)
                                                                                                                                    aTCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAg            <  2:3474041/148‑1 (MQ=255)
                                                                                                                                     tCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGc           >  2:2232654/1‑148 (MQ=255)
                                                                                                                                     tCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGc           >  1:806044/1‑148 (MQ=255)
                                                                                                                                     tCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGc           <  1:3301588/148‑1 (MQ=255)
                                                                                                                                     tCGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGc           >  2:467304/1‑148 (MQ=255)
                                                                                                                                      cGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGcc          >  1:2526331/1‑148 (MQ=255)
                                                                                                                                      cGGACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGcc          >  1:2526414/1‑148 (MQ=255)
                                                                                                                                       ggACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCt         <  1:2584063/148‑1 (MQ=255)
                                                                                                                                        gACCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCtt        >  2:256681/1‑148 (MQ=255)
                                                                                                                                         aCCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCttt       <  2:2659851/148‑1 (MQ=255)
                                                                                                                                         aCCCGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCttt       >  2:1796228/1‑148 (MQ=255)
                                                                                                                                            cGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCtt    <  2:2928484/148‑1 (MQ=255)
                                                                                                                                            cGGAATAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCtt    >  2:2644788/1‑148 (MQ=255)
                                                                                                                                              gaataaATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTtat  >  2:2451800/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCTGTTGCAGGAAAAACGCAAAACCATCGGTATGCGTGTGCCTTCTAACCCTATCGCTCAGGCGTTACTTGAAGCATTGGGCGAACCGATGCTTTCCACTTCGCTAATGCTGCCAGGCAGCGAATTTACCGAATCGGACCCGGAAGAAATTAAAGATCGTCTGGAAAAACAGGTAGATTTGATTATTCACGGCGGCTATCTCGGGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTAT  >  CP009273/1318684‑1318973

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: