Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,680,642 C→A 100% E68D (GAG→GAT fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,680,6420CA96.8% 583.6 / 9.0 190E68D (GAG→GATfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (91/93);  total (94/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 192 total aligned reads displayed.

ATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGC  >  CP009273/1680497‑1680786
                                                                                                                                                 |                                                                                                                                                
  tCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCtt                                                                                                                                              <  1:799312/148‑1 (MQ=255)
  tCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCtt                                                                                                                                              <  1:2613603/148‑1 (MQ=255)
   cACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTc                                                                                                                                             <  2:1185176/148‑1 (MQ=255)
    aCTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCa                                                                                                                                            <  1:1281981/148‑1 (MQ=255)
    aCTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCa                                                                                                                                            <  1:1282596/148‑1 (MQ=255)
    aCTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCCGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCa                                                                                                                                            <  1:563091/148‑1 (MQ=255)
     cTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAg                                                                                                                                           <  1:3515544/148‑1 (MQ=255)
       ggCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGac                                                                                                                                         >  1:1095993/1‑148 (MQ=255)
        gCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGaca                                                                                                                                        >  1:2409364/1‑148 (MQ=255)
        gCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGaca                                                                                                                                        <  2:1471978/148‑1 (MQ=255)
        gCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGaca                                                                                                                                        <  2:690950/148‑1 (MQ=255)
         cccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaa                                                                                                                                       <  2:3042352/148‑1 (MQ=255)
         cccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaa                                                                                                                                       >  1:1923676/1‑148 (MQ=255)
           cGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                     <  1:719099/148‑1 (MQ=255)
           cGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGacaaca                                                                                                                                     <  2:2036016/148‑1 (MQ=255)
              ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGc                                                                                                                                  <  1:2874152/148‑1 (MQ=255)
               tAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGcc                                                                                                                                 <  1:691239/148‑1 (MQ=255)
               tAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGcc                                                                                                                                 >  2:3380561/1‑148 (MQ=255)
               tAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGcc                                                                                                                                 >  2:1613340/1‑148 (MQ=255)
                aGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGATTTATCTTCAGACAACAAGCCt                                                                                                                                <  1:3162954/148‑1 (MQ=255)
                 gCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTa                                                                                                                               >  1:3114880/1‑148 (MQ=255)
                  ccAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTaa                                                                                                                              >  2:1718800/1‑148 (MQ=255)
                    aGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAAt                                                                                                                            <  2:773630/148‑1 (MQ=255)
                     gCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATc                                                                                                                           >  1:3261971/1‑148 (MQ=255)
                      cACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCt                                                                                                                          >  1:261057/1‑148 (MQ=255)
                      cACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCt                                                                                                                          >  1:2061394/1‑148 (MQ=255)
                       aCTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCtt                                                                                                                         <  2:3273083/148‑1 (MQ=255)
                        cTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTc                                                                                                                        <  1:3250449/148‑1 (MQ=255)
                        cTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTc                                                                                                                        >  1:1042468/1‑148 (MQ=255)
                          tCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCAt                                                                                                                      <  2:2826945/148‑1 (MQ=255)
                          tCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCAt                                                                                                                      <  2:2262723/148‑1 (MQ=255)
                            gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTa                                                                                                                    <  2:3165976/148‑1 (MQ=255)
                              tcatgttcatgtTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAAc                                                                                                                  >  1:2555713/1‑148 (MQ=255)
                                 tgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACttt                                                                                                               <  2:2262469/148‑1 (MQ=255)
                                  gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                              >  2:2178677/1‑148 (MQ=255)
                                   ttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTg                                                                                                             <  1:2199185/148‑1 (MQ=255)
                                   ttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTg                                                                                                             >  1:3026699/1‑148 (MQ=255)
                                   ttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTg                                                                                                             >  1:1252492/1‑148 (MQ=255)
                                         ttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTg                                                                                                       >  1:1702400/1‑148 (MQ=255)
                                         ttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAAATTTTGCCGCTg                                                                                                       >  1:2124827/1‑148 (MQ=255)
                                          tACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGc                                                                                                      >  1:232603/1‑148 (MQ=255)
                                          tACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGc                                                                                                      >  1:3280711/1‑148 (MQ=255)
                                              ttGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGc                                                                                                  <  2:839826/148‑1 (MQ=255)
                                                  gTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTgg                                                                                              <  2:3213216/148‑1 (MQ=255)
                                                      cGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAg                                                                                          <  2:3317362/148‑1 (MQ=255)
                                                        gAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcg                                                                                        >  1:548694/1‑148 (MQ=255)
                                                         aGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgcgc                                                                                       <  1:49998/148‑1 (MQ=255)
                                                           ccGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCa                                                                                     <  2:2201429/148‑1 (MQ=255)
                                                               tCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAgcgc                                                                                 >  2:401854/1‑148 (MQ=255)
                                                               tCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAgcgc                                                                                 <  1:1392473/148‑1 (MQ=255)
                                                                  gCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATg                                                                              >  1:3170494/1‑148 (MQ=255)
                                                                      gATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATc                                                                          >  1:1659440/1‑148 (MQ=255)
                                                                       aTAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCa                                                                         >  1:1546735/1‑148 (MQ=255)
                                                                       aTAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCa                                                                         >  1:249821/1‑148 (MQ=255)
                                                                           ccAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGa                                                                     <  1:1749697/148‑1 (MQ=255)
                                                                             aGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAgg                                                                   <  1:2344219/148‑1 (MQ=255)
                                                                              gCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGt                                                                  <  1:197122/148‑1 (MQ=255)
                                                                                  gAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGc                                                              <  2:1090170/148‑1 (MQ=255)
                                                                                  gAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGc                                                              >  1:881837/1‑148 (MQ=255)
                                                                                     ttCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtt                                                           >  1:1459079/1‑148 (MQ=255)
                                                                                       cgtcgtCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCTCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtttt                                                         <  2:2095521/148‑1 (MQ=255)
                                                                                         tcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                       >  1:3549611/1‑148 (MQ=255)
                                                                                          cgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                      >  2:2106131/1‑148 (MQ=255)
                                                                                          cgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                      >  1:1978965/1‑148 (MQ=255)
                                                                                          cgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                      <  1:3124344/148‑1 (MQ=255)
                                                                                            tcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCg                                                    >  2:1459625/1‑148 (MQ=255)
                                                                                               tGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCg                                                 <  2:1401207/148‑1 (MQ=255)
                                                                                                 cTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCTCTTTATCTTCAGACAACAAGCCTAAATCTTCATTATCTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                               <  1:41088/148‑1 (MQ=255)
                                                                                                    tCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCCGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATg                                            >  1:1402894/1‑148 (MQ=255)
                                                                                                    tCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATg                                            <  1:2806130/148‑1 (MQ=255)
                                                                                                     ccTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGc                                           >  1:564841/1‑148 (MQ=255)
                                                                                                       tGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCgg                                         <  1:3339697/148‑1 (MQ=255)
                                                                                                         ccAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGaa                                       <  2:1633947/148‑1 (MQ=255)
                                                                                                         ccAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGaa                                       <  2:261728/148‑1 (MQ=255)
                                                                                                              aCTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCt                                  >  2:3554907/1‑148 (MQ=255)
                                                                                                               cTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTc                                 <  1:2632951/148‑1 (MQ=255)
                                                                                                               cTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTc                                 >  1:2383866/1‑148 (MQ=255)
                                                                                                                ttCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                >  1:931012/1‑148 (MQ=255)
                                                                                                                ttCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTACACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTcc                                >  2:1443750/1‑148 (MQ=255)
                                                                                                                 tCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCa                               >  2:482796/1‑148 (MQ=255)
                                                                                                                     ccgccgCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcga                           >  2:2916165/1‑148 (MQ=255)
                                                                                                                        ccgccTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcgagcg                        >  2:1611295/1‑148 (MQ=255)
                                                                                                                         cgccTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGt                       <  2:1274637/148‑1 (MQ=255)
                                                                                                                           ccggACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                     <  2:483098/145‑1 (MQ=255)
                                                                                                                            cTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGa                    <  2:1095993/148‑1 (MQ=255)
                                                                                                                             tGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAg                   >  1:3238791/1‑148 (MQ=255)
                                                                                                                              gACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGt                  >  2:2417133/1‑148 (MQ=255)
                                                                                                                              gACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGt                  >  1:419689/1‑148 (MQ=255)
                                                                                                                               aCGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGtt                 <  2:1947377/148‑1 (MQ=255)
                                                                                                                               aCGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGtt                 <  2:145368/148‑1 (MQ=255)
                                                                                                                                 gAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtg               <  2:2457489/148‑1 (MQ=255)
                                                                                                                                 gAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtg               <  2:2930957/148‑1 (MQ=255)
                                                                                                                                 gAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtg               <  1:2922457/148‑1 (MQ=255)
                                                                                                                                  aaTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt              <  2:1339210/148‑1 (MQ=255)
                                                                                                                                   aTGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgtg             <  1:1057778/148‑1 (MQ=255)
                                                                                                                                      gcgcTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgcgc          >  2:1553970/1‑148 (MQ=255)
                                                                                                                                         ctcgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGcccc       >  2:2025216/1‑148 (MQ=255)
                                                                                                                                         ctcgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGcccc       >  1:1577858/1‑148 (MQ=255)
                                                                                                                                          tcgctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCca      >  2:96390/1‑148 (MQ=255)
                                                                                                                                              tttATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGc  <  2:1140072/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGC  >  CP009273/1680497‑1680786

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: