Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,252,258 C→A 67.3% 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,2580CA67.3% 341.0 / 305.5 283D251Y (GAT→TAT) tdcGL‑serine dehydratase 3, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (40/53);  new base A (94/96);  total (134/149)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

GTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGC  >  CP009273/3252109‑3252407
                                                                                                                                                     |                                                                                                                                                     
gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATc                                                                                                                                                       <  2:3988357/150‑1 (MQ=255)
  aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGc                                                                                                                                                     <  2:638369/150‑1 (MQ=255)
   ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCt                                                                                                                                                    >  2:2125794/1‑150 (MQ=255)
      aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGa                                                                                                                                                 <  1:554998/150‑1 (MQ=255)
        ccAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCCTCGGCTCGTTAGAGATGTTATCGCTGGtaa                                                                                                                                               >  1:2464336/1‑150 (MQ=255)
         cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGaaaa                                                                                                                                              >  1:4680719/1‑149 (MQ=255)
         cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAAc                                                                                                                                              >  1:2128479/1‑150 (MQ=255)
         cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAAc                                                                                                                                              >  1:3515859/1‑150 (MQ=255)
         cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAAc                                                                                                                                              >  1:3515692/1‑150 (MQ=255)
          aGTACTGCCTGAATAATTCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAAcc                                                                                                                                             <  1:2183113/150‑1 (MQ=255)
            tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAg                                                                                                                                           <  1:4279306/150‑1 (MQ=255)
             aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGc                                                                                                                                          >  2:4427280/1‑150 (MQ=255)
             aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGc                                                                                                                                          <  1:2272630/150‑1 (MQ=255)
              cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCt                                                                                                                                         <  1:4463099/150‑1 (MQ=255)
               tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTg                                                                                                                                        >  1:1290432/1‑150 (MQ=255)
                     aataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGacgacg                                                                                                                                  <  2:1927358/150‑1 (MQ=255)
                     aataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGacgacg                                                                                                                                  <  2:3343742/150‑1 (MQ=255)
                        aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAg                                                                                                                               >  1:1096874/1‑150 (MQ=255)
                            ccGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCt                                                                                                                           <  2:2586132/150‑1 (MQ=255)
                               cacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacg                                                                                                                        >  2:17245/1‑150 (MQ=255)
                               cacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacg                                                                                                                        >  2:3604749/1‑150 (MQ=255)
                               cacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCAACGGATCGGTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCAacg                                                                                                                        >  1:1201080/1‑150 (MQ=255)
                                  gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggca                                                                                                                     <  2:114837/150‑1 (MQ=255)
                                      tgTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTacggcacggc                                                                                                                 >  1:3080539/2‑150 (MQ=255)
                                       gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGgcg                                                                                                                >  1:4665065/1‑150 (MQ=255)
                                         tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGcgcg                                                                                                              >  1:4422919/1‑150 (MQ=255)
                                         tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGcgcg                                                                                                              <  1:2384318/150‑1 (MQ=255)
                                         tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGcgcg                                                                                                              <  1:2184520/150‑1 (MQ=255)
                                         tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGAGTTTTATTCACTAACCGACAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATTTTATAGATGGAAAGCAGCTGACGACGCAGCGCTACGGCACGGcgcg                                                                                                              <  1:1137658/150‑1 (MQ=38)
                                          aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCgg                                                                                                             >  2:4511519/1‑150 (MQ=255)
                                             cGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGcac                                                                                                          <  1:2993987/150‑1 (MQ=255)
                                                    gTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGc                                                                                                   <  1:2599516/150‑1 (MQ=255)
                                                    gTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGc                                                                                                   <  1:856243/150‑1 (MQ=255)
                                                        ccACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGAc                                                                                               <  1:1368175/150‑1 (MQ=255)
                                                                  gccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCAc                                                                                     <  1:2690918/150‑1 (MQ=255)
                                                                     aGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGcc                                                                                  <  1:696426/150‑1 (MQ=255)
                                                                      gCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCt                                                                                 <  2:4080295/150‑1 (MQ=255)
                                                                       cTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCtt                                                                                <  1:4755960/150‑1 (MQ=255)
                                                                       cTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCtt                                                                                >  2:2253613/1‑150 (MQ=255)
                                                                            ttttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtg                                                                           <  1:1406525/150‑1 (MQ=255)
                                                                            ttttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtg                                                                           >  2:1577487/1‑150 (MQ=255)
                                                                            ttttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtg                                                                           >  1:1010687/1‑150 (MQ=255)
                                                                             tttCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAgtgt                                                                          <  2:2571624/150‑1 (MQ=255)
                                                                              ttcttcACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGtt                                                                         >  1:1591263/1‑150 (MQ=255)
                                                                                cttcACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCa                                                                       >  2:2410612/1‑150 (MQ=255)
                                                                                 ttcACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCAt                                                                      <  2:3780840/150‑1 (MQ=255)
                                                                                   cACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATcc                                                                    <  1:333657/150‑1 (MQ=255)
                                                                                      tAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCAc                                                                 >  2:3084484/1‑150 (MQ=255)
                                                                                       aaCCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACg                                                                >  1:1435284/1‑150 (MQ=255)
                                                                                       aaCCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACg                                                                >  2:2037965/1‑150 (MQ=255)
                                                                                        aCCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCGGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGt                                                               >  1:3678273/1‑150 (MQ=255)
                                                                                              aGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATa                                                         >  1:4370051/1‑150 (MQ=255)
                                                                                              aGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATa                                                         >  1:3130896/1‑150 (MQ=255)
                                                                                                cgcgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATAcc                                                       <  1:1340692/150‑1 (MQ=255)
                                                                                                  cgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGt                                                     >  2:3105360/1‑150 (MQ=255)
                                                                                                  cgTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGt                                                     <  2:3481301/150‑1 (MQ=255)
                                                                                                    tACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCa                                                   <  2:4689918/150‑1 (MQ=255)
                                                                                                    tACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCa                                                   <  2:2239025/150‑1 (MQ=255)
                                                                                                     aCATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCAt                                                  <  1:244300/150‑1 (MQ=255)
                                                                                                     aCATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCAt                                                  <  1:3812082/150‑1 (MQ=255)
                                                                                                      cATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACAATGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATg                                                 >  1:3488286/1‑150 (MQ=255)
                                                                                                          ttGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATc                                             <  1:2907989/150‑1 (MQ=255)
                                                                                                           tGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCa                                            <  1:1212412/150‑1 (MQ=255)
                                                                                                           tGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCa                                            <  1:1994876/150‑1 (MQ=255)
                                                                                                           tGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCa                                            <  1:2812497/150‑1 (MQ=255)
                                                                                                            gATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCAc                                           <  1:1353124/150‑1 (MQ=255)
                                                                                                            gATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCAc                                           >  2:397288/1‑150 (MQ=255)
                                                                                                             agccAGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACt                                          <  2:2767044/148‑1 (MQ=255)
                                                                                                                 aGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGcc                                      <  1:1091605/150‑1 (MQ=255)
                                                                                                                 aGTCGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGcc                                      >  1:834094/1‑150 (MQ=255)
                                                                                                                    cGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                                   >  2:731657/1‑150 (MQ=255)
                                                                                                                    cGATGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGa                                   >  2:3036474/1‑150 (MQ=255)
                                                                                                                      aTGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATa                                 >  1:3234670/1‑150 (MQ=255)
                                                                                                                      aTGACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATa                                 <  1:3552504/150‑1 (MQ=255)
                                                                                                                        gACATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACg                               >  1:1790896/1‑150 (MQ=255)
                                                                                                                          cATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACggg                             >  1:3816742/1‑150 (MQ=255)
                                                                                                                          cATTCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACggg                             >  2:2922716/1‑150 (MQ=255)
                                                                                                                             tCATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCaa                          <  1:1944619/150‑1 (MQ=255)
                                                                                                                              cATCGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCaaa                         >  1:432603/1‑150 (MQ=255)
                                                                                                                                 cGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAAcc                      >  2:1212327/1‑150 (MQ=255)
                                                                                                                                 cGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAAcc                      <  2:496756/150‑1 (MQ=255)
                                                                                                                                 cGGATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAAcc                      <  1:1026321/150‑1 (MQ=255)
                                                                                                                                  tgATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCg                     >  2:4231396/2‑150 (MQ=255)
                                                                                                                                  ggATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCg                     >  2:793152/1‑150 (MQ=255)
                                                                                                                                  ggATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCg                     >  2:2143415/1‑150 (MQ=255)
                                                                                                                                   gATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCgg                    >  2:1991093/1‑150 (MQ=255)
                                                                                                                                   gATCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCgg                    >  2:4301892/1‑150 (MQ=255)
                                                                                                                                     tCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCg                  <  1:2733664/150‑1 (MQ=255)
                                                                                                                                     tCGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCg                  <  1:3851858/150‑1 (MQ=255)
                                                                                                                                      cGTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGt                 <  1:1703201/150‑1 (MQ=255)
                                                                                                                                       gTTAGAGATGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTc                >  1:4566044/1‑150 (MQ=255)
                                                                                                                                        ttAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCa               >  2:4411530/1‑150 (MQ=255)
                                                                                                                                           gagaTGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAAtt            <  2:1591272/150‑1 (MQ=255)
                                                                                                                                               tGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCg        >  2:1259299/1‑150 (MQ=255)
                                                                                                                                               tGTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCg        <  2:2700775/150‑1 (MQ=255)
                                                                                                                                                gTTATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGc       >  1:4041147/1‑150 (MQ=255)
                                                                                                                                                 ttATAGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCt      >  1:581326/1‑150 (MQ=255)
                                                                                                                                                  tataGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCtt     <  2:15958/150‑1 (MQ=255)
                                                                                                                                                   aTCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCttt    <  1:765164/150‑1 (MQ=255)
                                                                                                                                                   aTCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCttt    <  1:3202182/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGGATCGTTAGAGATGTTATCGCTGGAAACCAGCTGACGACGCAGCGCTACGGCACGGCGCGGCACATTGAGCGGACCAGGCAGCACGCCTTCAGTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGC  >  CP009273/3252109‑3252407

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: