Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,012 A→G 14.0% 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,0120AG14.0% 1286.2 / 166.8 483intergenic (+175/+48)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (230/185);  new base G (25/43);  total (255/228)
Fisher's exact test for biased strand distribution p-value = 5.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 495 total aligned reads displayed.

CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102
                                                                                          |                                                                                          
cGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACa                                                                                            >  1:1702921/1‑91 (MQ=255)
  aTGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACaaa                                                                                          <  1:278374/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtt                                                                                        <  1:975430/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtt                                                                                        <  1:1142487/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  1:3234470/1‑91 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAtt                                                                                        >  1:2736778/1‑91 (MQ=255)
      aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttg                                                                                      <  1:2209378/91‑3 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:3433997/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:3253685/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:1686517/91‑4 (MQ=255)
        aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTa                                                                                    <  1:1886403/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  1:2133954/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  1:298760/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  1:2760805/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                 <  1:4123055/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                 <  1:2654441/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                 <  1:53311/91‑1 (MQ=255)
            gaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCt                                                                                <  1:1261363/91‑1 (MQ=255)
             aCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtt                                                                               >  1:772431/1‑91 (MQ=255)
              cGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCttt                                                                              >  1:1558363/1‑91 (MQ=255)
               gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttt                                                                             >  1:2023944/1‑91 (MQ=255)
                    gACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                        <  1:1014646/91‑1 (MQ=255)
                     aCAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaata                                                                       >  1:3417724/1‑91 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    <  1:1217897/91‑1 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTATtaatataa                                                                    >  1:3976614/1‑91 (MQ=255)
                        gCGTCTAAAATCTCACAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAACAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    >  1:3543775/1‑91 (MQ=255)
                         cGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataat                                                                   >  1:1834038/1‑91 (MQ=255)
                          gTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataata                                                                  <  1:3367110/91‑1 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 <  1:2797134/91‑1 (MQ=255)
                            ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  1:1070257‑M2/38‑91 (MQ=255)
                            cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  1:244662/1‑91 (MQ=255)
                            cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                <  1:834877/91‑1 (MQ=255)
                            cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  1:1737596/1‑91 (MQ=255)
                             taaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  1:2184348‑M2/37‑91 (MQ=255)
                             TAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  1:2883587‑M1/1‑91 (MQ=255)
                              aaaatctcataacaggaaaatattatgagattttaATCATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              >  1:3303717‑M2/36‑91 (MQ=255)
                               AAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTATTATAATACCTGt                                                             >  1:3199888‑M1/1‑91 (MQ=255)
                                 atctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                           >  1:2045091‑M2/33‑91 (MQ=255)
                                 ATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                           >  1:2567724‑M1/1‑91 (MQ=255)
                                  TCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTc                                                          <  1:103772‑M1/91‑1 (MQ=255)
                                     cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                       >  1:2276035‑M2/29‑91 (MQ=255)
                                     cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATCCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                       >  1:3679851‑M2/29‑91 (MQ=255)
                                     cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAACCGTACAAATTCGTACGCTTTTTTAATATAGTACCTGTTTCACt                                                       >  1:4149520‑M2/29‑91 (MQ=255)
                                      ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                      <  1:3606864‑M1/91‑1 (MQ=255)
                                         ATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATg                                                   <  1:711841‑M1/91‑1 (MQ=255)
                                          tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                  <  1:232757/91‑1 (MQ=255)
                                              ttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                              >  1:2471820/1‑91 (MQ=255)
                                               tCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtt                                             >  1:2091729/1‑91 (MQ=255)
                                                ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                            >  1:3268069/1‑91 (MQ=255)
                                                  tGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTa                                          <  1:2867532/91‑1 (MQ=255)
                                                   gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTaa                                         <  1:284406/91‑1 (MQ=255)
                                                   gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTaa                                         >  1:2522085/1‑91 (MQ=255)
                                                    ttatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAg                                        <  1:1429591‑M2/78‑1 (MQ=255)
                                                        gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                    >  1:289284‑M2/10‑91 (MQ=255)
                                                        gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                    <  1:3933359‑M2/82‑1 (MQ=255)
                                                         agattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCt                                   >  1:2378393‑M2/9‑91 (MQ=255)
                                                          gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                  <  1:868883‑M2/84‑1 (MQ=255)
                                                          gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                  <  1:2677566‑M2/84‑1 (MQ=255)
                                                           attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 <  1:146736‑M2/85‑1 (MQ=255)
                                                           attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                 >  1:3140235‑M2/7‑91 (MQ=255)
                                                             tttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcat                               <  1:2934917‑M2/87‑1 (MQ=255)
                                                               taATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatca                             >  1:795334‑M2/3‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  1:576321‑M2/2‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  1:3075876‑M2/90‑1 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  1:2302338‑M2/90‑1 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  1:2696117‑M2/2‑91 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           <  1:3456235/91‑1 (MQ=255)
                                                                    ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                        <  1:3335315/91‑1 (MQ=255)
                                                                    ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                        >  1:3174919/1‑91 (MQ=255)
                                                                    ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCGGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                        <  1:3883530/91‑1 (MQ=255)
                                                                     taaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACt                       <  1:1331776/91‑1 (MQ=255)
                                                                      aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                      >  1:3485505/1‑91 (MQ=255)
                                                                       acacAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGCATCTTCATCATTAACTGa                     >  1:1397568/1‑91 (MQ=255)
                                                                       acacAAGAAAAAAGCGTACAAATTCATACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGa                     <  1:3865722/91‑1 (MQ=255)
                                                                        ctcaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                    <  1:1654061/89‑1 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  >  1:911048/1‑91 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  <  1:3630054/91‑1 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                  <  1:3739474/91‑1 (MQ=255)
                                                                          caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTCCATCATTAACTGAAAc                  <  1:965328/91‑1 (MQ=255)
                                                                            ttaaaaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                >  1:1113681/3‑91 (MQ=38)
                                                                            aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                >  1:1938282/1‑91 (MQ=255)
                                                                             gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc               >  1:4020933/1‑91 (MQ=255)
                                                                             gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc               >  1:1796448/1‑91 (MQ=255)
                                                                                  aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa          >  1:3017939/1‑91 (MQ=255)
                                                                                   aGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         >  1:613091/1‑91 (MQ=255)
                                                                                   aGCGTACAAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         <  1:1306126/91‑1 (MQ=255)
                                                                                   aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         <  1:1363965/91‑1 (MQ=255)
                                                                                   aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         >  1:253853/1‑91 (MQ=255)
                                                                                   aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat         <  1:4083616/91‑1 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  1:1983987/1‑91 (MQ=255)
                                                                                     cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa       >  1:3830047/1‑91 (MQ=255)
                                                                                     cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa       >  1:1629113/1‑91 (MQ=255)
                                                                                     cGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa       <  1:3708669/91‑1 (MQ=255)
                                                                                      gTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTCCAGTaataat      >  1:339106/1‑91 (MQ=255)
                                                                                      gTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      >  1:449845/1‑91 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      <  1:1149374/91‑1 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:1175184/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:2955121/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:407111/1‑91 (MQ=255)
                                                                                          aaaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc  >  1:1483274/1‑91 (MQ=255)
                                                                                          |                                                                                          
CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102

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: