Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,087,754 C→A 41.7% H31N (CAT→AAT)  rhaS → transcriptional activator of rhaBAD and rhaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,087,7540CA41.7% 315.3 / inf 280H31N (CAT→AAT) rhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (64/99);  new base A (59/58);  total (123/157)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

AAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGT  >  CP009273/4087608‑4087901
                                                                                                                                                  |                                                                                                                                                   
   aaCCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatc                                                                                                                                               <  1:3627534/150‑1 (MQ=255)
       tGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGa                                                                                                                                           >  1:3673522/1‑150 (MQ=255)
           cTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGAtttt                                                                                                                                       <  2:3947708/150‑1 (MQ=255)
             cGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCa                                                                                                                                     <  1:2300008/150‑1 (MQ=255)
               tGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATg                                                                                                                                   <  1:895255/150‑1 (MQ=255)
               tGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATg                                                                                                                                   >  1:2538198/1‑150 (MQ=255)
               tGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATg                                                                                                                                   >  2:4327375/1‑150 (MQ=255)
                gATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGa                                                                                                                                  <  2:3035578/150‑1 (MQ=255)
                gATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGa                                                                                                                                  >  2:1134276/1‑150 (MQ=255)
                gATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTACGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGa                                                                                                                                  >  2:443999/1‑150 (MQ=255)
                 aTTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGaa                                                                                                                                 >  2:656376/1‑150 (MQ=255)
                   tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGTCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAt                                                                                                                               <  2:3782473/150‑1 (MQ=255)
                   tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAt                                                                                                                               >  1:1528238/1‑150 (MQ=255)
                    aCTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAtt                                                                                                                              <  2:3943533/150‑1 (MQ=255)
                     cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg                                                                                                                             <  2:3171389/150‑1 (MQ=255)
                     cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg                                                                                                                             <  2:2015572/150‑1 (MQ=255)
                     cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg                                                                                                                             <  2:399918/150‑1 (MQ=255)
                     cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATtg                                                                                                                             <  1:42448/150‑1 (MQ=255)
                     cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATtg                                                                                                                             <  1:2019810/150‑1 (MQ=255)
                       aTTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtgtg                                                                                                                           <  1:3461019/150‑1 (MQ=255)
                       aTTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtgtg                                                                                                                           <  1:1747229/150‑1 (MQ=255)
                       aTTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATtgtg                                                                                                                           <  2:2899473/150‑1 (MQ=255)
                       aTTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATtgtg                                                                                                                           <  2:1538536/150‑1 (MQ=255)
                         ttCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGAt                                                                                                                         <  2:848207/150‑1 (MQ=255)
                          tCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGAtt                                                                                                                        <  2:87926/150‑1 (MQ=255)
                           cGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTg                                                                                                                       >  2:3515039/1‑150 (MQ=255)
                            gCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGt                                                                                                                      <  1:4502178/150‑1 (MQ=255)
                            gCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGt                                                                                                                      >  2:1643745/1‑150 (MQ=255)
                            gCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGt                                                                                                                      <  2:106996/150‑1 (MQ=255)
                            gCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGt                                                                                                                      <  1:1704780/150‑1 (MQ=255)
                              cGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCg                                                                                                                    >  1:276638/1‑150 (MQ=255)
                              cGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCg                                                                                                                    <  2:4587828/150‑1 (MQ=255)
                                 gtTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAAc                                                                                                                 >  1:3130542/1‑150 (MQ=255)
                                     acgacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATgg                                                                                                             <  1:234976/150‑1 (MQ=255)
                                     acgacATCAGGAGGCCAGTATGACAGTATTACATAGTGTGGATTTTTTACCGTCAGGTAACGCGTCCGTGGCGATAGAACCCCGTCTCCCGCAGCCGGATTTTCCTGAAAATCATCATGATTTTAATGAAATAGTGATTGTCGAACATgg                                                                                                             >  2:926570/1‑150 (MQ=255)
                                       gacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCa                                                                                                           >  2:3058658/1‑150 (MQ=255)
                                       gacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCa                                                                                                           <  1:2022669/150‑1 (MQ=255)
                                         cATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACg                                                                                                         <  2:1769067/150‑1 (MQ=255)
                                         cATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACg                                                                                                         >  2:1790626/1‑150 (MQ=255)
                                          aTCAGGAGGCCAGTATGACCGTATTAAATAGTGTGGATTTTTTTCCGTATGGTAACGCGTCCGTGGCGATATAACCCCGGCTAACGCAGGCGGATTTTCATGAAAAACATCATGATTTTCATGAAATTGTGAGTGTCGAACATGTCACgg                                                                                                        >  2:2342067/1‑150 (MQ=255)
                                            cAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGt                                                                                                      <  1:2756705/150‑1 (MQ=255)
                                            cAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGt                                                                                                      >  1:3391541/1‑150 (MQ=255)
                                                aggCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTc                                                                                                  >  2:1863658/1‑150 (MQ=255)
                                                 ggCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCa                                                                                                 <  2:2593754/150‑1 (MQ=255)
                                                  gCCAGTATGACCGTATTACCTAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAt                                                                                                <  2:3453624/150‑1 (MQ=255)
                                                        aTGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGttt                                                                                          >  1:1466395/1‑150 (MQ=255)
                                                         tGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCTAACATGGCACGGGTATTCATGTGTTTa                                                                                         <  2:3837976/150‑1 (MQ=255)
                                                         tGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTa                                                                                         >  2:2923143/1‑150 (MQ=255)
                                                         tGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTa                                                                                         >  2:2670521/1‑150 (MQ=255)
                                                           aCCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAAt                                                                                       >  1:3369197/1‑150 (MQ=255)
                                                              gTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATggg                                                                                    >  2:3183646/1‑150 (MQ=255)
                                                                aTTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCa                                                                                  >  2:3584573/1‑150 (MQ=255)
                                                                   aCATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGcc                                                                               <  1:4620472/150‑1 (MQ=255)
                                                                   aCATAGCGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGcc                                                                               <  1:2187678/150‑1 (MQ=255)
                                                                          gtgGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATAcc                                                                        >  2:1083406/1‑150 (MQ=255)
                                                                           tgGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCa                                                                       <  1:1392181/150‑1 (MQ=255)
                                                                            ggATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCAt                                                                      <  1:4206305/150‑1 (MQ=255)
                                                                                     tCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTgg                                                             <  2:1015606/150‑1 (MQ=255)
                                                                                     tCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTgg                                                             >  1:558502/1‑150 (MQ=255)
                                                                                       cGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCa                                                           <  1:2970589/150‑1 (MQ=255)
                                                                                       cGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCa                                                           >  1:2625799/1‑150 (MQ=255)
                                                                                        gTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCAc                                                          >  1:2005437/1‑150 (MQ=255)
                                                                                         tCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACg                                                         <  1:4124957/150‑1 (MQ=255)
                                                                                          cTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAGGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACgg                                                        >  2:2611503/1‑150 (MQ=255)
                                                                                           tGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGt                                                       <  1:3617790/150‑1 (MQ=255)
                                                                                             gTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCt                                                     <  2:4452462/150‑1 (MQ=255)
                                                                                               aaCGCGTCCGTGGAGATAGAACCCCGTCTCCCGCAGGCGGATTTTCCTGAACATAATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGt                                                   <  1:1426876/150‑1 (MQ=255)
                                                                                                     tCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTa                                             <  2:4486819/150‑1 (MQ=255)
                                                                                                      ccGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTAc                                            <  1:1364895/150‑1 (MQ=255)
                                                                                                       cGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTAcg                                           <  1:1200443/150‑1 (MQ=255)
                                                                                                            cGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATc                                      <  1:103167/150‑1 (MQ=255)
                                                                                                             gATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCa                                     <  1:2941019/150‑1 (MQ=255)
                                                                                                              aTAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCAt                                    <  2:1085869/150‑1 (MQ=255)
                                                                                                              aTAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCAt                                    >  2:3306421/1‑150 (MQ=255)
                                                                                                              aTAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCAt                                    <  1:50313/150‑1 (MQ=255)
                                                                                                                aGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGa                                  >  1:2454121/1‑150 (MQ=255)
                                                                                                                 gAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGAt                                 <  1:1006421/150‑1 (MQ=255)
                                                                                                                 gAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTAATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGAt                                 <  1:107400/150‑1 (MQ=255)
                                                                                                                  aaCCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATc                                >  1:4039733/1‑150 (MQ=255)
                                                                                                                  aaCCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATc                                <  2:1904890/150‑1 (MQ=255)
                                                                                                                        ggCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATc                          >  1:731692/1‑150 (MQ=255)
                                                                                                                        ggCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATc                          >  2:3717491/1‑150 (MQ=255)
                                                                                                                         gCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCt                         <  2:449540/150‑1 (MQ=255)
                                                                                                                          cTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTg                        <  1:3891557/150‑1 (MQ=255)
                                                                                                                          cTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTg                        <  1:2047053/150‑1 (MQ=255)
                                                                                                                           tCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGt                       <  2:2500457/150‑1 (MQ=255)
                                                                                                                                cAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGaa                  <  1:1123031/150‑1 (MQ=255)
                                                                                                                                 aGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAAc                 <  1:34280/150‑1 (MQ=255)
                                                                                                                                   gCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACAt               <  2:61064/150‑1 (MQ=255)
                                                                                                                                    cGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATa              <  1:1968231/150‑1 (MQ=255)
                                                                                                                                     ggATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATAc             >  2:2943978/1‑150 (MQ=255)
                                                                                                                                     ggATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATAc             >  1:3341784/1‑150 (MQ=255)
                                                                                                                                       aTTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCg           >  2:1985929/1‑150 (MQ=255)
                                                                                                                                        ttttCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGa          <  2:2208676/150‑1 (MQ=255)
                                                                                                                                        ttttCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGa          >  1:3540926/1‑150 (MQ=255)
                                                                                                                                           tCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATaa       >  2:3010663/1‑150 (MQ=255)
                                                                                                                                            ccTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAAt      <  1:948200/150‑1 (MQ=255)
                                                                                                                                            ccTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAAt      >  1:633981/1‑150 (MQ=255)
                                                                                                                                             cTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATc     <  1:4014065/150‑1 (MQ=255)
                                                                                                                                              tGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCt    <  1:4187201/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGT  >  CP009273/4087608‑4087901

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: