Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,920:1 +T 7.9% intergenic (+343/‑47) 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,9201.T 7.9% 807.7 / 60.9 202intergenic (+343/‑47)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (81/105);  new base T (11/5);  total (92/110)
Fisher's exact test for biased strand distribution p-value = 6.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 211 total aligned reads displayed.

GATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGT  >  USA300TCH1516_ALE/576830‑577009
                                                                                           |                                                                                         
gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                          <  1:975487/91‑1 (MQ=255)
gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                          <  1:2001498/91‑1 (MQ=255)
 aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTAt                                                                                          >  2:3530207/1‑91 (MQ=255)
 aTAATTTGGAAGCGCCAGCAATATTATTAATGTTAGTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGGAtt                                                                                         <  2:1233005/91‑2 (MQ=37)
  tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATg                                                                                         <  2:1122907/91‑1 (MQ=255)
  tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATg                                                                                         >  2:2783130/1‑91 (MQ=255)
      ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                                     >  1:308843/1‑91 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTTTGTTAgt                                                                                    >  1:1943570/1‑91 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                    >  1:1679010/1‑91 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                    >  1:2423613/1‑91 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACGTGTATGTTAgt                                                                                    <  2:482053/91‑1 (MQ=255)
        ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   <  1:255445/91‑1 (MQ=255)
         gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  >  1:1346391/1‑91 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                 >  1:2895973/1‑91 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                                 <  1:921732/91‑5 (MQ=21)
           aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                >  2:1367473/1‑91 (MQ=255)
           aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                <  2:643553/91‑1 (MQ=255)
              gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                             >  2:1997862/1‑91 (MQ=255)
               ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttt                                                                            <  1:2781163/91‑1 (MQ=255)
               ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttt                                                                            >  1:1965916/1‑91 (MQ=255)
                cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCttcc                                                                           <  2:1179764/91‑3 (MQ=17)
                     atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                      <  1:2335261/91‑1 (MQ=255)
                     atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                      <  1:2335420/91‑1 (MQ=255)
                       attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                    <  1:2817584/91‑1 (MQ=255)
                        ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAtt                                                                   >  2:2262904/1‑91 (MQ=255)
                          attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTa                                                                 >  1:2072284/1‑91 (MQ=255)
                          attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTa                                                                 <  1:1438081/91‑1 (MQ=255)
                           ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                                <  1:1301224/91‑1 (MQ=255)
                             aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAtt                                                              >  1:318873/1‑91 (MQ=255)
                              aatgtATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCAAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                             <  2:3181888/87‑1 (MQ=255)
                               tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                            <  1:1997418/91‑1 (MQ=255)
                               tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                            >  1:70002/1‑91 (MQ=255)
                               tGTCATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                            <  1:3494685/91‑1 (MQ=255)
                                    ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       <  1:2264044/91‑1 (MQ=255)
                                    ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       <  1:3777481/91‑1 (MQ=255)
                                    ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       <  1:152220/91‑1 (MQ=255)
                                    ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       >  2:3109855/1‑91 (MQ=255)
                                     tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                                      >  1:3815753/1‑91 (MQ=255)
                                      cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                                >  2:1762714/1‑80 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                    >  2:1000075/1‑91 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                    <  1:1543177/91‑1 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                    >  2:3816857/1‑91 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATACTTAATTTCGCAGGTAc                                                    >  2:3714402/1‑91 (MQ=255)
                                         tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACtt                                                  >  2:2071032/1‑91 (MQ=255)
                                          aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                                 <  2:1725775/91‑1 (MQ=255)
                                           aaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCa                                                <  1:3562684/91‑1 (MQ=255)
                                            aTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCaa                                               >  2:2556739/1‑91 (MQ=255)
                                             ttCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAt                                              <  1:2880079/91‑1 (MQ=255)
                                               cAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTa                                            >  2:3282456/1‑91 (MQ=255)
                                                agagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTAt                                           >  2:1834737/1‑91 (MQ=255)
                                                 gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                          >  1:3689742/1‑91 (MQ=255)
                                                  agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                         >  1:1052749/1‑91 (MQ=255)
                                                  agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                         <  1:595591/91‑1 (MQ=255)
                                                    tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                                       >  1:1139661/1‑91 (MQ=255)
                                                    tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                                       <  2:3338227/91‑1 (MQ=255)
                                                       aaCCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTa                                    <  1:1243383/91‑1 (MQ=255)
                                                       aaCCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTa                                    >  1:2031218/1‑91 (MQ=255)
                                                         ccGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACg                                  >  2:1183884/1‑91 (MQ=255)
                                                          cGTGCGGAAGCGCCTCTAACAATACAGCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                                 <  2:1960403/91‑1 (MQ=255)
                                                          cGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                                 <  1:1397088/91‑1 (MQ=255)
                                                          cGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                                 >  2:2367661/1‑91 (MQ=255)
                                                           gTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTa                                >  2:338211/1‑91 (MQ=255)
                                                            tGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                               <  2:1522816/91‑1 (MQ=255)
                                                              cGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGtt                             >  2:1082101/1‑91 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           <  2:2531385/91‑1 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           <  2:723702/91‑1 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           >  2:28137/1‑91 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           <  1:2231391/91‑1 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           >  2:1535597/1‑91 (MQ=255)
                                                                 aaGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCaa                          <  1:3426629/91‑1 (MQ=255)
                                                                  aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                         <  2:1205679/91‑1 (MQ=255)
                                                                  aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                         <  2:299681/91‑1 (MQ=255)
                                                                         ctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag                  <  2:1711465/91‑1 (MQ=255)
                                                                          tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg                 >  1:1847370/1‑91 (MQ=255)
                                                                          tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg                 >  1:946374/1‑91 (MQ=255)
                                                                            aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag               <  2:1688427/91‑1 (MQ=255)
                                                                            aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag               <  1:2373459/91‑1 (MQ=255)
                                                                            aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag               <  1:2151357/91‑1 (MQ=255)
                                                                            aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag               <  1:274529/91‑1 (MQ=255)
                                                                             cAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGaggagg              <  1:946987/91‑1 (MQ=255)
                                                                              aaTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga             <  1:1022618/91‑1 (MQ=255)
                                                                              aaTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga             <  1:2451395/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)
                                                                                tACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTAAGTAGTTCAATCATTGAGGAGGaa            >  2:3341013/1‑91 (MQ=255)
                                                                                tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa           <  2:3182319/91‑1 (MQ=255)
                                                                                 aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          >  2:3356395/1‑91 (MQ=255)
                                                                                 aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          <  2:3591675/91‑1 (MQ=255)
                                                                                 aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          >  2:2278678/1‑91 (MQ=255)
                                                                                  cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg         <  1:2248604/91‑1 (MQ=255)
                                                                                  cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg         >  1:2644110/1‑91 (MQ=255)
                                                                                    aCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGAt       <  1:3260002/91‑1 (MQ=255)
                                                                                      ttGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCGGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCt     <  1:3643893/91‑1 (MQ=255)
                                                                                      ttGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCt     >  2:1285107/1‑91 (MQ=255)
                                                                                       tGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCtt    >  1:534831/1‑91 (MQ=255)
                                                                                       tGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCtt    >  2:392421/1‑91 (MQ=255)
                                                                                       tGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCtt    >  2:439738/1‑91 (MQ=255)
                                                                                        gTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTg   >  2:3814282/1‑91 (MQ=255)
                                                                                         gaTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGt  >  1:175012/2‑91 (MQ=255)
                                                                                           |                                                                                         
GATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGT  >  USA300TCH1516_ALE/576830‑577009

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: