Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,017 C→T 11.0% 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,0170CT11.0% 1922.2 / 163.7 669intergenic (+180/+43)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (307/286);  new base T (28/45);  total (337/332)
Fisher's exact test for biased strand distribution p-value = 3.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 675 total aligned reads displayed.

TAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11927‑12107
                                                                                          |                                                                                          
 aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                           >  2:664480/1‑91 (MQ=255)
 aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                           <  1:1812925/91‑1 (MQ=255)
  aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTACTAATAACACAAGAAGAAAGCGTACAAATTCGt                                                                                          >  1:3316850/1‑91 (MQ=255)
  aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                          >  2:2013224/1‑91 (MQ=255)
   aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTa                                                                                         <  1:1013313/91‑1 (MQ=255)
   aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTa                                                                                         <  2:3472570/91‑1 (MQ=255)
   aaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTa                                                                                         >  1:3186176/1‑91 (MQ=255)
    agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTAc                                                                                        <  1:1724926/91‑1 (MQ=255)
     gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                       >  2:2257103/1‑91 (MQ=255)
     gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                       >  2:3144255/1‑91 (MQ=255)
     gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                       <  1:2106884/91‑1 (MQ=255)
     gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATCTTAATAATAACAAAAGAAAAAAGCGTACACATTCGTACg                                                                                       >  2:3642039/1‑91 (MQ=255)
      agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                      <  1:2895199/91‑1 (MQ=255)
      agaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                      <  1:3851732/91‑1 (MQ=255)
      agaCGTTCTGACAGCGTCTAAAATCTCATAATATTATCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGc                                                                                      <  2:1511025/91‑1 (MQ=255)
         cGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCttt                                                                                   >  1:1945909/1‑91 (MQ=255)
           ttCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTATTAATAACACAAGAAAAAAGCGTACATATTCGTACGCttttt                                                                                 >  2:373671/1‑91 (MQ=255)
            tCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttttt                                                                                >  1:604428/1‑91 (MQ=255)
             cTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTta                                                                               <  2:3829818/91‑1 (MQ=255)
               gACTGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                             <  1:3131799/91‑1 (MQ=255)
               gACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaat                                                                             <  2:1833359/91‑1 (MQ=255)
                   gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                         <  2:1146771/91‑1 (MQ=255)
                   gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                         >  2:2665148/1‑91 (MQ=255)
                   gCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataa                                                                         >  2:799429/1‑91 (MQ=255)
                    ttattaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataat                                                                        >  1:36677‑M2/41‑91 (MQ=255)
                     gTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataata                                                                       >  2:3882053/1‑91 (MQ=255)
                      tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      >  2:1616438/1‑91 (MQ=255)
                      tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      >  2:1872697/1‑91 (MQ=255)
                      tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      <  1:1049946/91‑1 (MQ=255)
                         AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   <  1:2600597‑M1/91‑1 (MQ=255)
                          aaatctcataacaggaaaatattatgagatttttATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGt                                                                  >  2:454156‑M2/35‑91 (MQ=255)
                           AATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGtt                                                                 <  1:3575304‑M1/91‑1 (MQ=255)
                            ATCTCATAAGATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGttt                                                                >  1:3721965‑M1/1‑91 (MQ=255)
                             TCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTc                                                               <  2:650523‑M1/91‑1 (MQ=255)
                               tcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                             >  2:500640‑M2/30‑91 (MQ=255)
                                cataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                            >  2:2553814‑M2/29‑91 (MQ=255)
                                CATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACt                                                            >  2:3037874‑M1/1‑91 (MQ=255)
                                 ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                           >  1:1523014‑M1/1‑91 (MQ=255)
                                  TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  2:1425637‑M1/1‑91 (MQ=255)
                                  TAATAGTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  2:3584276‑M1/1‑91 (MQ=255)
                                    ATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCCGTTTCACTAATg                                                        >  2:277038‑M1/1‑91 (MQ=255)
                                     tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                       >  2:2991910/1‑91 (MQ=255)
                                     tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                       <  1:1640961/91‑1 (MQ=255)
                                       ttgtCCTGTTATGCGATTTTAATAATAACACAAGAAAAAAGCTTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAAc                                                     <  2:730586/88‑1 (MQ=255)
                                        tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                                    >  2:2344964/1‑91 (MQ=255)
                                        tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                                    >  1:1244495/1‑91 (MQ=255)
                                        tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                                    <  2:3427195/91‑1 (MQ=255)
                                         ttCCTGTTATGAGATTTTAATAATAACACACGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGt                                                   <  1:3766362/91‑1 (MQ=255)
                                          tCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtt                                                  >  1:1895546/1‑91 (MQ=255)
                                           ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                                 >  1:723101/1‑91 (MQ=255)
                                            atattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                                >  1:271267‑M2/17‑91 (MQ=255)
                                              gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTACTGAACGTTTTaa                                              >  2:2690113/1‑91 (MQ=255)
                                              gTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTaa                                              >  1:2337757/1‑91 (MQ=255)
                                                  tgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTAt                                          <  2:3011845‑M2/81‑1 (MQ=255)
                                                      attttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtc                                      <  1:28462‑M2/85‑1 (MQ=255)
                                                       ttttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                     <  2:238280‑M2/86‑1 (MQ=255)
                                                       ttgtaATGATAACGCAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATATGACCTGTTTCACTAATGAACGTTTTAAGTATCTtca                                     <  2:908055‑M2/86‑1 (MQ=255)
                                                          taATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatca                                  <  2:3244000‑M2/89‑1 (MQ=255)
                                                          taATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatca                                  >  1:3832939‑M2/3‑91 (MQ=255)
                                                          taATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTtcatca                                  >  2:2598743‑M2/3‑91 (MQ=255)
                                                             taataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTa                               >  1:3633215/1‑91 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  2:3697488/91‑1 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              >  2:1737421/1‑91 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  1:3455379/91‑1 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  1:1587431/91‑1 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAAGTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              >  2:2205140/1‑91 (MQ=255)
                                                               ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                             <  1:2389844/91‑1 (MQ=255)
                                                               ataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAAc                             <  1:526584/91‑1 (MQ=255)
                                                                 aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                           <  1:941158/91‑1 (MQ=255)
                                                                   cacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                         >  2:2085896/1‑91 (MQ=255)
                                                                   cacaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaa                         >  2:568749/1‑91 (MQ=255)
                                                                     caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                       <  1:2069260/91‑1 (MQ=255)
                                                                       aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                     >  1:3073236/1‑91 (MQ=255)
                                                                        gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc                    >  2:1965499/1‑91 (MQ=255)
                                                                         aaaaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACa                   <  1:899977/91‑1 (MQ=255)
                                                                           aaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGt                 <  1:2499307/91‑1 (MQ=255)
                                                                            aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                <  1:2540332/91‑1 (MQ=255)
                                                                            aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                <  1:2257103/91‑1 (MQ=255)
                                                                            aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACGGAAACTACAGta                <  1:2997422/91‑1 (MQ=255)
                                                                             aaGCGTACAAATTCGTACGCTTTCTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  2:3513084/1‑91 (MQ=255)
                                                                              aGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaat              >  1:1998949/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  2:104296/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             <  2:1248567/91‑1 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  2:3320888/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  2:3359217/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  1:134681/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             <  1:2294779/91‑1 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             <  2:1471483/91‑1 (MQ=255)
                                                                                 gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           >  1:1426982/1‑91 (MQ=255)
                                                                                 gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           >  1:1089489/1‑91 (MQ=255)
                                                                                  tACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAataata          >  1:2807420/1‑91 (MQ=255)
                                                                                  tACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAataata          >  1:3874367/1‑91 (MQ=255)
                                                                                   acgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  2:2128463/4‑91 (MQ=255)
                                                                                   aCAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  2:1351299/1‑91 (MQ=255)
                                                                                    cgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:2448362/3‑91 (MQ=255)
                                                                                    cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        <  1:3372627/91‑1 (MQ=255)
                                                                                    cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:986143/1‑91 (MQ=255)
                                                                                      aaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      >  1:2654211/1‑91 (MQ=255)
                                                                                         tCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTca   <  2:2887710/91‑1 (MQ=255)
                                                                                         tCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTca   <  1:454158/91‑1 (MQ=255)
                                                                                          |                                                                                          
TAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11927‑12107

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: