Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 803,528 C→T 13.0% intergenic (‑228/+42) USA300HOU_RS03920 ← / ← dtpT Putative lipid kinase/Di‑/tripeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE803,5280CT13.0% 618.6 / 63.6 229intergenic (‑228/+42)USA300HOU_RS03920/dtpTPutative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand):  ref base C (102/96);  new base T (17/13);  total (119/110)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 241 total aligned reads displayed.

AGTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAT  >  USA300TCH1516_ALE/803438‑803618
                                                                                          |                                                                                          
aGTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGc                                                                                            <  2:3023258/91‑1 (MQ=255)
 gTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGcc                                                                                           >  1:510782/1‑91 (MQ=255)
 gTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGcc                                                                                           >  1:460008/1‑91 (MQ=255)
  ttATAGCTGTCTATAATGAGCATATTAAGTTTTTTTTTATACTGATATCTTGAATTTAATTTATAGAAACCTATAAAAAAACAGTAAGTCa                                                                                          >  2:2509282/1‑91 (MQ=255)
  ttATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          <  1:3251896/91‑1 (MQ=255)
    ataGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCAtt                                                                                        <  1:2265435/91‑1 (MQ=255)
     taGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCAttt                                                                                       <  2:2423474/91‑1 (MQ=255)
      aGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTa                                                                                      >  2:3353090/1‑91 (MQ=255)
       gCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTaa                                                                                     >  2:1121893/1‑91 (MQ=255)
        cTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTaaa                                                                                    >  2:1702586/1‑91 (MQ=255)
        cTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTaaa                                                                                    >  2:2176032/1‑91 (MQ=255)
         tGTCTATAATGAGCATATTAAGTTTTGATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAAt                                                                                   >  2:173771/1‑91 (MQ=255)
             tataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAAAAAAATATAAAAAAACAGTAAGCCATTTAAATGACt                                                                               <  2:55164/91‑1 (MQ=255)
               taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                             <  1:923445/91‑1 (MQ=255)
               taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                             >  1:2414493/1‑91 (MQ=255)
                aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                            <  1:1091071/91‑1 (MQ=255)
                 aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACt                                                                           >  1:2241191/1‑91 (MQ=255)
                  tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                          >  2:2931053/1‑91 (MQ=255)
                  tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                          >  2:2931408/1‑91 (MQ=255)
                   gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                         >  2:1958068/1‑91 (MQ=255)
                   gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                         <  2:1139685/91‑1 (MQ=255)
                    aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTACTAGAAACCTATAAAAATCCAGTAAGCCATGTAAATGACTCACTGtt                                                                        >  2:2439442/1‑91 (MQ=255)
                    aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                        <  2:3170856/91‑1 (MQ=255)
                    aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTCACTGtt                                                                        <  1:3022635/91‑1 (MQ=38)
                    aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtt                                                                        <  1:3815657/91‑1 (MQ=255)
                     gCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGttt                                                                       <  1:1546844/91‑1 (MQ=255)
                      cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                      <  2:219947/91‑1 (MQ=255)
                       atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttttt                                                                     >  1:2501159/1‑91 (MQ=255)
                        tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                    <  2:1924861/91‑1 (MQ=255)
                        tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                    <  1:149502/91‑1 (MQ=255)
                         atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                      >  1:1329385/1‑88 (MQ=255)
                          ttAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGa                                                                  >  2:3592008/1‑91 (MQ=255)
                           tAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                 <  1:3155505/91‑1 (MQ=255)
                           tAAGTTTTTATTTATACTGATATCTTGAATTCAATTAATAGAAACCTATAAAAAAACCGCACGCCATTTAAAGGACTTACTGTTTTTTGaa                                                                 >  2:1243280/1‑91 (MQ=255)
                            aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAt                                                                <  2:1187753/91‑1 (MQ=255)
                             aGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAtt                                                               >  2:2097937/1‑91 (MQ=255)
                              gTTTTTATTTATACTGATATCTTGAGTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTa                                                              >  1:2003538/1‑91 (MQ=255)
                                     tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                       >  1:141537/1‑91 (MQ=255)
                                     tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                       <  1:3868249/91‑1 (MQ=255)
                                     tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                       <  2:266881/91‑1 (MQ=255)
                                      ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                      <  1:3208763/91‑1 (MQ=255)
                                      ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                      <  2:811478/91‑1 (MQ=255)
                                       tataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaacaa                                                     >  1:635739/1‑91 (MQ=255)
                                         taCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAata                                                   >  2:1157964/1‑91 (MQ=255)
                                         taCTGATATCGTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAGTGACTTACTGTTTTTTGAATTAGGCCAACAata                                                   <  2:299641/91‑1 (MQ=255)
                                           cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATAtt                                                 >  2:3736216/1‑91 (MQ=255)
                                             gATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTaa                                               >  1:1579539/1‑91 (MQ=255)
                                              atatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAAc                                              <  1:355649/91‑1 (MQ=255)
                                                atCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGt                                            <  2:460008/91‑1 (MQ=255)
                                                 tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                           <  2:1484300/91‑1 (MQ=255)
                                                 tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                           >  1:2537153/1‑91 (MQ=255)
                                                 tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                           >  1:2218229/1‑91 (MQ=255)
                                                 tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTAACTGTGTATTGAATTAGGCCAACAATAATAACCta                                           >  2:737879/1‑91 (MQ=255)
                                                  cTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCGATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtat                                          <  2:3398502/91‑1 (MQ=255)
                                                   ttGACTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATAACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                         >  2:3413212/1‑91 (MQ=255)
                                                   ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGGATTAGGCCAACAATATTAACGtata                                         >  2:3250787/1‑91 (MQ=255)
                                                   ttGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATAACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtata                                         >  2:3412720/1‑91 (MQ=255)
                                                    tGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATAc                                        >  2:400908/1‑91 (MQ=255)
                                                     gAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATAcc                                       >  1:2639405/1‑91 (MQ=255)
                                                       aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                     >  1:1162772/1‑91 (MQ=255)
                                                          taattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCa                                  <  2:845553/91‑1 (MQ=255)
                                                          taattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCa                                  >  2:1007203/1‑91 (MQ=255)
                                                              taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                              <  1:2176032/91‑1 (MQ=255)
                                                              taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                              >  2:1881818/1‑91 (MQ=255)
                                                              taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                              >  2:2715999/1‑91 (MQ=255)
                                                               aatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCt                             <  1:2714018/91‑1 (MQ=255)
                                                                atAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATcaaag                            <  1:3350916/91‑5 (MQ=255)
                                                                  aGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtg                          >  1:2950887/1‑91 (MQ=255)
                                                                   gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                         <  1:2809710/91‑1 (MQ=255)
                                                                   gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                         <  1:2132149/91‑1 (MQ=255)
                                                                   gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                         <  2:2413130/91‑1 (MQ=255)
                                                                   gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                         <  2:1588704/91‑1 (MQ=255)
                                                                     aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                       >  1:3088216/1‑91 (MQ=255)
                                                                      aCCTATAAAAAAACAGTAAGCCATTTAAATGAGTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTCTCATCGCTTtgatga                      <  1:3313923/91‑1 (MQ=255)
                                                                         tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                   >  2:2614995/1‑91 (MQ=255)
                                                                         tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                   >  1:3697158/1‑91 (MQ=255)
                                                                         tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                   >  2:140080/1‑91 (MQ=255)
                                                                          ataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaa                  >  1:1705847/1‑91 (MQ=255)
                                                                           taAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaaa                 <  1:2029139/91‑1 (MQ=255)
                                                                            aaaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                >  1:3544223/1‑91 (MQ=255)
                                                                             aaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg               <  1:3783638/91‑1 (MQ=255)
                                                                              aaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt              >  2:1042295/1‑91 (MQ=255)
                                                                                aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga            <  1:1776142/91‑1 (MQ=255)
                                                                                aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga            <  2:2351113/91‑1 (MQ=255)
                                                                                 aaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgag           <  1:812222/91‑1 (MQ=255)
                                                                                  aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGGCAACAATATTAACGTATACCTTTCATCGCTTTGTTGATTAAAGGTgaga          >  2:2101982/1‑91 (MQ=255)
                                                                                  aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgaga          >  1:680725/1‑91 (MQ=255)
                                                                                     gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg       >  1:1161782/1‑91 (MQ=255)
                                                                                     gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg       >  2:2305584/1‑91 (MQ=255)
                                                                                     gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTCGATGATTAAAGGTGAGATTg       >  2:3183561/1‑91 (MQ=255)
                                                                                     gTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg       <  2:894369/91‑1 (MQ=255)
                                                                                      tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc      >  1:3451669/1‑91 (MQ=255)
                                                                                      tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc      >  2:637616/1‑91 (MQ=255)
                                                                                      tAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc      >  1:649061/1‑91 (MQ=255)
                                                                                       aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     >  1:66571/1‑91 (MQ=255)
                                                                                       aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     <  2:654694/91‑1 (MQ=255)
                                                                                       aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     <  1:3438971/91‑1 (MQ=255)
                                                                                       aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     <  1:303775/91‑1 (MQ=255)
                                                                                         gCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa   <  2:3471778/91‑1 (MQ=255)
                                                                                         gCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa   >  1:2079577/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAT  >  USA300TCH1516_ALE/803438‑803618

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: