Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,252,258 C→A 100% D251Y (GAT→TAT)  tdcG ← L‑serine dehydratase 3, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,252,2580CA99.2% 866.6 / ‑2.2 260D251Y (GAT→TAT) tdcGL‑serine dehydratase 3, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (141/116);  total (143/117)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-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 264 total aligned reads displayed.

ATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCT  >  CP009273/3252112‑3252403
                                                                                                                                                  |                                                                                                                                                 
  aaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTg                                                                                                                                               <  1:3028156/149‑1 (MQ=255)
  aaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTg                                                                                                                                               <  2:3035899/149‑1 (MQ=255)
   aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTgg                                                                                                                                              >  1:715934/1‑149 (MQ=255)
      cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGaaa                                                                                                                                           <  1:2726006/149‑1 (MQ=255)
      cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGaaa                                                                                                                                           <  1:67974/149‑1 (MQ=255)
       aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAAc                                                                                                                                          >  2:3414073/1‑149 (MQ=255)
        gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAAcc                                                                                                                                         <  2:3959999/149‑1 (MQ=255)
        gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAAcc                                                                                                                                         <  1:1995179/149‑1 (MQ=255)
         tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCa                                                                                                                                        <  1:2779127/149‑1 (MQ=255)
          aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAg                                                                                                                                       <  2:2904918/149‑1 (MQ=255)
           cggCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGc                                                                                                                                      <  2:1322667/147‑1 (MQ=255)
           cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGc                                                                                                                                      <  2:4371592/149‑1 (MQ=255)
             gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTg                                                                                                                                    <  2:2406501/149‑1 (MQ=255)
               cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTgac                                                                                                                                  >  1:4272442/1‑149 (MQ=255)
               cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTgac                                                                                                                                  <  1:2564070/149‑1 (MQ=255)
               cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTgac                                                                                                                                  <  1:3895579/149‑1 (MQ=255)
                   ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGGTAGAGATGTTATAGCTGGAAACCAGCTGacgacg                                                                                                                              >  1:3607245/1‑149 (MQ=255)
                     aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCa                                                                                                                            >  2:1222045/1‑149 (MQ=255)
                        gCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTGAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAgcg                                                                                                                         >  1:2536719/1‑149 (MQ=255)
                             acgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacg                                                                                                                    <  2:292266/149‑1 (MQ=255)
                               gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggc                                                                                                                  >  1:2491700/1‑149 (MQ=255)
                                 acCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcac                                                                                                                >  2:3488454/1‑149 (MQ=255)
                                   cGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcacgg                                                                                                              <  2:2827997/149‑1 (MQ=255)
                                    gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcacggc                                                                                                             >  1:1712659/1‑149 (MQ=255)
                                    gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcacggc                                                                                                             >  2:3295685/1‑149 (MQ=255)
                                    gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcacggc                                                                                                             >  2:3690042/1‑149 (MQ=255)
                                      tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGgcgc                                                                                                           >  2:3943911/1‑149 (MQ=255)
                                      tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGgcgc                                                                                                           >  1:2882803/1‑149 (MQ=255)
                                       aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGcgcg                                                                                                          >  2:1538887/1‑149 (MQ=255)
                                           ggTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGcac                                                                                                      <  1:2321913/149‑1 (MQ=255)
                                              gCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACAtt                                                                                                   <  1:3339776/149‑1 (MQ=255)
                                               ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATTAAATTCATCTGATAGTTAGATATTTTATATCTGGAAACCAGCTGAAGACGCAGCGCTAAGGCACGGCGCGGCACATTg                                                                                                  >  1:2650527/1‑149 (MQ=255)
                                                     ccACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGa                                                                                            <  1:1659627/149‑1 (MQ=255)
                                                       aCGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGAcc                                                                                          <  1:3926701/149‑1 (MQ=255)
                                                        cgcgCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCa                                                                                         <  1:469405/149‑1 (MQ=255)
                                                           gcccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGTCATTCATCGGATCGTTAGAGATTTTATATCTTGAAACCAGCTTACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGgc                                                                                      >  1:838133/1‑149 (MQ=38)
                                                           gcccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGgc                                                                                      >  2:3871105/1‑149 (MQ=255)
                                                           gcccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGgc                                                                                      >  2:2137400/1‑149 (MQ=255)
                                                            cccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGgca                                                                                     <  1:4345728/149‑1 (MQ=255)
                                                            cccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGgca                                                                                     <  1:3551095/149‑1 (MQ=255)
                                                                   gCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGcc                                                                              <  2:725079/149‑1 (MQ=255)
                                                                   gCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGcc                                                                              <  1:3582039/149‑1 (MQ=255)
                                                                         ttttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgt                                                                        >  1:2390306/1‑149 (MQ=255)
                                                                         ttttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgt                                                                        <  1:2036042/149‑1 (MQ=255)
                                                                          tttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtg                                                                       >  2:2992388/1‑149 (MQ=255)
                                                                           ttcttcACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtgt                                                                      <  1:3928198/149‑1 (MQ=255)
                                                                             cttcACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTc                                                                    >  1:4397186/1‑149 (MQ=255)
                                                                                cACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCCGCACGCCTTCAGTGTTCATc                                                                 >  1:4434963/1‑149 (MQ=255)
                                                                                 aCTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATcc                                                                <  2:3765556/149‑1 (MQ=255)
                                                                                  cTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATccc                                                               <  1:3764303/149‑1 (MQ=255)
                                                                                  cTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATccc                                                               >  2:61301/1‑149 (MQ=255)
                                                                                   tAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCa                                                              <  2:4219083/149‑1 (MQ=255)
                                                                                    aaCCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCAc                                                             <  1:3815667/149‑1 (MQ=255)
                                                                                      ccGCCAGCGCGTACATGTTGATCCAGTCGATGACATTGATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGt                                                           <  1:730491/149‑1 (MQ=255)
                                                                                       cGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGtt                                                          <  2:631162/149‑1 (MQ=255)
                                                                                       cGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGtt                                                          <  2:1264355/149‑1 (MQ=255)
                                                                                        gCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTc                                                         >  1:2048845/1‑149 (MQ=255)
                                                                                         ccAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCa                                                        <  2:2310309/149‑1 (MQ=255)
                                                                                          cAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCaa                                                       <  1:322921/149‑1 (MQ=255)
                                                                                             cgcgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATAc                                                    <  1:3193254/149‑1 (MQ=255)
                                                                                              gcgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATAcc                                                   >  2:3112096/1‑149 (MQ=255)
                                                                                               cgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCg                                                  >  2:2366298/1‑149 (MQ=255)
                                                                                                gTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGt                                                 >  2:2743458/1‑149 (MQ=255)
                                                                                                gTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGt                                                 >  2:1617801/1‑149 (MQ=255)
                                                                                                 tACATGTTGATCCAGTCGATGACATTCATCGGATCGTTATAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACCTTGAGCGGACCCGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTc                                                >  1:3667862/1‑149 (MQ=255)
                                                                                                 tACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTc                                                <  1:3449689/149‑1 (MQ=255)
                                                                                                     tGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGc                                            <  1:1817308/149‑1 (MQ=255)
                                                                                                         gATCCAGTCGATGACAGTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCa                                        <  2:1088215/149‑1 (MQ=255)
                                                                                                          aTCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCAc                                       >  1:1430405/1‑149 (MQ=255)
                                                                                                             cAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTg                                    <  2:188179/149‑1 (MQ=255)
                                                                                                               gTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGcc                                  <  1:1058247/149‑1 (MQ=255)
                                                                                                                tCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCa                                 <  1:1966541/149‑1 (MQ=255)
                                                                                                                  gATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                               >  1:4413509/1‑149 (MQ=255)
                                                                                                                  gATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                               >  1:793491/1‑149 (MQ=255)
                                                                                                                  gATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                               >  1:793569/1‑149 (MQ=255)
                                                                                                                  gATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                               <  2:2882803/149‑1 (MQ=255)
                                                                                                                   aTGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGAt                              >  2:738990/1‑149 (MQ=255)
                                                                                                                     gACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATAc                            <  1:378220/149‑1 (MQ=255)
                                                                                                                     gACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATAc                            >  2:1930893/1‑149 (MQ=255)
                                                                                                                     gACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATAc                            >  2:1744461/1‑149 (MQ=255)
                                                                                                                      aCATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACg                           >  1:1514011/1‑149 (MQ=255)
                                                                                                                      aCATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACg                           >  1:4047806/1‑149 (MQ=255)
                                                                                                                      aCATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACg                           >  1:4048324/1‑149 (MQ=255)
                                                                                                                        aTTCATCGGATCGTTAGAGATGTTATAGCTGGACACCAGCTGACGACGCAGCGCGACGGCACGGCGCGGCCCAGTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCCCGTTCAAGACCGTCAGGCATCACTTGCCAGATACggg                         >  2:2669597/1‑149 (MQ=255)
                                                                                                                         ttCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGc                        <  2:2232305/149‑1 (MQ=255)
                                                                                                                         ttCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGc                        <  2:965716/149‑1 (MQ=255)
                                                                                                                           cATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCaa                      >  1:636905/1‑149 (MQ=255)
                                                                                                                             tCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCaaac                    <  1:3115614/149‑2 (MQ=255)
                                                                                                                             tCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCaaaa                    >  1:891400/1‑149 (MQ=255)
                                                                                                                               ggATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAAcc                  <  2:2486858/149‑1 (MQ=255)
                                                                                                                               ggATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAAcc                  <  2:93910/149‑1 (MQ=255)
                                                                                                                                     ttAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTc            >  1:1404226/1‑149 (MQ=255)
                                                                                                                                     ttAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTc            <  1:3488454/149‑1 (MQ=255)
                                                                                                                                     ttAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTc            >  2:2542444/1‑149 (MQ=255)
                                                                                                                                        gagaTGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAAt         >  1:609402/1‑149 (MQ=255)
                                                                                                                                        gagaTGTTATAGCTGGAAACCAACTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAAt         <  2:1375127/149‑1 (MQ=255)
                                                                                                                                             gTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCg    >  1:1208677/1‑149 (MQ=255)
                                                                                                                                              ttATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGc   >  2:994237/1‑149 (MQ=255)
                                                                                                                                              ttATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGc   >  1:4166137/1‑149 (MQ=255)
                                                                                                                                               tataGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCt  <  2:2536719/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCT  >  CP009273/3252112‑3252403

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: