Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,012 A→G 5.2% intergenic (+175/+48) USA300HOU_pUSA300HOUMR0014 → / ← USA300HOU_pUSA300HOUMR0015 recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073112,0120AG 5.2% 4703.4 / 141.1 1350intergenic (+175/+48)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (683/595);  new base G (20/51);  total (704/646)
Fisher's exact test for biased strand distribution p-value = 4.19e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 1382 total aligned reads displayed.

CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102
                                                                                          |                                                                                          
cGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACa                                                                                            >  2:4205963/1‑91 (MQ=255)
 gATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACaa                                                                                           >  2:1205862/1‑91 (MQ=255)
 gATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACaa                                                                                           >  1:4928031/1‑91 (MQ=255)
   tGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAt                                                                                         <  1:1084073/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  1:1490708/1‑91 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  2:2588302/1‑91 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttg                                                                                      <  2:1411079/91‑3 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttg                                                                                      <  1:1970463/91‑3 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                      <  1:2141914/91‑1 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                      <  1:2623938/91‑1 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:582410/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:3580712/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                     >  1:3246212/1‑91 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                   <  1:3169588/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  2:5009198/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 >  1:5251762/1‑91 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 <  2:2171078/91‑1 (MQ=255)
               gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttt                                                                             >  2:4806191/1‑91 (MQ=255)
               gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttt                                                                             >  1:2218686/1‑91 (MQ=255)
                ttCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCttttt                                                                            <  1:125697/91‑1 (MQ=255)
                 tCTGACAGCGTATAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttttt                                                                           <  2:5641045/91‑1 (MQ=255)
                  cTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTta                                                                          <  2:3007239/91‑1 (MQ=255)
                    gACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                        <  1:1154602/91‑1 (MQ=255)
                      cAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                      >  1:2508383/1‑91 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    >  2:4025899/1‑91 (MQ=255)
                         cGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataat                                                                   <  2:3900007/91‑1 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 >  1:555065/1‑91 (MQ=255)
                            cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                <  1:814921/91‑1 (MQ=255)
                             tAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  2:1829905/1‑91 (MQ=255)
                             tAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  1:3983323/1‑91 (MQ=255)
                             tAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  2:150128/1‑91 (MQ=255)
                               aaaTCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                             >  2:510751/1‑91 (MQ=255)
                                 aTCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                           <  2:3292591/91‑1 (MQ=255)
                                   ctcATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCa                                                         <  2:2660743/91‑1 (MQ=255)
                                    tcATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                        >  2:4378635/1‑91 (MQ=255)
                                       taataTTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     <  1:3470636/91‑1 (MQ=255)
                                       taataTTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     <  1:176415/91‑1 (MQ=255)
                                          tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                  <  2:1326017/91‑1 (MQ=255)
                                          aatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                  <  2:4116536/90‑1 (MQ=255)
                                               tCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtt                                             <  1:2333001/91‑1 (MQ=255)
                                                ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                            >  2:728540/1‑91 (MQ=255)
                                                ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                            >  1:1876725/1‑91 (MQ=255)
                                                ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                            >  1:923006/1‑91 (MQ=255)
                                                 cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                           <  2:3520006/91‑1 (MQ=255)
                                                  tGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACACTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTa                                          >  2:4703025/1‑91 (MQ=255)
                                                   gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTaa                                         >  2:3947011/1‑91 (MQ=255)
                                                    tgaTGAGATTGTAATAATAACACAAGAAATAAGCGTACATATTCTTTCGTTTTTTTAATATAATACCTGTTTCACTAATGATCGTTTTAAg                                        <  2:4477452/89‑1 (MQ=255)
                                                     tATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGt                                       <  1:4328013/91‑1 (MQ=255)
                                                      aTGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTa                                      <  1:5144202/91‑1 (MQ=255)
                                                       tGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTAt                                     >  2:3131860/1‑91 (MQ=255)
                                                       tGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTGAAGTAt                                     >  1:4725155/1‑91 (MQ=255)
                                                        gagaTTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                    <  2:162623/91‑1 (MQ=255)
                                                         agaTTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCt                                   <  2:4271648/91‑1 (MQ=255)
                                                          gaTTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                  <  1:1564528/91‑1 (MQ=255)
                                                          gaTTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                  <  1:379187/91‑1 (MQ=255)
                                                           aTTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 >  2:1125275/1‑91 (MQ=255)
                                                            ttttAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                <  2:910331/91‑1 (MQ=255)
                                                            ttttAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                <  2:288479/91‑1 (MQ=255)
                                                             tttAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcat                               >  2:536544/1‑91 (MQ=255)
                                                              ttAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatc                              <  1:4111048/91‑1 (MQ=255)
                                                                aataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  1:1404433/91‑1 (MQ=255)
                                                                aataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  2:5571889/91‑1 (MQ=255)
                                                                aataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  1:652766/1‑91 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           <  1:371065/91‑1 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           <  1:4183088/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:5463507/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:285960/91‑1 (MQ=255)
                                                                    ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATACAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                        >  2:5485268/1‑91 (MQ=255)
                                                                     taaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACt                       >  2:2946849/1‑91 (MQ=255)
                                                                     taaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACt                       >  2:2056012/1‑91 (MQ=255)
                                                                      aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                      >  1:5384769/1‑91 (MQ=255)
                                                                      aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                      >  1:4532918/1‑91 (MQ=255)
                                                                      aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                      <  2:2501974/91‑1 (MQ=255)
                                                                       acacAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGa                     <  2:4601477/91‑1 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATGTAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   <  1:4666660/91‑1 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  2:4932618/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  1:2701884/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  2:1742402/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  1:5022251/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   <  2:1895616/91‑1 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   <  2:619749/91‑1 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  >  2:3776537/1‑91 (MQ=255)
                                                                           aaGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACt                 >  1:4455974/1‑91 (MQ=255)
                                                                            aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                >  2:2330794/1‑91 (MQ=255)
                                                                                aaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGt            <  1:3879943/91‑1 (MQ=255)
                                                                                aaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTGTCACTAAGAAACGTTTTAAGTAACTTCATCATTAACTGAAACTACAGt            >  1:3368600/1‑91 (MQ=255)
                                                                                 aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta           >  1:3106394/1‑91 (MQ=255)
                                                                                 aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta           >  1:4832441/1‑91 (MQ=255)
                                                                                 aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta           <  2:2759675/91‑1 (MQ=255)
                                                                                 aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta           <  2:2436950/91‑1 (MQ=255)
                                                                                 aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta           >  2:1998801/1‑91 (MQ=255)
                                                                                  aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa          <  1:2997899/91‑1 (MQ=255)
                                                                                  aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa          >  2:3697998/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  2:3233924/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGGACGCTTTTTTAATATAATACCTGTTTCACTACTGAACGTTTTAAGTATCTTCATCATTACCTGCAACTACAGtaata        >  2:831024/1‑91 (MQ=255)
                                                                                     cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa       <  1:2427257/91‑1 (MQ=255)
                                                                                      gTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      >  1:721078/1‑91 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      <  1:1005182/91‑1 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      <  2:5208360/91‑1 (MQ=255)
                                                                                          aaaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc  <  1:1659266/91‑1 (MQ=255)
                                                                                          |                                                                                          
CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: