Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,512,693 T→A 100% K407* (AAG→TAG)  yjhG ← putative dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,512,6930TA99.5% 616.5 / ‑4.5 185K407* (AAG→TAG) yjhGputative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (76/108);  total (76/109)
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.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 185 total aligned reads displayed.

GCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGG  >  CP009273/4512548‑4512836
                                                                                                                                                 |                                                                                                                                               
 cGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGaa                                                                                                                                              <  1:2367497/148‑1 (MQ=255)
     ggATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGc                                                                                                                                          <  2:2954820/148‑1 (MQ=255)
       aTTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTg                                                                                                                                        <  1:2702748/148‑1 (MQ=255)
       aTTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTg                                                                                                                                        >  2:162081/1‑148 (MQ=255)
       aTTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTg                                                                                                                                        <  2:2323931/148‑1 (MQ=255)
       aTTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTg                                                                                                                                        <  1:2698569/148‑1 (MQ=255)
        tttGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGa                                                                                                                                       <  2:2081186/148‑1 (MQ=255)
             tCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCg                                                                                                                                  >  1:502260/1‑148 (MQ=255)
              cACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGt                                                                                                                                 >  1:1304671/1‑148 (MQ=255)
               aCCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGtt                                                                                                                                <  2:18294/148‑1 (MQ=255)
                ccGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTc                                                                                                                               >  1:1341997/1‑148 (MQ=255)
                  gAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCgg                                                                                                                             <  1:2067595/148‑1 (MQ=255)
                   aaCCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGa                                                                                                                            <  1:1398391/148‑1 (MQ=255)
                   aaCCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGa                                                                                                                            <  1:1458701/148‑1 (MQ=255)
                    aCCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAg                                                                                                                           >  1:39190/1‑148 (MQ=255)
                    aCCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAg                                                                                                                           <  1:2368198/148‑1 (MQ=255)
                      cTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTg                                                                                                                         >  2:305337/1‑148 (MQ=255)
                       ttCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGc                                                                                                                        <  1:808762/148‑1 (MQ=255)
                          tGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTcc                                                                                                                     <  1:2836094/148‑1 (MQ=255)
                            gcgcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCcca                                                                                                                   >  1:2537776/1‑148 (MQ=255)
                            gcgcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCcca                                                                                                                   >  1:3302972/1‑148 (MQ=255)
                             cgcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCccac                                                                                                                  >  2:2186204/1‑148 (MQ=255)
                              gcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCccacc                                                                                                                 <  1:2484379/148‑1 (MQ=255)
                                 atatTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGt                                                                                                              <  1:2206266/148‑1 (MQ=255)
                                      gCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAgg                                                                                                         <  2:1504112/148‑1 (MQ=255)
                                      gCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCATGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAgg                                                                                                         <  1:406069/148‑1 (MQ=255)
                                       cccACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGt                                                                                                        <  2:578715/148‑1 (MQ=255)
                                        ccACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGtt                                                                                                       <  2:954733/148‑1 (MQ=255)
                                        ccACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGtt                                                                                                       <  1:209724/148‑1 (MQ=255)
                                         cACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGttt                                                                                                      <  1:2993355/148‑1 (MQ=255)
                                          aCCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGtttt                                                                                                     <  2:148486/148‑1 (MQ=255)
                                            cGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttct                                                                                                   <  1:2933014/148‑1 (MQ=255)
                                              ggAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCGGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttcttt                                                                                                 >  2:890867/1‑148 (MQ=255)
                                               gAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttctttc                                                                                                >  1:615143/1‑148 (MQ=255)
                                               gAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttctttc                                                                                                >  1:2279393/1‑148 (MQ=255)
                                               gAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttctttc                                                                                                >  2:317201/1‑148 (MQ=255)
                                                  ggTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGc                                                                                             >  1:553555/1‑148 (MQ=255)
                                                      aTAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgc                                                                                         <  2:328323/148‑1 (MQ=255)
                                                      aTAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgc                                                                                         <  2:2533662/148‑1 (MQ=255)
                                                      aTAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgc                                                                                         <  2:2837699/148‑1 (MQ=255)
                                                        aGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgctg                                                                                       >  2:2505557/1‑148 (MQ=255)
                                                        aGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgctg                                                                                       >  1:3192158/1‑148 (MQ=255)
                                                        aGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgctg                                                                                       >  1:297202/1‑148 (MQ=255)
                                                        aGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgctg                                                                                       <  2:493271/148‑1 (MQ=255)
                                                           tGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCg                                                                                    <  2:2127126/148‑1 (MQ=255)
                                                            gAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCgg                                                                                   <  2:1741219/148‑1 (MQ=255)
                                                            gAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCgg                                                                                   <  2:1741028/148‑1 (MQ=255)
                                                             aGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGt                                                                                  <  2:397528/148‑1 (MQ=255)
                                                               gTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTaa                                                                                <  1:1642735/148‑1 (MQ=255)
                                                                   aTCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGt                                                                            >  1:3081754/1‑148 (MQ=255)
                                                                      cgcgcgCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCAt                                                                         <  2:1541916/148‑1 (MQ=255)
                                                                      cgcgagCTTTTGCTTGCTGAGGAGACATGATCACTTCGTCAGCGTTTATTTGTTCCTGATCGAGCAGTAGTTGCTAGATCCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGTTTTTCTTTGAGCGTTCTGCAGGTAACCGTCAt                                                                         <  2:1414023/148‑1 (MQ=255)
                                                                       gcgcgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATa                                                                        <  1:2028626/148‑1 (MQ=255)
                                                                       gcgcgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATa                                                                        <  1:3133716/148‑1 (MQ=255)
                                                                       gcgcgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATa                                                                        <  1:2028613/148‑1 (MQ=255)
                                                                        cgcgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATaa                                                                       >  2:494037/1‑148 (MQ=255)
                                                                           gcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGt                                                                    <  1:2619058/148‑1 (MQ=255)
                                                                           gcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGt                                                                    <  1:27840/148‑1 (MQ=255)
                                                                            cTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTc                                                                   >  1:2683446/1‑148 (MQ=255)
                                                                             ttttGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCt                                                                  <  2:2414225/148‑1 (MQ=255)
                                                                               ttgcttgcTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTc                                                                <  2:2406828/148‑1 (MQ=255)
                                                                                       tgcGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaaca                                                        <  1:2591927/148‑1 (MQ=255)
                                                                                       tgcGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaaca                                                        <  1:1526128/148‑1 (MQ=255)
                                                                                            agaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCg                                                   >  1:2422452/1‑148 (MQ=255)
                                                                                            agaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCg                                                   <  2:1721943/148‑1 (MQ=255)
                                                                                            agaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCg                                                   <  1:274056/148‑1 (MQ=255)
                                                                                              aCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAg                                                 <  2:3225564/148‑1 (MQ=255)
                                                                                              aCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAg                                                 >  1:1071991/1‑148 (MQ=255)
                                                                                              aCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAg                                                 >  1:925863/1‑148 (MQ=255)
                                                                                                 tGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCt                                              >  1:759751/1‑148 (MQ=255)
                                                                                                        ttCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                       >  2:206625/1‑148 (MQ=255)
                                                                                                            tCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCa                                   <  1:898018/148‑1 (MQ=255)
                                                                                                             cAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                  >  1:1514466/1‑148 (MQ=255)
                                                                                                             cAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                  >  1:77835/1‑148 (MQ=255)
                                                                                                             cAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                  >  2:3139618/1‑148 (MQ=255)
                                                                                                             cAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                  <  2:1050846/148‑1 (MQ=255)
                                                                                                              aGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                 <  2:681600/148‑1 (MQ=255)
                                                                                                              aGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                 <  2:2344891/148‑1 (MQ=255)
                                                                                                              aGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                 <  2:2847473/148‑1 (MQ=255)
                                                                                                              aGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGAGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                 <  2:2724694/148‑1 (MQ=255)
                                                                                                               gCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACa                                <  1:2634668/148‑1 (MQ=255)
                                                                                                                 gTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATg                              <  2:307615/148‑1 (MQ=255)
                                                                                                                  ttGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGa                             <  1:1466969/148‑1 (MQ=255)
                                                                                                                  ttGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGa                             >  2:1591227/1‑148 (MQ=255)
                                                                                                                  ttGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGa                             <  2:2709826/148‑1 (MQ=255)
                                                                                                                     aTTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACtt                          <  1:2886011/148‑1 (MQ=255)
                                                                                                                           tCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGca                    <  1:67547/148‑1 (MQ=255)
                                                                                                                               gATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCAcacc                >  2:408139/1‑148 (MQ=255)
                                                                                                                                 tCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCAcaccac              <  2:2674966/148‑1 (MQ=255)
                                                                                                                                    aGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTg           >  2:1693675/1‑148 (MQ=255)
                                                                                                                                    aGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTg           <  2:1438975/148‑1 (MQ=255)
                                                                                                                                     gCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGc          <  2:248306/148‑1 (MQ=255)
                                                                                                                                      cAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGcc         <  2:1384299/148‑1 (MQ=255)
                                                                                                                                      cAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGcc         >  1:1232305/1‑148 (MQ=255)
                                                                                                                                       aGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTTCTGCCGGTAACCGTCATAAAGTCTTCATTCAACAATCAGAGGCTGCGAAGATGCAACATGACTTCCGGCACACCACCTGCCa        >  1:1918815/1‑148 (MQ=255)
                                                                                                                                         gAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATa      <  1:1263536/148‑1 (MQ=255)
                                                                                                                                         gAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATa      >  1:1804796/1‑148 (MQ=255)
                                                                                                                                          aGTTGCTAGAACCGCTGACGCCGTTCGGAGCGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATaa     >  2:124202/1‑148 (MQ=255)
                                                                                                                                             tGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAgg  <  1:3278619/148‑1 (MQ=255)
                                                                                                                                             tGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAgg  <  2:782687/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGG  >  CP009273/4512548‑4512836

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