Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,012 A→G 10.1% 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,0120AG10.1% 1859.2 / 144.2 656intergenic (+175/+48)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (307/281);  new base G (13/54);  total (321/335)
Fisher's exact test for biased strand distribution p-value = 2.42e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 671 total aligned reads displayed.

CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102
                                                                                          |                                                                                          
cGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACa                                                                                            >  2:395192/1‑91 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtt                                                                                        <  2:2157695/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  1:745523/1‑91 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  2:2389844/1‑91 (MQ=255)
     tAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAttt                                                                                       <  1:2801391/91‑2 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                      >  2:664480/1‑91 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  2:2784931/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                     >  2:3042722/1‑91 (MQ=255)
        aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTa                                                                                    <  1:1013313/91‑1 (MQ=255)
        aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTa                                                                                    <  2:3472570/91‑1 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                   <  2:1012132/91‑1 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                   <  2:2247675/91‑1 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                   <  2:3250299/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                  <  2:1916494/91‑1 (MQ=255)
           agaCGTTTTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGACTTTGTACGc                                                                                 <  1:1945977/91‑1 (MQ=38)
           agaCGTTCTGACAGCGTCTAAAATTTCATAATATTTTCCTGTTATGAGATGTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 <  2:3605738/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                 <  1:3156764/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 <  2:1206217/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTATCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 <  2:1511025/91‑1 (MQ=255)
            gaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCt                                                                                <  1:3051016/91‑1 (MQ=255)
               gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttt                                                                             >  2:2671589/1‑91 (MQ=255)
                 ttgtgttattattaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttttt                                                                           >  1:1745224‑M2/49‑91 (MQ=255)
                  cTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTta                                                                          <  2:2587558/91‑1 (MQ=255)
                  cTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTta                                                                          <  2:3830990/91‑1 (MQ=255)
                    gACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                        <  2:1833359/91‑1 (MQ=255)
                     aCAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGCttttttcttg                                                                       <  1:3706215/91‑5 (MQ=37)
                     aCAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaata                                                                       >  1:1436732/1‑91 (MQ=255)
                      cAGCGTCTAAAATCTCATAATGTTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                      >  2:3648990/1‑91 (MQ=255)
                      cAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                      >  2:2491334/1‑91 (MQ=255)
                      cAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                      >  1:3788/1‑91 (MQ=255)
                      cAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatat                                                                      <  1:3735564/91‑1 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    >  2:799429/1‑91 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    <  1:62029/91‑1 (MQ=255)
                         cGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataat                                                                   <  1:48327/91‑1 (MQ=255)
                          tattaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCATTTTtaatataata                                                                  >  2:2601769‑M2/40‑91 (MQ=255)
                          gTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataata                                                                  <  2:376406/91‑1 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 >  2:1616438/1‑91 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 <  1:1049946/91‑1 (MQ=255)
                           attaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 >  2:2881738‑M2/39‑91 (MQ=255)
                            ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  2:3013183‑M2/38‑91 (MQ=255)
                            ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  1:2972878‑M2/38‑91 (MQ=255)
                               aaatctcataacaggaaaatattatgagatttttATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                             >  2:454156‑M2/35‑91 (MQ=255)
                               aaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                             >  1:574491‑M2/35‑91 (MQ=255)
                               AAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                             >  1:3852855‑M1/1‑91 (MQ=255)
                                AATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGtt                                                            <  1:3577912‑M1/91‑1 (MQ=255)
                                  tctcCTAATACTTTCCTGTTATGAGCTTTTACTAATAACACAAGAAAAAAGCGTACCACTTGGTACGCTTTTCTAACCTAATACATGTTTc                                                          >  2:2703809/1‑91 (MQ=255)
                                  TCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTc                                                          <  2:650523‑M1/91‑1 (MQ=255)
                                   CTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCa                                                         >  1:2482618‑M1/1‑91 (MQ=255)
                                       TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     <  1:2332467‑M1/91‑1 (MQ=255)
                                        aacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAAt                                                    >  2:142075‑M2/26‑91 (MQ=255)
                                        AATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAAt                                                    >  2:1104553‑M1/1‑91 (MQ=255)
                                        AATATTTTCCTGTTATGAGATTTTAATAACAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAAt                                                    >  1:3111214‑M1/1‑91 (MQ=255)
                                         ATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATg                                                   >  1:2831724‑M1/1‑91 (MQ=255)
                                         ATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATg                                                   >  2:1970443‑M1/1‑91 (MQ=255)
                                          tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                  >  2:2991784/1‑91 (MQ=255)
                                          tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                  <  1:1640961/91‑1 (MQ=255)
                                           atTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGaa                                                 <  1:930660/91‑1 (MQ=255)
                                            ttgtCCTGTTATGCGATTTTAATAATAACACAAGAAAAAAGCTTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAAc                                                <  2:730586/88‑1 (MQ=255)
                                             tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                               <  2:3427195/91‑1 (MQ=255)
                                             tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                               >  2:2344964/1‑91 (MQ=255)
                                               tCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtt                                             >  1:1895546/1‑91 (MQ=255)
                                                 cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                           >  2:3545579/1‑91 (MQ=255)
                                                 atattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                           >  1:271267‑M2/17‑91 (MQ=255)
                                                  tGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTa                                          <  1:2147186/91‑1 (MQ=255)
                                                       tgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTAt                                     >  1:2579872‑M2/11‑91 (MQ=255)
                                                           attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 >  2:3153916‑M2/7‑91 (MQ=255)
                                                           attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 <  1:408866‑M2/85‑1 (MQ=255)
                                                           attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 <  1:3326539‑M2/85‑1 (MQ=255)
                                                            ttttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                >  2:146821‑M2/6‑91 (MQ=255)
                                                            ttttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                <  2:1350255‑M2/86‑1 (MQ=255)
                                                            ttgtaATGATAACGCAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATATGACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                <  2:908055‑M2/86‑1 (MQ=255)
                                                              ttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatc                              <  1:2639806‑M2/88‑1 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           >  1:2018534/1‑91 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:3648990/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:3387884/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:1587431/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  2:3698135/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         >  2:1737421/1‑91 (MQ=255)
                                                                        cacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                    >  1:967821/1‑91 (MQ=255)
                                                                        cacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                    >  2:568749/1‑91 (MQ=255)
                                                                        cacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                    >  2:2085896/1‑91 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  >  1:2413679/1‑91 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  <  1:2520492/91‑1 (MQ=255)
                                                                            aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                >  1:2877502/1‑91 (MQ=255)
                                                                             taaaaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc               >  2:3329810/2‑91 (MQ=255)
                                                                             gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc               >  1:383538/1‑91 (MQ=255)
                                                                                aaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGt            <  1:1813520/91‑1 (MQ=255)
                                                                                   aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         <  1:2084437/91‑1 (MQ=255)
                                                                                   aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         >  1:2350445/1‑91 (MQ=255)
                                                                                    gCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  1:2273044/1‑91 (MQ=255)
                                                                                    gCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  1:1209365/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  1:134681/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        <  1:3271574/91‑1 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  1:351930/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  2:3359485/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAAAATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  2:2219427/1‑91 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      <  2:201445/91‑1 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      >  1:1089489/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  2:505950/1‑91 (MQ=255)
                                                                                          aaaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc  <  1:1536837/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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: