Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,211,658 G→A 6.5% intergenic (+651/+162) USA300HOU_RS11275 → / ← yidC hypothetical protein/Membrane protein insertase YidC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,211,6580GA 6.5% 698.7 / 23.2 211intergenic (+651/+162)USA300HOU_RS11275/yidChypothetical protein/Membrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base G (101/96);  new base A (7/7);  total (108/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 219 total aligned reads displayed.

GACATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAA  >  USA300TCH1516_ALE/2211568‑2211748
                                                                                          |                                                                                           
gACATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAg                                                                                             >  2:2372561/1‑91 (MQ=255)
gACATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAg                                                                                             >  2:2833339/1‑91 (MQ=255)
gACATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAg                                                                                             <  1:3127034/91‑1 (MQ=255)
  cATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGa                                                                                           >  1:1021063/1‑91 (MQ=255)
    ttGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATg                                                                                         >  2:307939/1‑91 (MQ=255)
    ttGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATg                                                                                         <  2:3832348/91‑1 (MQ=255)
    ttGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATg                                                                                         >  2:493983/1‑91 (MQ=255)
     tGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGTATGa                                                                                        >  1:533401/1‑91 (MQ=255)
      gCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGAc                                                                                       >  2:2172772/1‑91 (MQ=255)
       cGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACt                                                                                      <  1:3918880/91‑1 (MQ=255)
       cGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACt                                                                                      <  1:3372512/91‑1 (MQ=255)
        gATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACtt                                                                                     <  2:3133995/91‑1 (MQ=255)
        gATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACtt                                                                                     <  2:115225/91‑1 (MQ=255)
         aTGAAGAATGACAACTCTATTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACttt                                                                                    <  2:1391176/91‑1 (MQ=255)
           gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                  >  1:945120/1‑91 (MQ=255)
           gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                  <  1:1265168/91‑1 (MQ=255)
           gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                  >  1:2752809/1‑91 (MQ=255)
           gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                  >  2:1990859/1‑91 (MQ=255)
             agaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTcaca                                                                                >  2:1353900/1‑91 (MQ=255)
             agaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTcaca                                                                                <  2:205871/91‑1 (MQ=255)
              gaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAt                                                                               <  1:3779285/91‑1 (MQ=255)
               aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                              <  1:1573574/91‑1 (MQ=255)
               aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                              >  2:2864583/1‑91 (MQ=255)
                aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                             <  1:1990339/91‑1 (MQ=255)
                aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                             <  2:2618816/91‑1 (MQ=255)
                 tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                            >  1:293928/1‑91 (MQ=255)
                 tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                            >  1:1217986/1‑91 (MQ=255)
                 tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                            <  1:3440194/91‑1 (MQ=255)
                  gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                           >  1:624288/1‑91 (MQ=255)
                   aCAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaa                                                                          <  1:3326281/91‑1 (MQ=255)
                    cAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaaa                                                                         >  1:3491773/1‑91 (MQ=255)
                       ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                     >  2:2112277/1‑91 (MQ=255)
                       ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                     >  2:1519977/1‑91 (MQ=255)
                       ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                     >  2:518902/1‑91 (MQ=255)
                         ctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGtcatc                                                                   <  2:2332198/91‑1 (MQ=21)
                          tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCt                                                                  >  1:996754/1‑91 (MQ=255)
                          tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGtcatcc                                                                  <  1:229421/91‑2 (MQ=21)
                                taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCGAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                            <  2:3485628/91‑1 (MQ=255)
                                taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                            <  2:2781689/91‑1 (MQ=255)
                                   ccACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGt                                                         >  1:1720526/1‑91 (MQ=255)
                                   ccACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGt                                                         >  2:1211834/1‑91 (MQ=255)
                                    cACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGtt                                                        <  2:683022/91‑1 (MQ=255)
                                     aCTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttt                                                       <  1:304257/91‑1 (MQ=255)
                                      cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                      <  1:985375/91‑1 (MQ=255)
                                      cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                      <  2:1189645/91‑1 (MQ=255)
                                       tttttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCt                                                     >  1:567655/1‑91 (MQ=255)
                                          ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                  >  2:1680168/1‑91 (MQ=255)
                                          ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                  <  2:158450/91‑1 (MQ=255)
                                             cATACTGAAAGCACGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTAt                                               <  2:1237509/91‑1 (MQ=255)
                                             cATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTAt                                               >  1:1690566/1‑91 (MQ=255)
                                              aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                              <  2:708263/91‑1 (MQ=255)
                                              aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                              <  2:2406724/91‑1 (MQ=255)
                                                 cTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcaca                                           <  2:2563653/91‑1 (MQ=255)
                                                 cTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcaca                                           >  1:1808503/1‑91 (MQ=255)
                                                   gAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACata                                         >  2:2501469/1‑91 (MQ=255)
                                                     aaaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATAtt                                       <  2:3712232/91‑1 (MQ=255)
                                                      aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTg                                      <  2:2730562/91‑1 (MQ=255)
                                                      aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTg                                      <  1:1245909/91‑1 (MQ=255)
                                                       aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                     <  1:370999/91‑1 (MQ=255)
                                                        cTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGta                                    <  1:1604078/91‑1 (MQ=255)
                                                         tCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtat                                   >  2:84206/1‑91 (MQ=255)
                                                         tCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtat                                   <  1:1864282/91‑1 (MQ=255)
                                                           gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                 <  2:3273732/91‑1 (MQ=255)
                                                           gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                 >  2:1816971/1‑91 (MQ=255)
                                                            ttAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAt                                >  2:310355/1‑91 (MQ=255)
                                                            ttAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAt                                >  2:3310597/1‑91 (MQ=255)
                                                               atatTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgtg                             <  1:3907303/91‑1 (MQ=255)
                                                                 attattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATCTTGTATAATGTGac                           <  2:1857544/91‑1 (MQ=255)
                                                                 attattTCAAATAAAAACAGCAGTAGGATGACTTTCACATCTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGac                           <  2:1447387/91‑1 (MQ=255)
                                                                  ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                          >  1:74612/1‑91 (MQ=255)
                                                                    attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                        >  1:3580974/1‑91 (MQ=255)
                                                                     tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                       <  2:3922191/91‑1 (MQ=255)
                                                                     tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                       <  2:2354220/91‑1 (MQ=255)
                                                                      ttCAGATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                      <  2:2689392/91‑1 (MQ=255)
                                                                      ttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                      <  1:2209828/91‑1 (MQ=255)
                                                                      ttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                      <  1:1597155/91‑1 (MQ=255)
                                                                        cAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGt                    <  2:1754574/91‑1 (MQ=255)
                                                                          aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                  >  2:1179118/1‑91 (MQ=255)
                                                                            tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                <  2:922843/91‑1 (MQ=255)
                                                                            tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                >  1:2465573/1‑91 (MQ=255)
                                                                             agaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGc               >  1:2339208/3‑91 (MQ=11)
                                                                              aaaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt              >  1:73374/1‑91 (MQ=255)
                                                                               aaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTa             >  1:2633553/1‑91 (MQ=255)
                                                                               aaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTa             >  2:727967/1‑91 (MQ=14)
                                                                                aaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCGCATATTGTATAATGTGACACTAAGTTTCGCTAt            >  1:432940/1‑91 (MQ=255)
                                                                                 aCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTAtt           >  1:3216965/1‑91 (MQ=255)
                                                                                  cagcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCttcaa          <  1:1136824/91‑5 (MQ=255)
                                                                                  cagcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTg          >  1:283664/1‑91 (MQ=255)
                                                                                   agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa         >  1:2074786/1‑91 (MQ=255)
                                                                                   agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa         >  1:2558705/1‑91 (MQ=255)
                                                                                    gcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa        >  1:2892630/1‑91 (MQ=255)
                                                                                    gcagTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa        >  2:1667687/1‑91 (MQ=17)
                                                                                     cagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAg       >  2:2149931/1‑91 (MQ=255)
                                                                                      agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc      <  1:1040281/91‑1 (MQ=255)
                                                                                      agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc      >  2:1889431/1‑91 (MQ=255)
                                                                                       gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg     <  1:3572804/91‑1 (MQ=255)
                                                                                        tAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGa    >  2:2928951/1‑91 (MQ=255)
                                                                                         aGGATGACTTTCACAGTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaa   <  2:532464/91‑1 (MQ=255)
                                                                                          ggATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaaa  >  1:1871280/1‑91 (MQ=255)
                                                                                          ggATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTCAAGCGaaa  >  2:3443339/1‑91 (MQ=255)
                                                                                          |                                                                                           
GACATTGCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAA  >  USA300TCH1516_ALE/2211568‑2211748

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: