Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,012 A→G 12.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,0120AG12.0% 3319.5 / inf 1230intergenic (+175/+48)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (594/483);  new base G (29/120);  total (624/606)
Fisher's exact test for biased strand distribution p-value = 5.64e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 1267 total aligned reads displayed.

CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102
                                                                                          |                                                                                          
cGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACa                                                                                            >  1:1798566/1‑91 (MQ=255)
   tGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAt                                                                                         <  1:2919113/91‑1 (MQ=255)
   tGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAAt                                                                                         >  1:438751/1‑91 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtt                                                                                        <  2:183944/91‑1 (MQ=255)
    gTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtt                                                                                        <  2:2430825/91‑1 (MQ=255)
     tAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAttt                                                                                       <  1:1758130/91‑2 (MQ=255)
     tAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTc                                                                                       >  1:2595874/1‑91 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:1274333/91‑4 (MQ=255)
       aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAAtttgt                                                                                     <  1:743474/91‑4 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                   <  2:3821650/91‑1 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                   >  1:508527/1‑91 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTGTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                   <  1:3336431/91‑1 (MQ=255)
         agagACGTTCTGACAGCGTCTAAAATCGCATAATAGTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                   >  2:179163/1‑91 (MQ=255)
          tataCGTTATGACAGATTATAAAATCTAATAATATTTTCCTTTTATTATATTTTAATAATAACACAAGAAAAAAGCGTACAAATTAGTAct                                                                                  <  2:1420178/88‑2 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  2:486865/91‑1 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                  >  1:2077815/1‑91 (MQ=255)
          gagaCGTTCTGACAGCGTCTAAAATCTCATAATATGTTCCTGTTATGAGATTTTAATAATAGCACAAGAAAAAAGCGTACGAATTTGTACg                                                                                  <  2:789987/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                 <  1:1121664/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                 <  1:2331739/91‑1 (MQ=255)
           agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATGCGTACGc                                                                                 <  2:2694171/91‑1 (MQ=255)
                ttCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGCttttt                                                                            <  2:2082694/91‑1 (MQ=255)
                ttCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCttttt                                                                            >  1:1313880/1‑91 (MQ=255)
                   tGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaa                                                                         <  2:1169918/91‑1 (MQ=255)
                        gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    <  2:1686688/91‑1 (MQ=255)
                        gCGTCTAAAATCTCATAATATGTTCCTGTTATGAGATTTTAATAATAACACAAGAAAACAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    <  2:1340977/91‑1 (MQ=255)
                        attattaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                    >  1:2491269‑M2/42‑91 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 <  2:3724538/91‑1 (MQ=255)
                           tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 <  2:2854738/91‑1 (MQ=255)
                           attaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAACGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                 >  2:3205767‑M2/39‑91 (MQ=255)
                            ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                >  2:2150085‑M2/38‑91 (MQ=255)
                            cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                <  1:142050/91‑1 (MQ=255)
                             TAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               >  1:823102‑M1/1‑91 (MQ=255)
                             TAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCt                                                               <  2:275343‑M1/91‑1 (MQ=255)
                              aaaatctcataagaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              >  1:3036908‑M2/36‑91 (MQ=255)
                              AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              >  2:227293‑M1/1‑91 (MQ=255)
                              AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              >  2:3109650‑M1/1‑91 (MQ=255)
                              AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              <  2:2227889‑M1/91‑1 (MQ=255)
                              AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAACAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                              <  2:40533‑M1/91‑1 (MQ=255)
                               AAATCTCATAATATTTTCCTGTTATTTGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                             >  2:3418209‑M1/1‑91 (MQ=255)
                                AATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGtt                                                            <  2:3599505‑M1/91‑1 (MQ=255)
                                 ATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                           <  1:3516438‑M1/91‑1 (MQ=255)
                                 ATCTCATAATAGTGTCCTGTGATGGGATTTGAATAATAACACAAGAAAAAAGCGTACAAAGTCGTACGCTTTTGTAATATAATACCtgtgt                                                           >  2:712758‑M1/1‑89 (MQ=255)
                                    TCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                        >  1:2693074‑M1/1‑91 (MQ=255)
                                    TCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                        <  2:2932365‑M1/91‑1 (MQ=255)
                                      ataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                      >  1:3602201‑M2/28‑91 (MQ=255)
                                      ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAGGCTTTTTTAATATATTACCTGTTTCACTa                                                      <  1:1364427‑M1/91‑1 (MQ=255)
                                      ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                      >  1:1187387‑M1/1‑91 (MQ=255)
                                      ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                      <  1:1671534‑M1/91‑1 (MQ=255)
                                       taacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     >  2:162149‑M2/27‑91 (MQ=255)
                                       TAATATTTTCGTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     <  2:15146‑M1/91‑1 (MQ=255)
                                       TAATATTTTCCTGTTATGAGATTTTAATAGTCACACAAGAAAAAAGCGTACAACTTCATACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     >  1:3164112‑M1/1‑91 (MQ=255)
                                       TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     <  1:2871106‑M1/91‑1 (MQ=255)
                                       TAATACTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                     >  2:3646927‑M1/1‑91 (MQ=255)
                                            ttttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAAc                                                >  1:652239/1‑91 (MQ=255)
                                            ttttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAAc                                                <  1:2983586/91‑1 (MQ=255)
                                             tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                               >  2:657776/1‑91 (MQ=255)
                                              ttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                              <  1:1048833/91‑1 (MQ=255)
                                              ttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                              >  2:2621354/1‑91 (MQ=255)
                                              ttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                              <  2:1064682/91‑1 (MQ=255)
                                                ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                            <  1:2215241/91‑1 (MQ=255)
                                                 cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                           <  1:3570188/91‑1 (MQ=255)
                                                    ttatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAg                                        <  2:1267897‑M2/78‑1 (MQ=255)
                                                    ttatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAg                                        <  1:2753598‑M2/78‑1 (MQ=255)
                                                    ttatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACCAATGCACGTCTTAAg                                        >  2:1553822‑M2/14‑91 (MQ=255)
                                                     tatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGt                                       <  1:3413470‑M2/79‑1 (MQ=255)
                                                             tttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcat                               >  1:2033264‑M2/5‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAGGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  1:1948258‑M2/2‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  2:3139198‑M2/90‑1 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  2:3287855‑M2/2‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTcatcat                            >  1:947778‑M2/2‑91 (MQ=255)
                                                                aATAATAACACAAGAAAAAAGCGTACAAATGCGGACGCTTTTTTAATATAATACCTGTGTCACTAATGAACGTTTTAAGTATCTTcatcat                            <  2:934648‑M2/90‑1 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACACATTCGTCCGCTTTTTTAATATAATACCTGTTTCACTAGTGAACGTTTTAAGTATCTTCATCAtt                           >  2:202063/1‑91 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           >  1:1587108/1‑91 (MQ=255)
                                                                 ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                           >  2:3915355/1‑91 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  2:2356133/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:2151266/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         >  1:2794911/1‑91 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:3824815/91‑1 (MQ=255)
                                                                   aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                         <  1:563403/91‑1 (MQ=255)
                                                                     tatcacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACt                       <  2:2525791/88‑1 (MQ=255)
                                                                     taaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACt                       >  1:3907742/1‑91 (MQ=255)
                                                                      aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                      >  2:3518016/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  1:3770955/1‑91 (MQ=255)
                                                                         acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                   >  1:2282010/1‑91 (MQ=255)
                                                                          caAGAAAAAGGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTGTTAAGTATCTTCATCATTAACTGAAAc                  <  2:1016360/91‑1 (MQ=255)
                                                                            aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                >  2:2791559/1‑91 (MQ=255)
                                                                              aaaaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACa              <  1:2542076/91‑1 (MQ=255)
                                                                               aaaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAg             >  1:3263460/1‑91 (MQ=255)
                                                                                aaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGt            >  1:186630/1‑91 (MQ=255)
                                                                                    gCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  2:852572/1‑91 (MQ=255)
                                                                                    gCGTACATATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        <  2:604229/91‑1 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTCTCTTCAACATTAACTGAAACTACAGtaata        >  2:3031781/1‑91 (MQ=255)
                                                                                    gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata        >  2:1944602/1‑91 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat      <  1:3703726/91‑1 (MQ=255)
                                                                                      gTACAAATTCGTACGCTTTTTTAATATAAAACCTGATTCACTATTGAACGGTTTAAGTATCTTCATTATTAACTGCAACTACAGTaataat      >  2:3337558/1‑91 (MQ=255)
                                                                                       tACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAataata     >  1:713388/1‑91 (MQ=255)
                                                                                       tACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCATCATCATTAACTGAAACTACAGTAataata     >  2:1724812/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:1007023/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:2429138/1‑91 (MQ=255)
                                                                                         cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt   >  1:1512739/1‑91 (MQ=255)
                                                                                          |                                                                                          
CGATGTAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTC  >  CP000731/11922‑12102

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