Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,087,754 C→A 16.5% 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,7540CA16.5% 829.8 / 99.3 265H31N (CAT→AAT) rhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (113/108);  new base A (23/21);  total (136/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 287 total aligned reads displayed.

GGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTG  >  CP009273/4087605‑4087902
                                                                                                                                                     |                                                                                                                                                    
ggTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAc                                                                                                                                                      <  2:3026278/150‑1 (MQ=255)
 gTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAca                                                                                                                                                     >  1:3714380/1‑150 (MQ=255)
 gTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAca                                                                                                                                                     >  1:756153/1‑150 (MQ=255)
  tAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAaat                                                                                                                                                    >  1:1100610/1‑150 (MQ=255)
    aGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatca                                                                                                                                                  <  2:1596121/150‑1 (MQ=255)
    aGAACCTGACCTCGTGATTACTATTTCGACGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAaatca                                                                                                                                                  >  1:1374965/1‑150 (MQ=255)
     gAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcat                                                                                                                                                 >  1:727664/1‑150 (MQ=255)
      aaCCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatc                                                                                                                                                >  1:3053580/1‑150 (MQ=255)
      aaCCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatc                                                                                                                                                >  2:3308726/1‑150 (MQ=255)
        ccTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAcatcatcat                                                                                                                                              <  1:2082284/150‑1 (MQ=255)
         cTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATg                                                                                                                                             <  1:3567540/150‑1 (MQ=255)
          tGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGa                                                                                                                                            <  2:3448869/150‑1 (MQ=255)
           gACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGAt                                                                                                                                           <  1:3839349/150‑1 (MQ=255)
               tCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTc                                                                                                                                       >  2:3568343/1‑150 (MQ=255)
                  tGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATg                                                                                                                                    >  2:3067047/1‑150 (MQ=255)
                   gATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGa                                                                                                                                   <  1:2355244/150‑1 (MQ=255)
                   gATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGa                                                                                                                                   >  1:1439336/1‑150 (MQ=255)
                    aTTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGaa                                                                                                                                  >  2:945200/1‑150 (MQ=255)
                      tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAt                                                                                                                                >  1:2616508/1‑150 (MQ=255)
                      tACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACACCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAt                                                                                                                                <  1:3518161/150‑1 (MQ=255)
                       aCTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAtt                                                                                                                               >  1:1323906/1‑150 (MQ=255)
                       aCTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAAtt                                                                                                                               >  1:325842/1‑150 (MQ=255)
                       aCTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAAtt                                                                                                                               <  1:2852998/150‑1 (MQ=255)
                        cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg                                                                                                                              >  1:1213346/1‑150 (MQ=255)
                        cTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATtg                                                                                                                              <  2:3183472/150‑1 (MQ=255)
                         tATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATtgt                                                                                                                             >  2:2363098/1‑150 (MQ=255)
                             tCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGAtt                                                                                                                         <  1:4305576/150‑1 (MQ=255)
                                  gtgtTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGa                                                                                                                    <  1:2253817/150‑1 (MQ=255)
                                      tGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACAt                                                                                                                >  2:3386259/1‑150 (MQ=255)
                                      tGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACAt                                                                                                                >  1:3154381/1‑150 (MQ=255)
                                       gacgacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATg                                                                                                               <  1:4134147/150‑1 (MQ=255)
                                       gacgacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATg                                                                                                               <  1:3745301/150‑1 (MQ=255)
                                        acgacATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATgg                                                                                                              <  2:2737752/150‑1 (MQ=255)
                                            cATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACg                                                                                                          >  2:797926/1‑150 (MQ=255)
                                             aTCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACgg                                                                                                         <  1:1414591/150‑1 (MQ=255)
                                               cAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGt                                                                                                       >  1:186552/1‑150 (MQ=255)
                                                aggagtCCAGTATGAACGTATTACATAGTGTGGATTTATATACCTCTGGTAACGCGTCCGTGGCTATAGAAACACGGATCCCGCAGGCTTATTTTACTGAACATCATAAAGATTTTCATTAAAATGTGATTGTCGAACATGGCACGGGTa                                                                                                      >  2:1532774/1‑150 (MQ=255)
                                                   aggCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTc                                                                                                   >  1:2959851/1‑150 (MQ=255)
                                                     gCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAt                                                                                                 <  1:2535557/150‑1 (MQ=255)
                                                      ccAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtg                                                                                                >  1:215243/1‑150 (MQ=255)
                                                       cAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtgt                                                                                               <  1:4493357/150‑1 (MQ=255)
                                                       cAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtgt                                                                                               <  1:253997/150‑1 (MQ=255)
                                                       cAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtgt                                                                                               >  1:1462030/1‑150 (MQ=255)
                                                       cAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGAATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtgt                                                                                               >  1:2909806/1‑150 (MQ=255)
                                                        aGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCAtgtg                                                                                              >  1:3463294/1‑150 (MQ=255)
                                                         gTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATgtgt                                                                                             <  2:4560409/150‑1 (MQ=255)
                                                           aTGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGttt                                                                                           >  1:4618238/1‑150 (MQ=255)
                                                           aTGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGttt                                                                                           >  1:2663687/1‑150 (MQ=255)
                                                           aTGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGttt                                                                                           >  1:1497698/1‑150 (MQ=255)
                                                                 gTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATggg                                                                                     >  2:2616879/1‑150 (MQ=255)
                                                                   aTTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCa                                                                                   >  2:592458/1‑150 (MQ=255)
                                                                       cATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGccc                                                                               >  1:43501/1‑150 (MQ=255)
                                                                        aTAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCt                                                                              >  1:2204786/1‑150 (MQ=255)
                                                                        aTAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCt                                                                              >  2:1494381/1‑150 (MQ=255)
                                                                          aGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCtat                                                                            >  1:1686876/1‑150 (MQ=255)
                                                                             gtgGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATAcc                                                                         <  1:19722/150‑1 (MQ=255)
                                                                                      tttCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGt                                                                <  2:4350919/150‑1 (MQ=255)
                                                                                      tttCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGt                                                                <  2:3662054/150‑1 (MQ=255)
                                                                                        tCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTgg                                                              <  1:1293069/150‑1 (MQ=255)
                                                                                         ccGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGc                                                             <  1:1084600/150‑1 (MQ=255)
                                                                                            tCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTAGTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACg                                                          >  2:1111772/1‑150 (MQ=255)
                                                                                              tGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGt                                                        <  2:2278773/150‑1 (MQ=255)
                                                                                                  aaCGCGTCCGTGGCGATAGAACCCCGGCTCACTCAGGCGGATTTTCATGAACATCATCATGATTTTCATGAAATTGTTATTGTCTAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCACTATACCATCACCGGTGGCACGGTCTGt                                                    <  2:1711837/150‑1 (MQ=255)
                                                                                                       gTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGt                                               >  2:479385/1‑150 (MQ=255)
                                                                                                       gTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGt                                               >  2:33747/1‑150 (MQ=255)
                                                                                                          cGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTAcg                                            >  2:3825509/1‑150 (MQ=255)
                                                                                                            tGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTAcgcg                                          <  2:1401022/150‑1 (MQ=255)
                                                                                                            tGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTAcgcg                                          <  2:1272173/150‑1 (MQ=255)
                                                                                                              gCGATAGATCCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGAt                                        <  2:3817308/150‑1 (MQ=255)
                                                                                                                 aTAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCAt                                     <  1:3609824/150‑1 (MQ=255)
                                                                                                                    gAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGAt                                  <  2:1823519/150‑1 (MQ=255)
                                                                                                                     aaCCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATc                                 <  2:742630/150‑1 (MQ=255)
                                                                                                                     aaCCCCGGCTCCCGCAGGCGGATTTTCCTGAAAATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATc                                 >  2:3877365/1‑150 (MQ=255)
                                                                                                                       ccccGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCgg                               <  2:3147239/150‑1 (MQ=255)
                                                                                                                        cccGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGc                              >  1:1522350/1‑150 (MQ=255)
                                                                                                                         ccGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCCTGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCa                             >  1:3334825/1‑150 (MQ=255)
                                                                                                                         ccGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCa                             <  2:43501/150‑1 (MQ=255)
                                                                                                                          cGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCAt                            >  1:3121821/1‑150 (MQ=255)
                                                                                                                           ggCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATc                           >  1:2486897/1‑150 (MQ=255)
                                                                                                                           ggCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATc                           >  1:687457/1‑150 (MQ=255)
                                                                                                                             cTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTg                         <  2:3587428/150‑1 (MQ=255)
                                                                                                                             cTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTg                         <  2:3960678/150‑1 (MQ=255)
                                                                                                                              tCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGt                        <  1:1628484/150‑1 (MQ=255)
                                                                                                                               cccGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTa                       >  1:3026917/1‑150 (MQ=255)
                                                                                                                               cccGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTa                       >  1:494832/1‑150 (MQ=255)
                                                                                                                               cccGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTa                       >  1:2813247/1‑150 (MQ=255)
                                                                                                                                 cGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATTAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATg                     <  2:1970735/150‑1 (MQ=255)
                                                                                                                                 cGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATg                     <  1:106839/150‑1 (MQ=255)
                                                                                                                                    aGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAAc                  >  2:2534228/1‑150 (MQ=255)
                                                                                                                                    aGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGAGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAAc                  >  2:4105954/1‑150 (MQ=255)
                                                                                                                                      gCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACAt                <  1:86821/150‑1 (MQ=255)
                                                                                                                                      gCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACAt                <  2:2158288/150‑1 (MQ=255)
                                                                                                                                       cGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATa               <  2:3239234/150‑1 (MQ=255)
                                                                                                                                           ttttCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGa           >  1:522989/1‑150 (MQ=255)
                                                                                                                                              tCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATaa        <  2:2258044/150‑1 (MQ=255)
                                                                                                                                              tCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATaa        >  1:4432260/1‑150 (MQ=255)
                                                                                                                                                cTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATc      <  1:827170/150‑1 (MQ=255)
                                                                                                                                                 tGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCt     <  1:946333/150‑1 (MQ=255)
                                                                                                                                                  gAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCtg    >  1:4389308/1‑150 (MQ=255)
                                                                                                                                                   aaCATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCtgt   >  1:2806485/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTAAGAACCTGACCTCGTGATTACTATTTCGCCGTGTTGACGACATCAGGAGGCCAGTATGACCGTATTACATAGTGTGGATTTTTTTCCGTCTGGTAACGCGTCCGTGGCGATAGAACCCCGGCTCCCGCAGGCGGATTTTCCTGAACATCATCATGATTTTCATGAAATTGTGATTGTCGAACATGGCACGGGTATTCATGTGTTTAATGGGCAGCCCTATACCATCACCGGTGGCACGGTCTGTTTCGTACGCGATCATGATCGGCATCTGTATGAACATACCGATAATCTGTG  >  CP009273/4087605‑4087902

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