Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,732,549 G→T 26.1% D565Y (GAT→TAT)  malS → alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,732,5490GT26.1% 732.2 / 202.5 328D565Y (GAT→TAT) malSalpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (107/135);  new base T (39/47);  total (146/182)
Fisher's exact test for biased strand distribution p-value = 9.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

ATCGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATATG  >  CP009273/3732400‑3732698
                                                                                                                                                     |                                                                                                                                                     
  cGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGAt                                                                                                                                                     <  2:4445162/150‑1 (MQ=255)
   gCCACGGCTTCGATGCGATGATCAATGTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACATCTCGTCGCATGATa                                                                                                                                                    >  1:1685698/1‑150 (MQ=255)
    ccACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAc                                                                                                                                                   <  2:4092267/150‑1 (MQ=255)
     cACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                                  <  2:1161758/150‑1 (MQ=255)
     cACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                                  >  1:1459885/1‑150 (MQ=255)
     cACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                                  >  2:2615864/1‑150 (MQ=255)
       cGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCg                                                                                                                                                >  1:5257055/1‑150 (MQ=255)
          cTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCATGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCt                                                                                                                                             >  1:286198/1‑150 (MQ=255)
          cTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCt                                                                                                                                             <  2:382839/150‑1 (MQ=255)
          cTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCt                                                                                                                                             <  2:5231807/150‑1 (MQ=255)
             cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                          <  2:1084602/150‑1 (MQ=255)
             cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                          >  1:4547238/1‑150 (MQ=255)
             cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                          >  2:1490358/1‑150 (MQ=255)
              gatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTc                                                                                                                                         >  2:4926672/1‑150 (MQ=255)
               atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTcc                                                                                                                                        <  1:3930510/150‑1 (MQ=255)
                  cgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTg                                                                                                                                     <  1:5327954/150‑1 (MQ=255)
                  cgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTg                                                                                                                                     >  1:3254955/1‑150 (MQ=255)
                       atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAggg                                                                                                                                <  1:2233421/150‑1 (MQ=255)
                           aTTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATTGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCg                                                                                                                            >  1:895238/1‑150 (MQ=255)
                           aTTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGTGGCg                                                                                                                            <  1:2355929/150‑1 (MQ=255)
                             gtCGATTATCAGGAGCATGCGGCGAAAGAAGTCTCCTGTATTGCTCTTTTGGTTACGTCCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGAc                                                                                                                          <  2:2869701/149‑1 (MQ=255)
                               cGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACaa                                                                                                                        <  2:1885044/150‑1 (MQ=255)
                                 aTTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAag                                                                                                                      >  1:261444/1‑150 (MQ=255)
                                   tATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagca                                                                                                                    >  1:1449788/1‑150 (MQ=255)
                                   tATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGAGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagca                                                                                                                    <  1:3034484/150‑1 (MQ=255)
                                         gAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCagag                                                                                                              <  1:3825430/150‑1 (MQ=255)
                                          aGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGt                                                                                                             >  2:3745206/1‑150 (MQ=255)
                                             aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttat                                                                                                          >  2:3255917/1‑150 (MQ=255)
                                             aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttat                                                                                                          <  1:5018142/150‑1 (MQ=255)
                                              ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttatt                                                                                                         <  2:2113332/150‑1 (MQ=255)
                                              ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttatt                                                                                                         >  1:3145984/1‑150 (MQ=255)
                                              ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttatt                                                                                                         >  1:2778392/1‑150 (MQ=255)
                                              ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttatt                                                                                                         <  2:4510053/150‑1 (MQ=255)
                                                   cgAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAt                                                                                                    <  2:2669906/150‑1 (MQ=255)
                                                   cgAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAt                                                                                                    <  2:869785/150‑1 (MQ=255)
                                                    gAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                   <  2:4020990/150‑1 (MQ=255)
                                                    gAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                   <  1:179520/150‑1 (MQ=255)
                                                    gAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                   <  1:343213/150‑1 (MQ=255)
                                                     aaaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTa                                                                                                  >  1:1808145/1‑150 (MQ=255)
                                                      aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATGACTATTAg                                                                                                 >  2:4747222/1‑150 (MQ=255)
                                                       aGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGTCAGCAAATGGCGGAGAAATTGCAGGGATTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgc                                                                                                >  2:1766068/1‑150 (MQ=255)
                                                          aGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGcc                                                                                             <  2:5310859/150‑1 (MQ=255)
                                                          aGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGcc                                                                                             <  2:1358366/150‑1 (MQ=255)
                                                          aGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGcc                                                                                             <  1:2080126/150‑1 (MQ=255)
                                                             cGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAgg                                                                                          >  2:3369981/1‑150 (MQ=255)
                                                             cGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAgg                                                                                          <  2:5034534/150‑1 (MQ=255)
                                                               aCTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGgcg                                                                                        <  1:1909103/150‑1 (MQ=255)
                                                                 tGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGcgcg                                                                                      <  1:4428319/150‑1 (MQ=255)
                                                                   tCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGATATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGt                                                                                    <  2:3961734/150‑1 (MQ=255)
                                                                     tGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTAc                                                                                  <  1:3352590/150‑1 (MQ=255)
                                                                        cgcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACaaa                                                                               <  1:391773/150‑1 (MQ=255)
                                                                           aGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCt                                                                            >  1:3991458/1‑150 (MQ=255)
                                                                             aTGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCttt                                                                          <  1:3865561/150‑1 (MQ=255)
                                                                              tGGTTACGACCTGGCAGCAAATGGCGGTGATTTTGCAGGGTTTCAACGTGTTGATCTACCTCTCGTCGCATTATACCCGCCTGGTCCGTGAAGTTTGATACAAAGCAGCAGAGTTTTTACTATTAGCGCAATTCGCTGTTCAAATAtttt                                                                         <  2:614218/150‑1 (MQ=255)
                                                                                gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAt                                                                       <  1:1205337/150‑1 (MQ=255)
                                                                                gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAt                                                                       >  1:2444157/1‑150 (MQ=255)
                                                                                 aTACGACCTGGCAGCAAATGGCGGAGAAATTGCTGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGTGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATg                                                                      <  1:3682044/150‑1 (MQ=255)
                                                                                   aCGAACTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGt                                                                    <  2:4579208/150‑1 (MQ=255)
                                                                                       ccTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                                <  1:5074426/150‑1 (MQ=255)
                                                                                       ccTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                                >  2:1965633/1‑150 (MQ=255)
                                                                                        cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatga                                                               <  1:1250375/150‑1 (MQ=255)
                                                                                         tGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                              <  1:355644/150‑1 (MQ=255)
                                                                                          ggCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAAt                                                             <  2:2508866/150‑1 (MQ=255)
                                                                                             agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                          <  2:2973253/150‑1 (MQ=255)
                                                                                               caAATGGCGGAGATATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTcg                                                        <  2:5274922/150‑1 (MQ=255)
                                                                                                aaaTGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCACATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTcgc                                                       >  2:2096287/1‑150 (MQ=255)
                                                                                                 aaTGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTcgcg                                                      >  1:1622603/1‑150 (MQ=255)
                                                                                                    ggCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTcc                                                   <  1:1782000/150‑1 (MQ=255)
                                                                                                    ggCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTcc                                                   >  2:1027710/1‑150 (MQ=255)
                                                                                                     gCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCg                                                  >  2:3490831/1‑150 (MQ=255)
                                                                                                      cGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTccat                                                 <  2:1459835/150‑3 (MQ=255)
                                                                                                       ggAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                                <  2:3663066/150‑1 (MQ=255)
                                                                                                        gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTc                                               <  1:4595530/150‑1 (MQ=255)
                                                                                                          gaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCgg                                             >  1:4885423/1‑150 (MQ=255)
                                                                                                            aaTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTc                                           >  2:3881012/1‑150 (MQ=255)
                                                                                                             aTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTcc                                          <  1:5214925/150‑1 (MQ=255)
                                                                                                                  aGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAg                                     <  1:3008463/150‑1 (MQ=255)
                                                                                                                     gTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                  <  1:4366955/150‑1 (MQ=255)
                                                                                                                     gTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                  <  2:4804438/150‑1 (MQ=255)
                                                                                                                     gTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                                  >  2:4135327/1‑150 (MQ=255)
                                                                                                                      tttCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTc                                 >  1:646636/1‑150 (MQ=255)
                                                                                                                       ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCt                                <  1:143639/150‑1 (MQ=255)
                                                                                                                        tCAACGTGTTGAGCGACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTg                               <  2:4672014/150‑1 (MQ=255)
                                                                                                                         cAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGa                              >  2:4163164/1‑150 (MQ=255)
                                                                                                                         cAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGa                              <  2:3141977/150‑1 (MQ=255)
                                                                                                                          aaCGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                             >  1:1338948/1‑150 (MQ=255)
                                                                                                                            cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATcc                           >  1:4752613/1‑150 (MQ=255)
                                                                                                                              tgtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                         >  1:1035385/1‑150 (MQ=255)
                                                                                                                                ttGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTg                       <  1:3615930/150‑1 (MQ=255)
                                                                                                                                ttGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTg                       <  1:4424807/150‑1 (MQ=255)
                                                                                                                                 tGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGc                      >  2:1540493/1‑150 (MQ=255)
                                                                                                                                   aGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTCCAGGTTCTGATCCGCTGCaa                    >  1:40506/1‑150 (MQ=255)
                                                                                                                                   aGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCaa                    >  2:2439055/1‑150 (MQ=255)
                                                                                                                                         ctctCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTaca              >  2:758261/1‑150 (MQ=255)
                                                                                                                                         ctctCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTaca              <  2:4336208/150‑1 (MQ=255)
                                                                                                                                           ctcGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACg            >  2:1636839/1‑150 (MQ=255)
                                                                                                                                           ctcGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACg            >  1:1992297/1‑150 (MQ=255)
                                                                                                                                             cgtcgCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGtt          <  1:2677930/150‑1 (MQ=255)
                                                                                                                                                  cATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGat     >  2:144364/1‑150 (MQ=255)
                                                                                                                                                     gATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATATg  <  1:2925327/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATATG  >  CP009273/3732400‑3732698

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