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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,732,5490GT14.7% 799.3 / 76.6 258D565Y (GAT→TAT) malSalpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (105/115);  new base T (7/31);  total (112/146)
Fisher's exact test for biased strand distribution p-value = 6.73e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 262 total aligned reads displayed.

TCGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATAT  >  CP009273/3732401‑3732697
                                                                                                                                                    |                                                                                                                                                    
   ccACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATtata                                                                                                                                                   >  1:1539372/1‑149 (MQ=255)
    cACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAc                                                                                                                                                  >  2:746526/1‑149 (MQ=255)
     aCGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                                 <  2:1547039/149‑1 (MQ=255)
              atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTc                                                                                                                                        <  1:2705960/149‑1 (MQ=255)
               tgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTcc                                                                                                                                       <  2:3360461/149‑1 (MQ=255)
                 cgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGt                                                                                                                                     <  1:4009327/149‑1 (MQ=255)
                 cgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGt                                                                                                                                     <  1:3618586/149‑1 (MQ=255)
                  gatgatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTg                                                                                                                                    >  2:1597446/1‑149 (MQ=255)
                   atgatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGa                                                                                                                                   >  2:700999/1‑149 (MQ=255)
                    tgatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGaa                                                                                                                                  >  1:3543500/1‑149 (MQ=255)
                     gatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAg                                                                                                                                 <  1:2934268/149‑1 (MQ=255)
                      atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                                <  2:2627990/149‑1 (MQ=255)
                      atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                                <  2:2477391/149‑1 (MQ=255)
                      atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                                <  2:617257/149‑1 (MQ=255)
                        cAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgggg                                                                                                                              <  1:2601182/149‑1 (MQ=255)
                          aTTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACTTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGc                                                                                                                            >  2:3960233/1‑149 (MQ=255)
                          aTTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGc                                                                                                                            <  2:1654723/149‑1 (MQ=255)
                             tCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGAc                                                                                                                         <  1:1542611/149‑1 (MQ=255)
                              cGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACa                                                                                                                        >  2:1924387/1‑149 (MQ=255)
                               gATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACaa                                                                                                                       >  2:3152146/1‑149 (MQ=255)
                                aTTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACaaa                                                                                                                      <  2:2144700/149‑1 (MQ=255)
                                 ttATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAag                                                                                                                     >  1:543937/1‑149 (MQ=255)
                                 ttATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAag                                                                                                                     >  1:2662844/1‑149 (MQ=255)
                                       ggAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGcagcag                                                                                                               >  2:1806455/1‑149 (MQ=255)
                                         aGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCagag                                                                                                             <  1:2317936/149‑1 (MQ=255)
                                         aGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCagag                                                                                                             <  1:2018295/149‑1 (MQ=255)
                                           cAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtt                                                                                                           <  2:1056734/149‑1 (MQ=255)
                                           cAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtt                                                                                                           <  2:465727/149‑1 (MQ=255)
                                           cAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtt                                                                                                           <  2:566225/149‑1 (MQ=255)
                                            aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtta                                                                                                          >  2:2787180/1‑149 (MQ=255)
                                            aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtta                                                                                                          >  2:1421237/1‑149 (MQ=255)
                                             ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttat                                                                                                         >  2:1216527/1‑149 (MQ=255)
                                             ggcggcGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTCAAGGGGGCGACAAAGCAGCAGAGttat                                                                                                         >  1:44093/1‑149 (MQ=255)
                                              gcggcgAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGttatt                                                                                                        <  2:3610428/149‑1 (MQ=255)
                                                    aaaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                  >  1:2179913/1‑149 (MQ=255)
                                                    aaaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                                  >  1:3920998/1‑149 (MQ=255)
                                                      aGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAg                                                                                                <  2:3550311/149‑1 (MQ=255)
                                                      aGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAg                                                                                                <  2:1502646/149‑1 (MQ=255)
                                                         aGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgcgc                                                                                             >  1:3793273/1‑149 (MQ=255)
                                                           tCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCa                                                                                           <  1:2507302/149‑1 (MQ=255)
                                                                tGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTTAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGgcgc                                                                                      >  1:2344414/1‑149 (MQ=255)
                                                                  tCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCgg                                                                                    >  2:584516/1‑149 (MQ=255)
                                                                      gcgcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACa                                                                                >  1:1732902/1‑149 (MQ=255)
                                                                      gcgcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACa                                                                                >  1:638972/1‑149 (MQ=255)
                                                                        gcAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACaaa                                                                              >  2:1574597/1‑149 (MQ=255)
                                                                           gATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCt                                                                           <  2:2972712/149‑1 (MQ=255)
                                                                            aTGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtt                                                                          <  1:2996645/149‑1 (MQ=255)
                                                                             tGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCttt                                                                         <  2:3148720/149‑1 (MQ=255)
                                                                               gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTa                                                                       >  1:111526/1‑149 (MQ=255)
                                                                                aTACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAt                                                                      <  1:212356/149‑1 (MQ=255)
                                                                                  aCGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATgg                                                                    <  2:3458387/149‑1 (MQ=255)
                                                                                   cGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGt                                                                   <  1:682430/149‑1 (MQ=255)
                                                                                   cGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGt                                                                   <  1:51446/149‑1 (MQ=255)
                                                                                    gACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtg                                                                  <  2:3092002/149‑1 (MQ=255)
                                                                                    gACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtg                                                                  <  1:3162660/149‑1 (MQ=255)
                                                                                     aCCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtga                                                                 <  2:2178686/149‑1 (MQ=255)
                                                                                       cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                               >  2:601297/1‑149 (MQ=255)
                                                                                       cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                               >  2:2667451/1‑149 (MQ=255)
                                                                                       cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                               >  1:2887018/1‑149 (MQ=255)
                                                                                       cTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGtgatg                                                               >  1:3330042/1‑149 (MQ=255)
                                                                                         ggCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                             <  2:2679810/149‑1 (MQ=255)
                                                                                         ggCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                             >  2:1855566/1‑149 (MQ=255)
                                                                                         ggCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                             >  2:3534699/1‑149 (MQ=255)
                                                                                         ggCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                             <  2:2444142/149‑1 (MQ=255)
                                                                                            agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATcc                                                          <  1:559769/149‑1 (MQ=255)
                                                                                            agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATcc                                                          >  1:2578891/1‑149 (MQ=255)
                                                                                             gcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                         <  2:3864774/149‑1 (MQ=255)
                                                                                             gcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                         <  2:4071285/149‑1 (MQ=255)
                                                                                             gcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                         <  2:4071290/149‑1 (MQ=255)
                                                                                              caAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGATGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTc                                                        <  1:2536161/149‑1 (MQ=255)
                                                                                                    gCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTcc                                                  >  2:1948894/1‑149 (MQ=255)
                                                                                                      ggAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGt                                                <  1:2266510/149‑1 (MQ=255)
                                                                                                       gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                               <  2:1748679/149‑1 (MQ=255)
                                                                                                       gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                               <  2:3782539/149‑1 (MQ=255)
                                                                                                       gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                               >  1:1703452/1‑149 (MQ=255)
                                                                                                        agaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTc                                              <  1:1087278/149‑1 (MQ=255)
                                                                                                          aaaTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCgg                                            <  1:1458407/149‑1 (MQ=255)
                                                                                                             ttGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTcc                                         <  1:3462858/149‑1 (MQ=255)
                                                                                                                 aGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACa                                     >  1:3835538/1‑149 (MQ=255)
                                                                                                                  gggTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAg                                    >  1:1819075/1‑149 (MQ=255)
                                                                                                                    gTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGt                                  <  1:3562265/149‑1 (MQ=255)
                                                                                                                      ttCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTc                                <  1:1921847/149‑1 (MQ=255)
                                                                                                                        cAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTg                              <  1:209545/149‑1 (MQ=255)
                                                                                                                         aaCGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGa                             <  1:3024265/149‑1 (MQ=255)
                                                                                                                          aCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                            <  2:930117/149‑1 (MQ=255)
                                                                                                                          aCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                            <  2:2903976/149‑1 (MQ=255)
                                                                                                                             tgtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCg                         >  1:3037632/1‑149 (MQ=255)
                                                                                                                              gtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                        <  1:569828/149‑1 (MQ=255)
                                                                                                                              gtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                        >  1:3068298/1‑149 (MQ=255)
                                                                                                                                 gAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTGCGGTCGGACAGGTTCTGATCCGCTGc                     >  1:2441346/1‑149 (MQ=255)
                                                                                                                                  aGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCa                    >  1:1161649/1‑149 (MQ=255)
                                                                                                                                      aCCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGt                <  2:1178680/149‑1 (MQ=255)
                                                                                                                                       ccTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTa               <  1:1653101/149‑1 (MQ=255)
                                                                                                                                        ctctCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTGATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTac              >  2:2255314/1‑149 (MQ=255)
                                                                                                                                         tctcGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTaca             <  2:753375/149‑1 (MQ=255)
                                                                                                                                          ctcGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTacac            >  2:3499235/1‑149 (MQ=255)
                                                                                                                                              tcgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTc        <  2:3059581/149‑1 (MQ=255)
                                                                                                                                                 cATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGa     >  1:910411/1‑149 (MQ=255)
                                                                                                                                                  aTGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGat    <  1:875347/149‑1 (MQ=255)
                                                                                                                                                    gATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGatat  >  1:659169/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGATAT  >  CP009273/3732401‑3732697

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