Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,019,747 T→C 5.2% intergenic (+69/+44) USA300HOU_RS05030 → / ← ltaA Putative phosphoesterase/putative glycolipid permease LtaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,019,7470TC 5.2% 829.6 / 24.4 253intergenic (+69/+44)USA300HOU_RS05030/ltaAPutative phosphoesterase/putative glycolipid permease LtaA
Reads supporting (aligned to +/- strand):  ref base T (112/126);  new base C (4/9);  total (118/135)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 270 total aligned reads displayed.

TGAAACGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATAAAATA  >  USA300TCH1516_ALE/1019657‑1019837
                                                                                          |                                                                                          
 gAAACGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGta                                                                                           >  1:3754064/1‑91 (MQ=255)
 gAAACGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGta                                                                                           >  2:23235/1‑91 (MQ=255)
  aaaCGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGtaa                                                                                          <  2:63692/91‑1 (MQ=255)
  aaaCGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGtaa                                                                                          <  1:902665/91‑1 (MQ=255)
       tttAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATaaaa                                                                                     <  1:2309738/91‑1 (MQ=255)
         tAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAGTAAAAAt                                                                                   >  2:3502890/1‑91 (MQ=255)
          aaaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATg                                                                                  >  2:3378105/1‑91 (MQ=255)
          aaaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATg                                                                                  <  1:1432473/91‑1 (MQ=255)
           aaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGa                                                                                 >  2:529738/1‑91 (MQ=255)
           aaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGa                                                                                 >  1:1872344/1‑91 (MQ=255)
           aaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGa                                                                                 <  1:3323932/91‑1 (MQ=255)
           aaTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGa                                                                                 >  1:2524053/1‑91 (MQ=255)
            aTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGAc                                                                                <  1:3823876/91‑1 (MQ=255)
            aTTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTACATTCAATAAAAATGAc                                                                                <  1:966810/91‑1 (MQ=25)
              tAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTACATTCAATAAAAATGACCa                                                                              <  2:713820/91‑1 (MQ=25)
               aGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCAt                                                                             >  2:3149120/1‑91 (MQ=255)
                 cTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCAttt                                                                           <  1:1103720/91‑1 (MQ=255)
                  ttAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCAtttt                                                                          <  2:330681/91‑1 (MQ=255)
                  ttAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTACATTCAATAAAAATGACCAtttt                                                                          <  2:30970/91‑1 (MQ=21)
                    aaGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCt                                                                        >  1:2386922/1‑91 (MQ=255)
                    aaGTAACATAATAGTATTGTTAATCGTAGTACGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAACGACCATTTTCt                                                                        >  2:3018583/1‑91 (MQ=255)
                     aGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCtt                                                                       >  2:1881249/1‑91 (MQ=255)
                     aGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGATCATTTTTATTGAATGTAATAAAAATGACCATTTTCtt                                                                       >  2:2727959/1‑91 (MQ=255)
                      gTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCttt                                                                      <  2:321146/91‑1 (MQ=255)
                        aaCATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTAt                                                                    >  1:1318814/1‑91 (MQ=255)
                          cATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTAttt                                                                  >  1:2315205/1‑91 (MQ=255)
                          cATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTAttt                                                                  >  1:3433006/1‑91 (MQ=255)
                           ataataGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTAtttt                                                                 <  1:755373/91‑1 (MQ=255)
                            taataGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTa                                                                >  1:1752352/1‑91 (MQ=255)
                             aataGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTaa                                                               <  2:946584/91‑1 (MQ=255)
                             aataGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTaa                                                               >  1:816094/1‑91 (MQ=255)
                              ataGTATTGTTAATCGTAGTATGTTTCAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTCCTTTATTTTaaa                                                              >  2:3281457/1‑91 (MQ=255)
                               taGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTaaaa                                                             <  1:3410842/91‑1 (MQ=255)
                                aGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAAt                                                            <  2:3697385/91‑1 (MQ=255)
                                  tATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATAc                                                          <  1:2645064/91‑1 (MQ=255)
                                   aTTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACg                                                         >  1:1580546/1‑91 (MQ=255)
                                    ttGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGt                                                        >  1:2220149/1‑91 (MQ=255)
                                    ttGTGAATCGTAGTATGTTGGAATTAATAAGAAAATGTTCAGTTTCCGTGAATGGAATCCAAATGACAATGTGCGTTATTTCAAACGACGt                                                        >  2:3365614/1‑91 (MQ=255)
                                     tGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGtt                                                       <  1:1500017/91‑1 (MQ=255)
                                      gTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGttt                                                      <  2:1530503/91‑1 (MQ=255)
                                      gTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGttt                                                      >  1:1158583/1‑91 (MQ=255)
                                      gTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGttt                                                      >  1:854903/1‑91 (MQ=255)
                                         aaTCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTaa                                                   >  1:572154/1‑91 (MQ=255)
                                         aaTCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTaa                                                   >  2:2417373/1‑91 (MQ=255)
                                         aaTCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTaa                                                   <  1:1619827/91‑1 (MQ=255)
                                          aTCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAAc                                                  <  2:1238372/91‑1 (MQ=255)
                                           tCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATTACCATTTTCTTTATTTTAAAATACGTTTTAAcc                                                 >  1:3234772/1‑91 (MQ=255)
                                           tCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAAcc                                                 <  1:3705027/91‑1 (MQ=255)
                                           tCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAAcc                                                 <  1:1946950/91‑1 (MQ=255)
                                           tCGTAGTATGTTTGAATTAATAAGAAAATGGTCAATTTTATTGAATGTAATAAAAAGGACCATTTTCTTTATTTTAAAATACGTTTTAAcc                                                 <  2:3797361/91‑1 (MQ=255)
                                              tagtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACctta                                              <  1:456335/91‑1 (MQ=255)
                                               agtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttac                                             >  2:2006106/1‑91 (MQ=255)
                                               agtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttac                                             >  1:1717623/1‑91 (MQ=255)
                                               agtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttac                                             <  1:3468039/91‑1 (MQ=255)
                                                gtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttact                                            <  1:64403/91‑1 (MQ=255)
                                                gtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttact                                            >  2:2873026/1‑91 (MQ=255)
                                                gtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttact                                            <  1:1014558/91‑1 (MQ=255)
                                                gtaTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttact                                            >  2:2635118/1‑91 (MQ=255)
                                                 taTGTTTGGATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAGATGACCATTTTCTTTATTTTAAAATACGTTTTATCcttactt                                           <  2:443242/91‑1 (MQ=255)
                                                 taTGTTTGAATTAATAAGAGAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttactt                                           <  2:2800986/91‑1 (MQ=255)
                                                 taTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttactt                                           <  2:173393/91‑1 (MQ=255)
                                                  aTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttactta                                          <  2:458377/91‑1 (MQ=255)
                                                  aTGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACcttactta                                          <  2:2213988/91‑1 (MQ=255)
                                                     tttGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCt                                       <  2:2748142/91‑1 (MQ=255)
                                                      ttGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCtt                                      <  1:3024876/91‑1 (MQ=255)
                                                       tGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCttt                                     <  2:1676064/91‑1 (MQ=255)
                                                       tGAATTAATAAGAAAATGGTCATGATTATTGAATGTAATAAAAATGACCAGTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCttt                                     <  2:3671663/91‑1 (MQ=255)
                                                        gAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCtttt                                    <  2:2574358/91‑1 (MQ=255)
                                                         aaTTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTATATTACGTTTTAACCTTACTTAGCttttt                                   >  2:1836049/1‑91 (MQ=255)
                                                          aTTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTtc                                  <  2:900004/91‑1 (MQ=255)
                                                          aTTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTtc                                  <  2:2146241/91‑1 (MQ=255)
                                                          aTTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTtc                                  <  2:597938/91‑1 (MQ=255)
                                                              ataaGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTa                              >  2:486489/1‑91 (MQ=255)
                                                                aaGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTAtt                            <  2:1012368/91‑1 (MQ=255)
                                                                aaGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTAtt                            <  1:706727/91‑1 (MQ=255)
                                                                 aGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTAttt                           <  2:821088/91‑1 (MQ=255)
                                                                  gAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTa                          <  2:1938830/91‑1 (MQ=255)
                                                                  gAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTa                          >  1:2788531/1‑91 (MQ=255)
                                                                  gAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTa                          >  2:2681946/1‑91 (MQ=255)
                                                                  gAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTa                          >  1:2686125/1‑91 (MQ=255)
                                                                     aaTGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACta                       >  1:294884/1‑91 (MQ=255)
                                                                       tGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACtata                     <  1:887355/91‑1 (MQ=255)
                                                                        ggTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATaa                    <  1:3178561/91‑1 (MQ=255)
                                                                          tCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAg                  >  2:3328061/1‑91 (MQ=255)
                                                                          tCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAg                  >  1:1809098/1‑91 (MQ=255)
                                                                          tCATTTTTATTACATTCAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAg                  >  1:3613814/1‑91 (MQ=25)
                                                                              ttttATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGc              <  1:2122338/91‑1 (MQ=255)
                                                                               tttATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCt             <  1:2055788/91‑1 (MQ=255)
                                                                               tttATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCt             <  2:1926085/91‑1 (MQ=255)
                                                                                  aTTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTcc          <  1:3357953/91‑1 (MQ=255)
                                                                                   ttGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCa         >  1:308449/1‑91 (MQ=255)
                                                                                     gAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATa       >  2:2762992/1‑91 (MQ=255)
                                                                                     gAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATa       >  1:3723483/1‑91 (MQ=255)
                                                                                     gAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCAATATTTACTAGAAAGTAGCTTCCATa       >  2:3420340/1‑91 (MQ=255)
                                                                                      aaTGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATaa      >  1:3083146/1‑91 (MQ=255)
                                                                                       aTGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATaaa     >  2:3721529/1‑91 (MQ=255)
                                                                                        tGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATaaaa    <  2:161475/91‑1 (MQ=255)
                                                                                        tGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATaaaa    <  1:4158/91‑1 (MQ=255)
                                                                                         gTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATAAAAt   >  1:2877329/1‑91 (MQ=255)
                                                                                         gTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATAAAAt   >  2:1590322/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGAAACGTTTAAATTAGCTTAAGTAACATAATAGTATTGTTAATCGTAGTATGTTTGAATTAATAAGAAAATGGTCATTTTTATTGAATGTAATAAAAATGACCATTTTCTTTATTTTAAAATACGTTTTAACCTTACTTAGCTTTTTCTCTATTTACTATAAAGTAGCTTCCATAAAATA  >  USA300TCH1516_ALE/1019657‑1019837

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