Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,915 C→G 11.6% intergenic (+338/‑52) gltX → / → cysE Glutamate‑‑tRNA ligase/Serine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE576,9150CG11.6% 659.6 / 63.5 207intergenic (+338/‑52)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (77/104);  new base G (13/11);  total (90/117)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 209 total aligned reads displayed.

CTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATC  >  USA300TCH1516_ALE/576827‑577005
                                                                                        |                                                                                           
cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                           <  1:3836040/91‑1 (MQ=255)
 ttGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGt                                                                                         >  1:1179764/1‑91 (MQ=38)
  tGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGt                                                                                         <  1:893901/91‑1 (MQ=255)
  tGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAAGTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGt                                                                                         <  2:2352638/91‑1 (MQ=255)
   gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                      <  1:975487/91‑1 (MQ=255)
   gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                      <  1:2001498/91‑1 (MQ=255)
    aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTAt                                                                                      >  2:3530207/1‑91 (MQ=255)
    aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTAt                                                                                      >  1:2256346/1‑91 (MQ=255)
    aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAtt                                                                                     <  1:2124677/91‑2 (MQ=37)
    aTAATTTGGAAGCGCCAGCAATATTATTAATGTTAGTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGGAtt                                                                                     <  2:1233005/91‑2 (MQ=37)
     tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATg                                                                                     <  2:1122907/91‑1 (MQ=255)
         ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                                 >  1:308843/1‑91 (MQ=255)
          tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTTTGTTAgt                                                                                >  1:1943570/1‑91 (MQ=255)
          tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                >  1:2423613/1‑91 (MQ=255)
          tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                >  1:1679010/1‑91 (MQ=255)
          tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACGTGTATGTTAgt                                                                                <  2:482053/91‑1 (MQ=255)
            gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAGCAATACAACTTGTATGTTAgtgg                                                                              <  2:1415331/91‑1 (MQ=255)
             aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             >  2:1527783/1‑91 (MQ=255)
             aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             >  1:2895973/1‑91 (MQ=255)
             aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                             <  1:921732/91‑5 (MQ=21)
             aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                             <  1:363549/91‑5 (MQ=21)
              aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                            <  2:643553/91‑1 (MQ=255)
              aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                            <  1:145338/91‑1 (MQ=255)
              aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                            >  1:3212963/1‑91 (MQ=255)
                cgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCt                                                                          >  2:2025188/1‑91 (MQ=255)
                 gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                         >  2:1720223/1‑91 (MQ=255)
                 gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                         >  2:1997690/1‑91 (MQ=255)
                  ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttt                                                                        >  1:1965916/1‑91 (MQ=255)
                  ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttt                                                                        <  1:2781113/91‑1 (MQ=255)
                   cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCttcc                                                                       <  2:1179764/91‑3 (MQ=17)
                        atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                  <  2:2160939/91‑1 (MQ=255)
                        atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                  >  1:395772/1‑91 (MQ=255)
                        atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                  <  1:2335420/91‑1 (MQ=255)
                        atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                  <  1:2335261/91‑1 (MQ=255)
                             attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTa                                                             >  1:2072284/1‑91 (MQ=255)
                              ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTTCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                            <  1:2630976/91‑1 (MQ=255)
                                  tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                        >  1:70002/1‑91 (MQ=255)
                                  tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                        <  1:1997418/91‑1 (MQ=255)
                                   gTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCg                                                       <  2:3181315/91‑1 (MQ=255)
                                       ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                   <  1:3777481/91‑1 (MQ=255)
                                       ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                   >  2:428398/1‑91 (MQ=255)
                                       ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                   >  2:3109855/1‑91 (MQ=255)
                                         cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTa                                                 >  1:1628668/1‑91 (MQ=255)
                                          gatAAATTCAGAGTGGAACCGTGCGTAAGCGCCTCTAACAATACAACTTGTATGTTAGTTGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                <  2:264786/89‑1 (MQ=255)
                                          gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                >  2:3816857/1‑91 (MQ=255)
                                          gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                >  2:1000075/1‑91 (MQ=255)
                                             aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                             >  1:1298872/1‑91 (MQ=255)
                                                ttCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAt                                          <  1:2880079/91‑1 (MQ=255)
                                                 tCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAtt                                         >  2:1347896/1‑91 (MQ=255)
                                                 tCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAtt                                         <  1:2179008/91‑1 (MQ=255)
                                                    gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                      <  1:2812229/91‑1 (MQ=255)
                                                    gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                      >  1:3689742/1‑91 (MQ=255)
                                                    gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                      >  2:3260002/1‑91 (MQ=255)
                                                     agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                     >  1:1052749/1‑91 (MQ=255)
                                                     agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                     <  1:595591/91‑1 (MQ=255)
                                                     agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                     >  1:1529027/1‑91 (MQ=255)
                                                       tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                                   <  2:3338353/91‑1 (MQ=255)
                                                       tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                                   >  1:1139661/1‑91 (MQ=255)
                                                        ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  1:3180700/91‑1 (MQ=255)
                                                        ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  1:2692263/91‑1 (MQ=255)
                                                         gAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAttt                                 <  1:1310189/91‑1 (MQ=255)
                                                         gAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAttt                                 <  2:1798623/91‑1 (MQ=255)
                                                          aaCCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTa                                >  1:2031218/1‑91 (MQ=255)
                                                          aaCCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAATTCAATTATGTATTTa                                >  2:3704093/1‑91 (MQ=255)
                                                            ccGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACg                              >  2:1183884/1‑91 (MQ=255)
                                                            ccGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACg                              <  1:3517262/91‑1 (MQ=255)
                                                             cGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                             >  2:2367661/1‑91 (MQ=255)
                                                              gTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTa                            >  1:3425998/1‑91 (MQ=255)
                                                              gTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTa                            >  2:338211/1‑91 (MQ=255)
                                                               tGCGGAAGCGCCTCTATCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                           <  2:3715674/91‑1 (MQ=255)
                                                               tGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                           <  2:1757337/91‑1 (MQ=255)
                                                                gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                          <  2:2794/91‑1 (MQ=255)
                                                                 cGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGtt                         >  2:1082101/1‑91 (MQ=255)
                                                                  ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                        <  1:2518515/91‑1 (MQ=255)
                                                                  ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                        <  1:1347896/91‑1 (MQ=255)
                                                                  ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                        >  2:2265608/1‑91 (MQ=255)
                                                                   gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                       <  2:723702/91‑1 (MQ=255)
                                                                   gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                       >  2:28137/1‑91 (MQ=255)
                                                                   gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                       <  2:2531385/91‑1 (MQ=255)
                                                                   gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                       >  2:1535597/1‑91 (MQ=255)
                                                                     aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                     <  2:299681/91‑1 (MQ=255)
                                                                     aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                     <  2:1205679/91‑1 (MQ=255)
                                                                       cgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTACTTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCa                   >  2:1884183/1‑91 (MQ=255)
                                                                          ctctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTg                >  1:3476091/1‑91 (MQ=255)
                                                                          cactAACATACAAGTTGTATTGTTAGAGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTg                >  2:3348753/3‑91 (MQ=25)
                                                                           actAACATACAAGTTGTATTGTTAGAGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTga               >  1:224525/2‑91 (MQ=37)
                                                                             tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg             >  1:946374/1‑91 (MQ=255)
                                                                             tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg             >  1:1847370/1‑91 (MQ=255)
                                                                               aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag           <  2:1688427/91‑1 (MQ=255)
                                                                               aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag           <  1:2373459/91‑1 (MQ=255)
                                                                                cAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGaggagg          <  1:946987/91‑1 (MQ=255)
                                                                                 aaTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga         <  1:1022618/91‑1 (MQ=255)
                                                                                  aTACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga         >  1:2407749/1‑91 (MQ=255)
                                                                                  aTACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga         >  2:3873495/1‑91 (MQ=255)
                                                                                  aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa        <  1:886031/91‑1 (MQ=255)
                                                                                   tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa       <  2:3181982/91‑1 (MQ=255)
                                                                                    aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt      <  2:3591675/91‑1 (MQ=255)
                                                                                     cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCGTTGAGGAGGAAATg     <  2:3293090/91‑1 (MQ=255)
                                                                                        cTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATc  <  1:2447630/91‑1 (MQ=255)
                                                                                        cTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATc  <  1:1427062/91‑1 (MQ=255)
                                                                                        |                                                                                           
CTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATC  >  USA300TCH1516_ALE/576827‑577005

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