Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,234 C→T 32.9% 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,2340CT32.9% 308.2 / 214.8 213E27K (GAA→AAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (61/82);  new base T (26/44);  total (87/126)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 213 total aligned reads displayed.

CAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTG  >  CP009273/647088‑647381
                                                                                                                                                  |                                                                                                                                                   
 aGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGAAACGCATTTAATAATATTTACGGCAATTCCACCGCACGCGATAAGCTATTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGc                                                                                                                                                 >  2:4283735/1‑150 (MQ=255)
  gcgcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCa                                                                                                                                                <  2:2158428/150‑1 (MQ=255)
    gcCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAgg                                                                                                                                              <  2:391515/150‑1 (MQ=255)
       aGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAAt                                                                                                                                           >  2:3372603/1‑150 (MQ=255)
        gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                          >  1:3610491/1‑150 (MQ=255)
        gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAAtt                                                                                                                                          >  1:2540115/1‑150 (MQ=255)
        gTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAAtt                                                                                                                                          <  1:729603/150‑1 (MQ=255)
         tCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTg                                                                                                                                         >  2:3806897/1‑150 (MQ=255)
           cTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTggg                                                                                                                                       <  1:1562197/150‑1 (MQ=255)
              cACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGa                                                                                                                                    >  2:749284/1‑150 (MQ=255)
               aCCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGAt                                                                                                                                   >  2:2575041/1‑150 (MQ=255)
                ccGCGGACTGATTCACTGATAGCAACGAATTTAATAATATTTCCGGAAATTCCACCGCACGCGATAAGCTTTTCATCGCGGTTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGAtt                                                                                                                                  <  1:3006329/150‑1 (MQ=255)
                ccGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGAtc                                                                                                                                  <  2:1195735/150‑2 (MQ=255)
                ccGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGAtc                                                                                                                                  <  2:936295/150‑2 (MQ=255)
                 cgcgGACGGATTCACTGATTGCAACGCTTTTAATAATATTTCCGGAAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAACTACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAttt                                                                                                                                 <  1:1386387/150‑1 (MQ=255)
                 cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAttt                                                                                                                                 <  2:3379615/150‑1 (MQ=255)
                 cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGAttt                                                                                                                                 <  1:1386346/150‑1 (MQ=255)
                 cgcgGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGAttt                                                                                                                                 <  2:2112464/150‑1 (MQ=255)
                     gacggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGc                                                                                                                             >  2:2343155/1‑150 (MQ=255)
                      acggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGcc                                                                                                                            >  2:2730230/1‑150 (MQ=255)
                        ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                          >  1:441502/1‑150 (MQ=255)
                        ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                          >  2:768425/1‑150 (MQ=255)
                        ggaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAt                                                                                                                          >  1:2845759/1‑150 (MQ=255)
                         gaTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGACTTGGGCGATTTCCGCCAtt                                                                                                                         >  1:4208085/1‑150 (MQ=255)
                          aTTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCAttt                                                                                                                        >  1:3435065/1‑150 (MQ=255)
                               cTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTa                                                                                                                   <  1:2008841/150‑1 (MQ=255)
                                tGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAt                                                                                                                  >  2:2384193/1‑150 (MQ=255)
                                tGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAt                                                                                                                  <  2:116401/150‑1 (MQ=255)
                                   tAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                               <  1:1225198/150‑1 (MQ=255)
                                   tAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                               <  1:2170220/150‑1 (MQ=255)
                                   tAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                               <  2:51075/150‑1 (MQ=255)
                                   tAGCAACGCATTTAATAATATTTCCGGCAAGTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAtttt                                                                                                               <  2:1506651/150‑1 (MQ=255)
                                    aGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTAttttg                                                                                                              <  2:4127152/150‑2 (MQ=255)
                                    aGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTc                                                                                                              <  2:681747/150‑1 (MQ=255)
                                       aaCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGa                                                                                                           >  2:377167/1‑150 (MQ=255)
                                       aaCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGa                                                                                                           >  2:3990155/1‑150 (MQ=255)
                                        aCGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGaa                                                                                                          <  2:1794690/150‑1 (MQ=255)
                                         cGCATTTAATAATATTTCCGGCAATTCCACAGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAAc                                                                                                         <  2:2930477/150‑1 (MQ=255)
                                          gCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACg                                                                                                        >  2:1684307/1‑150 (MQ=255)
                                            aTTTAATAATATGTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtt                                                                                                      <  2:3540428/150‑1 (MQ=255)
                                             tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGttt                                                                                                     >  1:3021358/1‑150 (MQ=255)
                                             tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGttt                                                                                                     >  1:2127823/1‑150 (MQ=255)
                                             tttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGttt                                                                                                     <  2:1152680/150‑1 (MQ=255)
                                              ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                    <  2:3056040/150‑1 (MQ=255)
                                              ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                    <  2:877226/150‑1 (MQ=255)
                                              ttAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGtttt                                                                                                    <  1:3863461/150‑1 (MQ=255)
                                               taataaTATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTa                                                                                                   <  1:1025378/150‑1 (MQ=255)
                                               taataaTATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTa                                                                                                   <  1:2870706/150‑1 (MQ=255)
                                                aataatATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCACCATTTTTTTATTTTCTGAACGTTTTAc                                                                                                  <  1:847175/150‑1 (MQ=255)
                                                 ataataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcg                                                                                                 >  1:928465/1‑150 (MQ=255)
                                                   aataTTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTAcgcg                                                                                               >  2:2400790/1‑150 (MQ=255)
                                                     tatTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCggg                                                                                             <  2:1034069/150‑1 (MQ=255)
                                                      atTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATATTTTTATTTTCTGAAAGTTTTACGCGGGt                                                                                            >  1:2555976/1‑150 (MQ=255)
                                                          ccGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGaaa                                                                                        >  2:293911/1‑150 (MQ=255)
                                                           cGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACACGGGTGaaaa                                                                                       >  2:4035525/1‑150 (MQ=255)
                                                              cAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                    >  1:883279/1‑150 (MQ=255)
                                                              cAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                    >  1:604055/1‑150 (MQ=255)
                                                              cAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAatat                                                                                    >  1:3265811/1‑150 (MQ=255)
                                                                 ggccACCGCCCGCGATCAGCTTTTCCTCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAACTCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTGTCTGAACGTTTTACGCGGGTGAAAATATCAt                                                                                 <  1:1161336/148‑1 (MQ=255)
                                                                  tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                                >  1:3061594/1‑150 (MQ=255)
                                                                  tCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAtt                                                                                <  1:438558/150‑1 (MQ=255)
                                                                     aCCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGcc                                                                             <  1:1740248/150‑1 (MQ=255)
                                                                      ccGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCg                                                                            <  1:2478766/150‑1 (MQ=255)
                                                                       cGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGa                                                                           <  2:3174341/150‑1 (MQ=255)
                                                                       cGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGa                                                                           >  2:2378033/1‑150 (MQ=255)
                                                                           cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                       <  1:1787968/150‑1 (MQ=255)
                                                                           cgcgATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                       >  2:2513426/1‑150 (MQ=255)
                                                                               aTAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataata                                                                   <  1:1539139/150‑1 (MQ=255)
                                                                                 aaGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                                 >  1:2904904/1‑150 (MQ=255)
                                                                                   gCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                               <  2:105909/150‑1 (MQ=255)
                                                                                     ttttCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                             <  1:398400/150‑1 (MQ=255)
                                                                                         cATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATcc                                                         >  2:2339976/1‑150 (MQ=255)
                                                                                         cATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATcc                                                         <  1:2232618/150‑1 (MQ=255)
                                                                                              cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                    <  1:3615264/150‑1 (MQ=255)
                                                                                              cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                                    <  1:2575485/150‑1 (MQ=255)
                                                                                               gggTTACGTTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAgg                                                   <  2:1179682/150‑1 (MQ=255)
                                                                                               gggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAgg                                                   >  1:2430890/1‑150 (MQ=255)
                                                                                               gggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCGGAAgg                                                   >  1:2062482/1‑150 (MQ=255)
                                                                                                ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGt                                                  <  2:3837141/150‑1 (MQ=255)
                                                                                                 gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                 <  2:2540115/150‑1 (MQ=255)
                                                                                                 gTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTg                                                 <  2:3610491/150‑1 (MQ=255)
                                                                                                     cGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtg                                             <  2:1382718/150‑1 (MQ=255)
                                                                                                       gTAATAAATACTTCGACTGTGGTGTCAACGCTCAACTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgtt                                           <  2:1429822/150‑1 (MQ=255)
                                                                                                         aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                         <  2:3220702/150‑1 (MQ=255)
                                                                                                         aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                         <  2:1192625/150‑1 (MQ=255)
                                                                                                          ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATgttgtt                                        <  1:2384193/150‑1 (MQ=255)
                                                                                                                aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                  <  2:2801179/150‑1 (MQ=255)
                                                                                                                aCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                                  <  1:1247671/150‑1 (MQ=255)
                                                                                                                  ttCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAttt                                <  1:2343033/150‑1 (MQ=255)
                                                                                                                  ttCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAttt                                <  1:3372603/150‑1 (MQ=255)
                                                                                                                    cGACTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGc                              <  1:116548/150‑1 (MQ=255)
                                                                                                                    cGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGc                              >  1:3542882/1‑150 (MQ=255)
                                                                                                                      aCTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCaa                            >  2:236916/1‑150 (MQ=255)
                                                                                                                       cTGTGGTGTCAACGCTCAAATCATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAAc                           <  1:1454113/150‑1 (MQ=255)
                                                                                                                        tgtgGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                          <  2:2371694/150‑1 (MQ=255)
                                                                                                                                cAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt                  <  2:902884/150‑1 (MQ=255)
                                                                                                                                      tCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg            >  2:2587777/1‑150 (MQ=255)
                                                                                                                                        aaaTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAATGTGAATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGtt          >  2:1543199/1‑150 (MQ=255)
                                                                                                                                            cATTTTTATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaa      <  2:1940424/150‑1 (MQ=255)
                                                                                                                                            cATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCaa      >  1:2384540/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGCGCCAGTCCTTCACCGCGGACGGATTCACTGATAGCAACGCATTTAATAATATTTCCGGCAATTCCACCGCACGCGATAAGCTTTTCATCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTG  >  CP009273/647088‑647381

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: