Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,946,267 G→T 13.0% I1032I (ATC→ATA recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,946,2670GT13.0% 975.9 / 91.3 293I1032I (ATC→ATArecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (110/145);  new base T (15/23);  total (125/168)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 297 total aligned reads displayed.

TTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTC  >  CP009273/2946118‑2946416
                                                                                                                                                     |                                                                                                                                                     
ttATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCg                                                                                                                                                       >  1:19614/1‑150 (MQ=255)
  ataGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTAt                                                                                                                                                     <  1:1547804/150‑1 (MQ=255)
        gAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGTgg                                                                                                                                               >  2:4380746/1‑150 (MQ=255)
        gAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCgg                                                                                                                                               <  1:3582247/150‑1 (MQ=255)
        gAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCgg                                                                                                                                               <  2:3969590/150‑1 (MQ=255)
        gAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCgg                                                                                                                                               >  1:1397456/1‑150 (MQ=255)
         aGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGt                                                                                                                                              <  2:625687/150‑1 (MQ=255)
          gCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGtt                                                                                                                                             <  1:1461002/150‑1 (MQ=255)
           cAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTc                                                                                                                                            <  2:3656414/150‑1 (MQ=255)
           cAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTc                                                                                                                                            <  1:798251/150‑1 (MQ=255)
           cAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTc                                                                                                                                            <  1:2495510/150‑1 (MQ=255)
            aGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCa                                                                                                                                           <  1:3730150/150‑1 (MQ=255)
            aGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCa                                                                                                                                           <  1:2931743/150‑1 (MQ=255)
                aataaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTaa                                                                                                                                       >  2:1551225/1‑150 (MQ=255)
                 ataaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAAt                                                                                                                                      <  2:1045583/150‑1 (MQ=255)
                    aCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCgg                                                                                                                                   >  1:2642810/1‑150 (MQ=255)
                     cGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGc                                                                                                                                  >  1:10693/1‑150 (MQ=255)
                      gCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCa                                                                                                                                 <  1:1180113/150‑1 (MQ=255)
                         cTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGAt                                                                                                                              >  2:2219128/1‑150 (MQ=255)
                         cTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGAt                                                                                                                              >  1:455364/1‑150 (MQ=255)
                          ttCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATa                                                                                                                             <  2:3484826/150‑1 (MQ=255)
                                gCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTc                                                                                                                       >  2:3349805/1‑150 (MQ=255)
                                   gAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCAt                                                                                                                    <  2:3094040/150‑1 (MQ=255)
                                     acacCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCAtct                                                                                                                  <  1:2147464/150‑1 (MQ=255)
                                       acCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTcc                                                                                                                >  2:2603490/1‑150 (MQ=255)
                                       acCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTcc                                                                                                                >  1:4138401/1‑150 (MQ=255)
                                        ccAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCa                                                                                                               >  1:2004025/1‑150 (MQ=255)
                                         cAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCAc                                                                                                              >  2:2346540/1‑150 (MQ=255)
                                            gTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTg                                                                                                           >  1:2884156/1‑150 (MQ=255)
                                             tCGATACAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGt                                                                                                          <  2:3159461/150‑1 (MQ=255)
                                             tCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGt                                                                                                          <  2:1958511/150‑1 (MQ=255)
                                             tCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGt                                                                                                          <  2:892025/150‑1 (MQ=255)
                                              cGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGtt                                                                                                         <  2:1275342/150‑1 (MQ=255)
                                                aTAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTa                                                                                                       <  1:232840/150‑1 (MQ=255)
                                                  aaaGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTAtt                                                                                                     >  2:2375638/1‑150 (MQ=255)
                                                    aGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgc                                                                                                   <  2:3580573/150‑1 (MQ=255)
                                                     gCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcg                                                                                                  >  1:3280007/1‑150 (MQ=255)
                                                              cATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAg                                                                                         <  2:776331/150‑1 (MQ=255)
                                                               aTGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGt                                                                                        <  2:3323692/150‑1 (MQ=255)
                                                                 gCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTg                                                                                      <  2:4036741/150‑1 (MQ=255)
                                                                  ccACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGa                                                                                     >  1:1828235/1‑150 (MQ=255)
                                                                  ccACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGa                                                                                     <  1:2676612/150‑1 (MQ=255)
                                                                   cACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGAc                                                                                    >  2:1009280/1‑150 (MQ=255)
                                                                     cGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCGTCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTc                                                                                  >  2:4521979/1‑150 (MQ=255)
                                                                      gTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCa                                                                                 <  1:3395556/150‑1 (MQ=255)
                                                                                gAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGcc                                                                       >  2:2192217/1‑150 (MQ=255)
                                                                                gAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGcc                                                                       <  2:180264/150‑1 (MQ=255)
                                                                                  aCTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCgg                                                                     >  1:862746/1‑150 (MQ=255)
                                                                                     ccAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGttt                                                                  <  2:3301944/150‑1 (MQ=255)
                                                                                     ccAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGttt                                                                  <  2:29709/150‑1 (MQ=255)
                                                                                     ccAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGttt                                                                  <  2:3907552/150‑1 (MQ=255)
                                                                                      cAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTc                                                                 <  1:4262727/150‑1 (MQ=255)
                                                                                      cAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTc                                                                 >  2:4257660/1‑150 (MQ=255)
                                                                                       aGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCa                                                                >  1:2417832/1‑150 (MQ=255)
                                                                                       aGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCa                                                                <  2:1376512/150‑1 (MQ=255)
                                                                                        gcggcgGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCAt                                                               <  1:4300222/150‑1 (MQ=255)
                                                                                          ggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTg                                                             <  1:709968/150‑1 (MQ=255)
                                                                                          ggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTg                                                             <  2:4087029/150‑1 (MQ=255)
                                                                                            cggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTgag                                                           >  2:2054047/1‑150 (MQ=255)
                                                                                                cAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGt                                                       <  1:2805420/150‑1 (MQ=255)
                                                                                                 aGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTg                                                      >  2:2042555/1‑150 (MQ=255)
                                                                                                 aGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTg                                                      >  1:1535475/1‑150 (MQ=255)
                                                                                                   ccTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGcc                                                    >  1:3898670/1‑150 (MQ=255)
                                                                                                   ccTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGcc                                                    <  2:3770045/150‑1 (MQ=255)
                                                                                                    cTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCt                                                   <  1:2440560/150‑1 (MQ=255)
                                                                                                    cTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCt                                                   <  1:4464715/150‑1 (MQ=255)
                                                                                                     tGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTg                                                  <  2:551069/150‑1 (MQ=255)
                                                                                                     tGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTg                                                  <  2:550930/150‑1 (MQ=255)
                                                                                                     tGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTg                                                  <  2:656943/150‑1 (MQ=255)
                                                                                                        ggATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTggag                                               <  1:2051047/150‑1 (MQ=255)
                                                                                                         gATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTggagg                                              <  1:2720239/150‑1 (MQ=255)
                                                                                                         gATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTggagg                                              >  2:3440931/1‑150 (MQ=255)
                                                                                                          aTAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTggagga                                             <  2:22494/150‑1 (MQ=255)
                                                                                                            aGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACa                                           <  1:3586808/150‑1 (MQ=255)
                                                                                                            aGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACa                                           >  1:3993222/1‑150 (MQ=255)
                                                                                                              cGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGc                                         >  1:716971/1‑150 (MQ=255)
                                                                                                              cGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGc                                         >  2:1106546/1‑150 (MQ=255)
                                                                                                                ggTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCg                                       >  2:1107884/1‑150 (MQ=255)
                                                                                                                    aaaCTGGCGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGAt                                   <  1:2752526/150‑1 (MQ=255)
                                                                                                                    aaaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGAt                                   >  1:2119878/1‑150 (MQ=255)
                                                                                                                     aaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATc                                  >  2:1486851/1‑150 (MQ=255)
                                                                                                                       cTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCa                                >  2:275444/1‑150 (MQ=255)
                                                                                                                       cTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCa                                <  1:1343534/150‑1 (MQ=255)
                                                                                                                        tGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCAc                               >  1:3844012/1‑150 (MQ=255)
                                                                                                                        tGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCAc                               <  2:1725939/150‑1 (MQ=255)
                                                                                                                         ggCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACt                              <  2:712145/150‑1 (MQ=255)
                                                                                                                         ggCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACt                              >  2:506932/1‑150 (MQ=255)
                                                                                                                         ggCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACt                              <  2:455364/150‑1 (MQ=255)
                                                                                                                           cGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCg                            <  2:3280007/150‑1 (MQ=255)
                                                                                                                           cGGATTAACGTATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCg                            <  1:1216548/150‑1 (MQ=255)
                                                                                                                           cGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCg                            >  1:2552736/1‑150 (MQ=255)
                                                                                                                           cGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCg                            >  2:2491258/1‑150 (MQ=255)
                                                                                                                           cGGATTAACGGATCAAGCTGACTGGCTATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCg                            <  1:3502689/150‑1 (MQ=255)
                                                                                                                               ttAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCa                        <  2:2429786/150‑1 (MQ=255)
                                                                                                                               ttAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCa                        <  1:1039260/150‑1 (MQ=255)
                                                                                                                                 aaCGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAAt                      <  1:2192217/150‑1 (MQ=255)
                                                                                                                                        cAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGtt               >  2:283126/1‑150 (MQ=255)
                                                                                                                                          aGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTcc             >  1:4111552/1‑150 (MQ=255)
                                                                                                                                             tGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCAc          >  1:1012838/1‑150 (MQ=255)
                                                                                                                                                    cGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGAtt   >  1:3168360/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTC  >  CP009273/2946118‑2946416

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: