Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 12,017 C→T 12.5% 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,0170CT12.5% 3518.3 / inf 1266intergenic (+180/+43)USA300HOU_pUSA300HOUMR0014/USA300HOU_pUSA300HOUMR0015recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (598/508);  new base T (79/79);  total (678/588)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 1276 total aligned reads displayed.

TAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11927‑12107
                                                                                          |                                                                                          
 aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAATAAAGCGTACAAATTCg                                                                                           >  2:3028808/1‑91 (MQ=255)
 aaaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCg                                                                                           >  2:3431492/1‑91 (MQ=255)
  aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                          >  2:3102569/1‑91 (MQ=255)
  aaaGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGt                                                                                          <  2:3537715/91‑1 (MQ=255)
    agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                        <  1:3059548/91‑1 (MQ=255)
    agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                        >  1:508527/1‑91 (MQ=255)
    agagACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTGTAATAATAACACAAGAAAAAAGCGTACAAATTCGTAc                                                                                        <  1:3336431/91‑1 (MQ=255)
     gagaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACg                                                                                       <  1:1725087/91‑1 (MQ=255)
      agaCGTTCTGACAGCGTCTAAAATCTCATGATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGc                                                                                      <  2:153946/91‑1 (MQ=38)
       gaCGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGCt                                                                                     <  2:1825000/91‑1 (MQ=255)
         cGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCttt                                                                                   >  1:1114286/1‑91 (MQ=255)
          gTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCtttt                                                                                  >  2:422042/1‑91 (MQ=255)
              tGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACGAATTTGTACGCttttttct                                                                              <  2:1602647/91‑3 (MQ=38)
              tGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaa                                                                              >  1:2547095/1‑91 (MQ=255)
                aCAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaata                                                                            >  1:3838749/1‑91 (MQ=255)
                     gTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTtaatataata                                                                       >  1:466821/1‑91 (MQ=255)
                      tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      <  1:3717987/91‑1 (MQ=255)
                      tCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAc                                                                      <  2:3724538/91‑1 (MQ=255)
                       ttaaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     >  1:2764686‑M2/38‑91 (MQ=255)
                       cTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATAcc                                                                     <  2:1044236/91‑1 (MQ=255)
                         aaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   >  1:641404‑M2/36‑91 (MQ=255)
                         aaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   >  2:3346740‑M2/36‑91 (MQ=255)
                         aaaatctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   >  1:393024‑M2/36‑91 (MQ=255)
                         AAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTg                                                                   <  1:684344‑M1/91‑1 (MQ=255)
                          aaatctcataacaggaaaatattatgaggttttaATAATAACACAAGCAAAAAGCGAAAAAATTCGTACGCTTTTTTCATATACTATCTGt                                                                  >  2:860743‑M2/35‑91 (MQ=255)
                           AATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGtt                                                                 <  2:2217126‑M1/91‑1 (MQ=255)
                             tctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTc                                                               >  2:705924‑M2/32‑91 (MQ=255)
                             TCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTc                                                               <  2:3147495‑M1/91‑1 (MQ=255)
                              ctcataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCa                                                              >  2:2269334‑M2/31‑91 (MQ=255)
                              CTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCa                                                              <  2:2849021‑M1/91‑1 (MQ=255)
                              CTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCa                                                              <  2:3831545‑M1/91‑1 (MQ=255)
                               TCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCAc                                                             <  2:148539‑M1/91‑1 (MQ=255)
                                 ataacaggaaaatattatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                           >  1:3602201‑M2/28‑91 (MQ=255)
                                 ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                           >  2:1716343‑M1/1‑91 (MQ=255)
                                 ATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTa                                                           <  2:619458‑M1/91‑1 (MQ=255)
                                  TAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  2:804347‑M1/1‑91 (MQ=255)
                                  TAATACTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTaa                                                          >  2:3646927‑M1/1‑91 (MQ=255)
                                   AATATTTTCCTGTTATGAGATTTTAATAATTACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAAt                                                         >  2:322448‑M1/1‑91 (MQ=255)
                                   AATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACCAATTCGTACGCTTTTTTAATACAATACCTGTTACCCTAAt                                                         <  2:1207056‑M1/91‑1 (MQ=255)
                                     tatTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                       <  1:3405410/91‑1 (MQ=255)
                                     tatTGTCCTTTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAGATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGa                                                       <  2:2340767/91‑1 (MQ=255)
                                        tttCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACg                                                    >  2:657776/1‑91 (MQ=255)
                                           ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                                 <  1:1295365/91‑1 (MQ=255)
                                           ccTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGttt                                                 <  1:2215241/91‑1 (MQ=255)
                                            cTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGtttt                                                <  1:3570188/91‑1 (MQ=255)
                                             tGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTa                                               >  1:557498/1‑91 (MQ=255)
                                               ttatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAg                                             >  1:1391570‑M2/14‑91 (MQ=255)
                                                tatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGt                                            >  2:496691‑M2/13‑91 (MQ=255)
                                                tatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGt                                            <  2:12950‑M2/79‑1 (MQ=255)
                                                tatgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGt                                            <  2:1131217‑M2/79‑1 (MQ=255)
                                                 atgagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTa                                           >  1:1542829‑M2/12‑91 (MQ=255)
                                                   gagattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATc                                         >  1:1046555‑M2/10‑91 (MQ=255)
                                                     gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                       <  1:3777232‑M2/84‑1 (MQ=255)
                                                     gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                       <  1:298643‑M2/84‑1 (MQ=255)
                                                     gattttaATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCtt                                       >  1:2706161‑M2/8‑91 (MQ=255)
                                                            ataataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCAtt                                >  2:2773649/1‑91 (MQ=255)
                                                             taataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTa                               >  1:262623/1‑91 (MQ=255)
                                                             taataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTa                               >  1:294270/1‑91 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  1:2151266/91‑1 (MQ=255)
                                                              aataaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTaa                              <  2:1769803/91‑1 (MQ=255)
                                                                 aaCACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTg                           <  2:1376936/91‑1 (MQ=255)
                                                                    acaAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGaaa                        >  1:2282010/1‑91 (MQ=255)
                                                                     caAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAAc                       <  1:2433832/91‑1 (MQ=255)
                                                                      aaGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACt                      >  1:2891150/1‑91 (MQ=255)
                                                                      aaGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGCTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACt                      >  2:2743346/1‑91 (MQ=255)
                                                                       aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                     <  2:1001037/91‑1 (MQ=255)
                                                                       aGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTa                     >  1:72283/1‑91 (MQ=255)
                                                                        gAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTAc                    >  1:1466964/1‑91 (MQ=255)
                                                                          aaaaaGCGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAg                  >  1:3263460/1‑91 (MQ=255)
                                                                          aaaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAg                  >  1:2509177/1‑91 (MQ=255)
                                                                          aaaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAg                  <  2:3604115/91‑1 (MQ=255)
                                                                           aaaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATTATTAACTGAAACTACAGt                 >  2:2917174/1‑91 (MQ=255)
                                                                            aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                <  2:2595093/91‑1 (MQ=255)
                                                                            aaaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGta                <  1:1926749/91‑1 (MQ=255)
                                                                             aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  1:1341646/1‑91 (MQ=255)
                                                                             aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  2:2673192/1‑91 (MQ=255)
                                                                             aaGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaa               >  2:2069576/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             >  1:283092/1‑91 (MQ=255)
                                                                               gCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaata             <  1:1834483/91‑1 (MQ=255)
                                                                                cGTACGAATTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGtaataa            >  1:2033169/1‑91 (MQ=255)
                                                                                 gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           <  1:1254844/91‑1 (MQ=255)
                                                                                 gTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTaataat           >  2:1997255/1‑91 (MQ=255)
                                                                                   acgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  1:3500835/4‑91 (MQ=255)
                                                                                   acgaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  1:790915/4‑91 (MQ=255)
                                                                                   aCAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  2:2126760/1‑91 (MQ=255)
                                                                                   aCAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         >  1:1794758/1‑91 (MQ=255)
                                                                                   aCAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATAc         <  1:117964/91‑1 (MQ=255)
                                                                                    cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  1:1007023/1‑91 (MQ=255)
                                                                                    cAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACt        >  2:329477/1‑91 (MQ=255)
                                                                                     gaaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc       >  1:2360448/2‑91 (MQ=255)
                                                                                     aaaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTc       <  2:2417765/91‑1 (MQ=255)
                                                                                      aaTTTGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      >  2:1678531/1‑91 (MQ=255)
                                                                                      aaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      <  1:2085981/91‑1 (MQ=255)
                                                                                      aaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      >  1:1645217/1‑91 (MQ=255)
                                                                                      aaTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      >  2:1914983/1‑91 (MQ=255)
                                                                                      aaTTCGTACGCTTTTTTAATATAATACCTGTTTCACCAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCg      <  1:870420/91‑1 (MQ=255)
                                                                                       aTTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGt     <  2:938395/91‑1 (MQ=255)
                                                                                        ttCGTACGCTTTTTTAAGATAATACCTGTTTCACTAATGAACGTTTTAAGTAGCTTCATCATTAACTGAAACTACAGTAATAATACTCGTc    <  2:249698/91‑1 (MQ=255)
                                                                                         tCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTca   <  1:3180154/91‑1 (MQ=255)
                                                                                          cGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTTATCATTAACTGAAACTACAGTAATAATACTCGTcat  >  2:3917119/1‑90 (MQ=255)
                                                                                          |                                                                                          
TAAAAGAGACGTTCTGACAGCGTCTAAAATCTCATAATATTTTCCTGTTATGAGATTTTAATAATAACACAAGAAAAAAGCGTACAAATTCGTACGCTTTTTTAATATAATACCTGTTTCACTAATGAACGTTTTAAGTATCTTCATCATTAACTGAAACTACAGTAATAATACTCGTCAA  >  CP000731/11927‑12107

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: