Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,570,803 G→A 7.9% I159I (ATC→ATT hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,570,8030GA 7.9% 992.2 / 53.8 279I159I (ATC→ATThsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (118/139);  new base A (8/14);  total (126/153)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 289 total aligned reads displayed.

AATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATAC  >  CP009273/4570654‑4570952
                                                                                                                                                     |                                                                                                                                                     
aaTTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTggg                                                                                                                                                       <  2:453858/150‑1 (MQ=255)
  ttttCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGAt                                                                                                                                                     >  1:3147469/1‑150 (MQ=255)
   tttCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGAtt                                                                                                                                                    <  1:630827/150‑1 (MQ=255)
        aTTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGtat                                                                                                                                               <  2:2142441/150‑1 (MQ=255)
          tAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGAGTTGTGGGATTTGCTCAAAACGTGCGTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTGTTGTTCGGCCCGTGGTGGGCTTGGTatat                                                                                                                                             >  2:1636500/1‑150 (MQ=255)
           aaCTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATAtt                                                                                                                                            >  2:240206/1‑150 (MQ=255)
           aaCTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGAATTGGTATAtt                                                                                                                                            >  1:4076133/1‑150 (MQ=255)
             cTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTa                                                                                                                                          <  1:4938154/150‑1 (MQ=255)
              tGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTat                                                                                                                                         >  1:1966322/1‑150 (MQ=255)
                cgcCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTatca                                                                                                                                       <  2:3273744/150‑1 (MQ=255)
                    ccccAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTatcaaatc                                                                                                                                   >  2:1000331/1‑150 (MQ=255)
                     cccAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTatcaaatca                                                                                                                                  <  1:3388653/150‑1 (MQ=255)
                      caaaTACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTatat                                                                                                                                             >  2:3931145/3‑138 (MQ=255)
                       cAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTAtcaaatcaaa                                                                                                                                <  2:780658/150‑1 (MQ=255)
                          tACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGACAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCt                                                                                                                             >  2:2707369/1‑150 (MQ=255)
                           aCCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCtt                                                                                                                            <  1:171975/150‑1 (MQ=255)
                            ccGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTg                                                                                                                           <  2:1269423/150‑1 (MQ=255)
                               cTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCg                                                                                                                        <  1:3540154/150‑1 (MQ=255)
                                   aCGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCtt                                                                                                                    >  2:837444/1‑150 (MQ=255)
                                   aCGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCtt                                                                                                                    >  2:2728057/1‑150 (MQ=255)
                                       aaaCGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAata                                                                                                                <  2:1024947/150‑1 (MQ=255)
                                          cGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAAtatta                                                                                                             <  1:2390384/150‑1 (MQ=255)
                                           gTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAAtattat                                                                                                            >  1:4802442/1‑150 (MQ=255)
                                            tttCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATattatt                                                                                                           >  2:1690056/1‑150 (MQ=255)
                                               cAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAat                                                                                                        <  2:1343841/150‑1 (MQ=255)
                                                      tgtgGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCa                                                                                                 >  1:4829771/1‑150 (MQ=255)
                                                        tgGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCAcc                                                                                               >  1:2266605/1‑150 (MQ=255)
                                                         gggATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACca                                                                                              <  2:2429725/150‑1 (MQ=255)
                                                         gggATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACca                                                                                              >  1:1107863/1‑150 (MQ=255)
                                                         gggATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACca                                                                                              <  1:1869414/150‑1 (MQ=255)
                                                          ggATTTGCTCAGAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACcag                                                                                             >  2:3540014/1‑150 (MQ=255)
                                                             tttGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACcagcag                                                                                          >  1:2770381/1‑150 (MQ=255)
                                                                gCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGaaa                                                                                       <  1:838649/150‑1 (MQ=255)
                                                                gCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGaaa                                                                                       >  1:3675261/1‑150 (MQ=255)
                                                                  tCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGt                                                                                     <  1:1003421/150‑1 (MQ=255)
                                                                   cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGTTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGtg                                                                                    >  1:767085/1‑150 (MQ=255)
                                                                   cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGtg                                                                                    <  1:2600687/150‑1 (MQ=255)
                                                                      aaCGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGtgatg                                                                                 <  1:3327201/150‑1 (MQ=255)
                                                                      aaCGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGtgatg                                                                                 <  1:5101054/150‑1 (MQ=255)
                                                                      aaCGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGtgatg                                                                                 >  2:3798539/1‑150 (MQ=255)
                                                                            cTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAAttt                                                                           >  1:2478314/1‑150 (MQ=255)
                                                                                ggTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGtt                                                                       >  1:4105127/1‑150 (MQ=255)
                                                                                 gTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGttt                                                                      <  2:4054930/150‑1 (MQ=255)
                                                                                   gCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCg                                                                    <  1:1477491/150‑1 (MQ=255)
                                                                                   gCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCg                                                                    <  1:3841660/150‑1 (MQ=255)
                                                                                     tGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGAt                                                                  <  1:910111/150‑1 (MQ=255)
                                                                                      gTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATa                                                                 >  1:4857918/1‑150 (MQ=255)
                                                                                      gTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATa                                                                 <  1:1443178/150‑1 (MQ=255)
                                                                                         tACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAag                                                              >  1:762964/1‑150 (MQ=255)
                                                                                          aCCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAaga                                                             <  1:397696/150‑1 (MQ=255)
                                                                                                  cagcagCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGAttt                                                     <  2:3612773/150‑1 (MQ=255)
                                                                                                    gcagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtg                                                   <  2:2878715/150‑1 (MQ=255)
                                                                                                    gcagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtg                                                   >  1:4675629/1‑150 (MQ=255)
                                                                                                    gcagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtg                                                   <  2:4839085/150‑1 (MQ=255)
                                                                                                     cagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgt                                                  <  1:1620974/150‑1 (MQ=255)
                                                                                                      agcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgtg                                                 <  2:4486793/150‑1 (MQ=255)
                                                                                                      agcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgtg                                                 >  1:4485936/1‑150 (MQ=255)
                                                                                                      agcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgtg                                                 >  2:223522/1‑150 (MQ=255)
                                                                                                      agcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgtg                                                 >  2:3458993/1‑150 (MQ=255)
                                                                                                       gcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGa                                                <  1:361532/150‑1 (MQ=255)
                                                                                                         gTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGaaa                                              <  2:4360858/150‑1 (MQ=255)
                                                                                                          tATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAAt                                             >  1:5119096/1‑150 (MQ=255)
                                                                                                             cGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAg                                          <  2:2648195/150‑1 (MQ=255)
                                                                                                             cGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAg                                          >  2:5054639/1‑150 (MQ=255)
                                                                                                               aGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCa                                        >  1:2122693/1‑150 (MQ=255)
                                                                                                                    ttCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATaaa                                   >  1:295206/1‑150 (MQ=255)
                                                                                                                    ttCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATaaa                                   <  2:4317728/150‑1 (MQ=255)
                                                                                                                       aGCGATGATTTTTTGTTCGGCATGTGGTGGGATTGGTATATTTATCAAATCAAAGCGTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCTATAAAGAGAAAATTTTGTGAAATGAGCAATAAAAcc                                <  1:657158/150‑1 (MQ=255)
                                                                                                                       aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAAcc                                <  2:4880395/150‑1 (MQ=255)
                                                                                                                         cGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAg                              >  1:4244359/1‑150 (MQ=255)
                                                                                                                          gatgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGa                             >  2:1367540/1‑150 (MQ=255)
                                                                                                                          gatgatTTTTTGTTCGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGa                             >  2:193512/1‑150 (MQ=255)
                                                                                                                          gatgatTTTTTGTTCGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGa                             >  2:193366/1‑150 (MQ=255)
                                                                                                                           atgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGaa                            <  2:2822903/150‑1 (MQ=255)
                                                                                                                            tgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGaaa                           <  1:2643341/150‑1 (MQ=255)
                                                                                                                              atTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAat                         >  2:3717871/1‑150 (MQ=255)
                                                                                                                               ttttttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAata                        <  1:973458/150‑1 (MQ=255)
                                                                                                                                tttttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAatat                       >  1:4144475/1‑150 (MQ=255)
                                                                                                                                  tttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATTAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATaa                     <  1:3994770/150‑1 (MQ=255)
                                                                                                                                  tttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATaa                     <  1:1355172/150‑1 (MQ=255)
                                                                                                                                      ttCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttt                 >  2:473644/1‑150 (MQ=255)
                                                                                                                                      ttCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttt                 <  2:3579053/150‑1 (MQ=255)
                                                                                                                                        cGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt               >  1:3479708/1‑150 (MQ=255)
                                                                                                                                        cGGCAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt               <  2:1535652/150‑1 (MQ=255)
                                                                                                                                         ggCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTc              <  2:591606/150‑1 (MQ=255)
                                                                                                                                           cAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg            <  2:3620071/150‑1 (MQ=255)
                                                                                                                                           cAAGTGGTGGAATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg            <  2:3682129/150‑1 (MQ=255)
                                                                                                                                            aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAgg           <  1:3919394/150‑1 (MQ=255)
                                                                                                                                            aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAgg           >  1:5050866/1‑150 (MQ=255)
                                                                                                                                            aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAgg           >  2:137171/1‑150 (MQ=255)
                                                                                                                                               tggtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACg        >  1:288418/1‑150 (MQ=255)
                                                                                                                                               tggtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACg        <  2:2529859/150‑1 (MQ=255)
                                                                                                                                                ggtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt       >  2:1833495/1‑150 (MQ=255)
                                                                                                                                                ggtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt       >  1:1277383/1‑150 (MQ=255)
                                                                                                                                                ggtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt       <  1:1824653/150‑1 (MQ=255)
                                                                                                                                                  tggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTaa     >  2:31173/1‑150 (MQ=255)
                                                                                                                                                    ggATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATa   <  1:4404693/150‑1 (MQ=255)
                                                                                                                                                    ggATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATa   <  1:677792/150‑1 (MQ=255)
                                                                                                                                                     gATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATac  >  1:4408412/1‑150 (MQ=255)
                                                                                                                                                     gATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATac  >  1:3078954/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATAC  >  CP009273/4570654‑4570952

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: