Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,017 C→T 17.6% intergenic (+180/+43) USA300HOU_pUSA300HOUMR0014 → / ← USA300HOU_pUSA300HOUMR0015 recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073112,0170CT17.6% 1260.8 / 232.6 515intergenic (+180/+43)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (236/188);  new base T (62/29);  total (298/217)
Fisher's exact test for biased strand distribution p-value = 3.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 518 total aligned reads displayed.

AAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11928‑12107
                                                                                         |                                                                                          
aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAACAACACAAGAAAAAAGCGTACAAATTCg                                                                                           >  1:334461/1‑91 (MQ=255)
 aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                          >  1:1051206/1‑91 (MQ=255)
 aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                          >  1:3168994/1‑91 (MQ=255)
   agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                        <  1:1983476/91‑1 (MQ=255)
   agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                        <  1:3584872/91‑1 (MQ=255)
    gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTGCGAATTTGTACg                                                                                       <  1:3197435/91‑1 (MQ=255)
    gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                       >  1:1668669/1‑91 (MQ=255)
     agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                      <  1:72326/91‑1 (MQ=255)
     agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                      <  1:2017133/91‑1 (MQ=255)
     agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                      <  1:2654441/91‑1 (MQ=255)
       aCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtt                                                                                    <  1:907199/91‑1 (MQ=255)
         gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTAAAAATTCGTACGCtttt                                                                                  >  1:1386978/1‑91 (MQ=255)
         gTTCGGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCgttt                                                                                  >  1:3273849/1‑91 (MQ=255)
           tCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGCtttttt                                                                                <  1:3890796/91‑1 (MQ=255)
              gACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                             <  1:1014646/91‑1 (MQ=255)
                cAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                           >  1:2562673/1‑91 (MQ=255)
                  attattaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                         >  1:2883633‑M2/42‑91 (MQ=255)
                   cGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataat                                                                        >  1:1834038/1‑91 (MQ=255)
                     tCTAAAATCTCATCATATTTTCCTGTTATGAGATTTTAATAAGGACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      <  1:4147339/91‑1 (MQ=255)
                     tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      <  1:2797134/91‑1 (MQ=255)
                      ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     >  1:1070257‑M2/38‑91 (MQ=255)
                      cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     >  1:1737596/1‑91 (MQ=255)
                      cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     >  1:3500724/1‑91 (MQ=255)
                      cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     >  1:244662/1‑91 (MQ=255)
                       TAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                                    <  1:4111127‑M1/91‑1 (MQ=255)
                        AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   >  1:253325‑M1/1‑91 (MQ=255)
                        AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   >  1:287988‑M1/1‑91 (MQ=255)
                          AATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGtt                                                                 <  1:1274120‑M1/91‑1 (MQ=255)
                           atctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                                >  1:2045091‑M2/33‑91 (MQ=255)
                              tcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                             >  1:2297044‑M2/30‑91 (MQ=255)
                              tcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                             >  1:1109110‑M2/30‑91 (MQ=255)
                               cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                            >  1:769583‑M2/29‑91 (MQ=255)
                               cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATCCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                            >  1:3679851‑M2/29‑91 (MQ=255)
                                 TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  1:591986‑M1/1‑91 (MQ=255)
                                 TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  1:236990‑M1/1‑91 (MQ=255)
                                  AATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAAt                                                         >  1:932253‑M1/1‑91 (MQ=255)
                                   ATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATg                                                        >  1:1599221‑M1/1‑91 (MQ=255)
                                        ttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                                   >  1:2471820/1‑91 (MQ=255)
                                        aaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                                   >  1:1319956‑M2/20‑91 (MQ=255)
                                         tCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtt                                                  >  1:2091729/1‑91 (MQ=255)
                                          ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                                 >  1:1347597/1‑91 (MQ=255)
                                          ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                                 <  1:3797720/91‑1 (MQ=255)
                                           cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                                >  1:1883183/1‑91 (MQ=255)
                                           cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                                >  1:3136494/1‑91 (MQ=255)
                                           cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                                <  1:1134061/91‑1 (MQ=255)
                                             gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTaa                                              <  1:284406/91‑1 (MQ=255)
                                                  gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                         <  1:4051860‑M2/82‑1 (MQ=255)
                                                  gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                         >  1:289284‑M2/10‑91 (MQ=255)
                                                  gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                         <  1:3878750‑M2/82‑1 (MQ=255)
                                                    gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                       <  1:1926909‑M2/84‑1 (MQ=255)
                                                    gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                       <  1:2677566‑M2/84‑1 (MQ=255)
                                                     attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                      <  1:146736‑M2/85‑1 (MQ=255)
                                                      ttttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                     <  1:277779‑M2/86‑1 (MQ=255)
                                                       tttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcat                                    <  1:2934917‑M2/87‑1 (MQ=255)
                                                       tttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcat                                    >  1:85411‑M2/5‑91 (MQ=255)
                                                          aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                                 <  1:3075876‑M2/90‑1 (MQ=255)
                                                          aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                                 >  1:3442142‑M2/2‑91 (MQ=255)
                                                            taataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCAGTTTCACTAATGAACGTTTTAAGTATCTTCATCAttc                               >  1:1793186/1‑90 (MQ=255)
                                                             aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  1:1405492/91‑1 (MQ=255)
                                                             aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  1:487591/91‑1 (MQ=255)
                                                                 acacAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGa                          <  1:3372379/91‑1 (MQ=255)
                                                                   acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                        >  1:3618767/1‑91 (MQ=255)
                                                                    caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                       >  1:566072/1‑91 (MQ=255)
                                                                    caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                       <  1:3739474/91‑1 (MQ=255)
                                                                     aaGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACt                      >  1:2202702/1‑91 (MQ=255)
                                                                     aaGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACt                      >  1:2325049/1‑91 (MQ=255)
                                                                      aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAcca                     >  1:1630566/1‑89 (MQ=255)
                                                                       gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc                    >  1:791862/1‑91 (MQ=255)
                                                                        aaaaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACa                   >  1:4012852/1‑91 (MQ=255)
                                                                        aaaaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACa                   >  1:438698/1‑91 (MQ=255)
                                                                           aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                >  1:1701849/1‑91 (MQ=255)
                                                                           aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                >  1:623813/1‑91 (MQ=255)
                                                                            aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  1:1221473/1‑91 (MQ=255)
                                                                            aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  1:3017939/1‑91 (MQ=255)
                                                                             aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat              <  1:2393591/91‑1 (MQ=255)
                                                                              gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  1:3960271/1‑91 (MQ=255)
                                                                              gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  1:4082436/1‑91 (MQ=255)
                                                                               cGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            >  1:3457034/1‑91 (MQ=255)
                                                                               cGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            >  1:3245202/1‑91 (MQ=255)
                                                                               cGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            >  1:2588993/1‑91 (MQ=255)
                                                                               cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            <  1:3708669/91‑1 (MQ=255)
                                                                               cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            <  1:2719143/91‑1 (MQ=255)
                                                                                gTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTCCAGTaataat           >  1:339106/1‑91 (MQ=255)
                                                                                gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           <  1:1149374/91‑1 (MQ=255)
                                                                                gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           >  1:3320068/1‑91 (MQ=255)
                                                                                gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           <  1:2756185/91‑1 (MQ=255)
                                                                                 tACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAataata          >  1:969801/1‑91 (MQ=255)
                                                                                  acgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  1:2507202/4‑91 (MQ=255)
                                                                                  acgaaTTTGTACGCTTTTCTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  1:2976558/4‑91 (MQ=255)
                                                                                   cgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:3419703/3‑91 (MQ=255)
                                                                                   cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:1176022/1‑91 (MQ=255)
                                                                                   cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:2955121/1‑91 (MQ=255)
                                                                                    aaaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc       >  1:1483274/1‑91 (MQ=255)
                                                                                     aaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAATTATCTTCATCATTCACTGAAACTACAGTAATAACACTCg      >  1:3852838/1‑91 (MQ=255)
                                                                                      aTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGt     >  1:3648300/1‑91 (MQ=255)
                                                                                      aTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGt     >  1:2231808/1‑91 (MQ=255)
                                                                                      aTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGt     <  1:2143191/91‑1 (MQ=255)
                                                                                       ttCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTc    >  1:1504142/1‑91 (MQ=255)
                                                                                       ttCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTc    >  1:8243/1‑91 (MQ=255)
                                                                                        tCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTca   <  1:1436063/91‑1 (MQ=255)
                                                                                         |                                                                                          
AAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11928‑12107

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