Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,234 C→T 100% E27K (GAA→AAA)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,2340CT96.0% 800.2 / 17.8 250E27K (GAA→AAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base T (119/120);  total (122/128)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 251 total aligned reads displayed.

GCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAA  >  CP009273/647090‑647378
                                                                                                                                                |                                                                                                                                                
    cAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAgg                                                                                                                                           >  2:2735459/1‑148 (MQ=255)
      gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGaa                                                                                                                                         <  2:2116579/148‑1 (MQ=255)
        ccTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                       <  1:1673210/148‑1 (MQ=255)
         cTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTg                                                                                                                                      >  1:607962/1‑148 (MQ=255)
           tCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAAGTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTggg                                                                                                                                    <  1:2152706/148‑1 (MQ=255)
              ccGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGa                                                                                                                                 >  1:3314483/1‑148 (MQ=255)
               cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAt                                                                                                                                <  2:1276846/148‑1 (MQ=255)
               cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAt                                                                                                                                <  2:3051095/148‑1 (MQ=255)
                 cggacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAttt                                                                                                                              >  1:1691016/1‑148 (MQ=255)
                    acggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCg                                                                                                                           >  1:3511855/1‑148 (MQ=255)
                      ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGcc                                                                                                                         <  2:1090096/148‑1 (MQ=255)
                       gaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCa                                                                                                                        >  1:877768/1‑148 (MQ=255)
                       gaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCa                                                                                                                        >  1:3401578/1‑148 (MQ=255)
                        aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                       >  2:2509313/1‑148 (MQ=255)
                            aCTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttttt                                                                                                                   >  2:655944/1‑148 (MQ=255)
                             cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                  >  2:1619702/1‑148 (MQ=255)
                             cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAtttttt                                                                                                                  >  2:53019/1‑148 (MQ=255)
                              tGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCAttttttt                                                                                                                 <  2:3212030/148‑1 (MQ=255)
                               gATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTa                                                                                                                <  1:1584075/148‑1 (MQ=255)
                                 tAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtt                                                                                                              >  2:1742105/1‑148 (MQ=255)
                                   gCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                            >  1:2223175/1‑148 (MQ=255)
                                    cAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTc                                                                                                           >  1:1609778/1‑148 (MQ=255)
                                     aaCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCt                                                                                                          <  1:890347/148‑1 (MQ=255)
                                      aCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTg                                                                                                         >  1:1171083/1‑148 (MQ=255)
                                       cGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGa                                                                                                        <  2:2429683/148‑1 (MQ=255)
                                          aTTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                     <  2:1629547/148‑1 (MQ=255)
                                           tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGt                                                                                                    >  1:1741618/1‑148 (MQ=255)
                                            ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                   <  1:2692517/148‑1 (MQ=255)
                                            ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                   >  2:422225/1‑148 (MQ=255)
                                            ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                   >  1:3217628/1‑148 (MQ=255)
                                            ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                   <  1:3086541/148‑1 (MQ=255)
                                               ataataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTa                                                                                                <  1:797059/148‑1 (MQ=255)
                                                 aataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcg                                                                                              <  1:2887679/148‑1 (MQ=255)
                                                   tatTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcgcg                                                                                            <  2:1027139/148‑1 (MQ=255)
                                                     tttCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCggg                                                                                          <  1:1032347/148‑1 (MQ=255)
                                                      gtCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGt                                                                                         <  2:1363762/147‑1 (MQ=255)
                                                           gCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGaaaa                                                                                    >  1:169489/1‑148 (MQ=255)
                                                              aTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                     >  2:2301668/1‑128 (MQ=255)
                                                              aTTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                 >  1:2722916/1‑148 (MQ=255)
                                                                    ccGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGc                                                                           >  1:2609790/1‑148 (MQ=255)
                                                                     cGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                          >  1:2269217/1‑148 (MQ=255)
                                                                     cGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                          <  2:2836570/148‑1 (MQ=255)
                                                                        aCGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGaa                                                                       <  1:3125831/148‑1 (MQ=255)
                                                                         cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAAc                                                                      <  1:1349306/148‑1 (MQ=255)
                                                                          gcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACa                                                                     <  1:226738/148‑1 (MQ=255)
                                                                           cgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                    <  2:1536839/148‑1 (MQ=255)
                                                                            gATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACata                                                                   <  1:1864046/148‑1 (MQ=255)
                                                                             aTAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                                  >  2:1860503/1‑148 (MQ=255)
                                                                             aTAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                                  <  1:991387/148‑1 (MQ=255)
                                                                              tAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                                 >  2:327690/1‑148 (MQ=255)
                                                                              tAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                                 >  1:101145/1‑148 (MQ=255)
                                                                               aaGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataata                                                                >  1:3195019/1‑148 (MQ=255)
                                                                                aGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACAtaataa                                                               >  1:1639683/1‑148 (MQ=255)
                                                                                 gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                              >  2:969220/1‑148 (MQ=255)
                                                                                 gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                              <  1:3075449/148‑1 (MQ=255)
                                                                                 gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                              <  1:301701/148‑1 (MQ=255)
                                                                                 gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                              <  1:774202/148‑1 (MQ=255)
                                                                                  cTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaat                                                             >  2:1100370/1‑148 (MQ=255)
                                                                                  cTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaat                                                             <  2:450268/148‑1 (MQ=255)
                                                                                     ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTTGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                          <  2:2610945/148‑1 (MQ=255)
                                                                                     ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                          <  1:2959977/148‑1 (MQ=255)
                                                                                     ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                          <  1:1223473/148‑1 (MQ=255)
                                                                                     ttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                          >  1:2576405/1‑148 (MQ=255)
                                                                                       cATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                        <  1:3301036/148‑1 (MQ=255)
                                                                                        aTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATc                                                       <  1:3559016/148‑1 (MQ=255)
                                                                                         tCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATcc                                                      >  2:2367408/1‑148 (MQ=255)
                                                                                         tCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATcc                                                      <  1:2330616/148‑1 (MQ=255)
                                                                                          cgcgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                     >  2:619607/1‑148 (MQ=255)
                                                                                            cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGa                                                   <  1:2033813/148‑1 (MQ=255)
                                                                                              ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                 >  1:2534896/1‑148 (MQ=255)
                                                                                              ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                 <  2:3401578/148‑1 (MQ=255)
                                                                                              ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                 <  2:1271448/148‑1 (MQ=255)
                                                                                               gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAgg                                                <  2:3035374/148‑1 (MQ=255)
                                                                                                 tACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                              <  2:3232312/148‑1 (MQ=255)
                                                                                                      taataaATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgt                                         >  2:1236433/1‑148 (MQ=255)
                                                                                                       aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgtt                                        >  1:671718/1‑148 (MQ=255)
                                                                                                        ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                       >  2:1175861/1‑148 (MQ=255)
                                                                                                        ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                       <  1:195026/148‑1 (MQ=255)
                                                                                                         taaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                      >  2:158430/1‑148 (MQ=255)
                                                                                                         taaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                      >  2:3118588/1‑148 (MQ=255)
                                                                                                          aaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATgttgtt                                     >  1:1329382/1‑148 (MQ=255)
                                                                                                           aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                 <  2:2090627/135‑1 (MQ=255)
                                                                                                           aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTa                                    >  2:2318951/1‑148 (MQ=255)
                                                                                                            ataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTAt                                   <  2:1299234/148‑1 (MQ=255)
                                                                                                             taCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                  >  1:2804434/1‑148 (MQ=255)
                                                                                                             taCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                  >  1:942820/1‑148 (MQ=255)
                                                                                                             taCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                                  >  1:2426149/1‑148 (MQ=255)
                                                                                                             taCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCATGAAGGTGCATGTTGTTATc                                  >  2:2352341/1‑148 (MQ=255)
                                                                                                               cTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGa                                <  2:746584/148‑1 (MQ=255)
                                                                                                                     cTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCa                          <  1:1742105/148‑1 (MQ=255)
                                                                                                                         ggtgTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGa                      <  1:2751184/148‑1 (MQ=255)
                                                                                                                             tCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgt                  <  2:928881/148‑1 (MQ=255)
                                                                                                                                aCGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt               >  1:1887864/1‑148 (MQ=255)
                                                                                                                                aCGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt               <  2:1928649/148‑1 (MQ=255)
                                                                                                                                 cGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATgttgtt              <  2:2426149/148‑1 (MQ=255)
                                                                                                                                    tCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCaa           <  1:2686678/148‑1 (MQ=255)
                                                                                                                                       aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGt        <  2:1809314/148‑1 (MQ=255)
                                                                                                                                       aaTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGt        <  2:867440/148‑1 (MQ=255)
                                                                                                                                         tCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTg      >  2:3318108/1‑148 (MQ=255)
                                                                                                                                             ttttATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaaa  <  1:2187348/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAA  >  CP009273/647090‑647378

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