Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 24,137 G→A 7.2% intergenic (+389/‑41) purA → / → USA300TCH1516_00018 Adenylosuccinate synthetase/tRNA‑Glu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE24,1370GA 7.2% 711.2 / 27.1 224intergenic (+389/‑41)purA/USA300TCH1516_00018Adenylosuccinate synthetase/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base G (83/124);  new base A (11/6);  total (94/130)
Fisher's exact test for biased strand distribution p-value = 7.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

AATGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAACAC  >  USA300TCH1516_ALE/24047‑24226
                                                                                          |                                                                                         
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                           <  2:900854/91‑1 (MQ=255)
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                           <  1:2495194/91‑1 (MQ=255)
  tGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCa                                                                                         <  2:1199975/91‑1 (MQ=255)
  tGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCGTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGCCATGACa                                                                                         <  2:583839/91‑1 (MQ=255)
   gAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCaa                                                                                        <  2:2268105/91‑1 (MQ=255)
     aTAGCCTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGCCATGACAAGa                                                                                      <  1:3296287/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAACATATCTTGTCATGGCAAGa                                                                                      >  2:2439886/1‑91 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGa                                                                                      <  2:3153218/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGa                                                                                      <  1:498261/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGa                                                                                      <  1:2580650/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGCCATGACAAGa                                                                                      <  2:3800627/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGCCATGACAAGa                                                                                      <  2:2740143/91‑1 (MQ=255)
      tAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGAt                                                                                     <  2:1609496/91‑1 (MQ=255)
      tAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGAt                                                                                     <  1:641656/91‑1 (MQ=255)
       aGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGAtt                                                                                    >  2:1232139/1‑91 (MQ=255)
        gACTGCCAATAAAATGCAATTTTTCACTTATAACATTATTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTg                                                                                   <  2:1915330/91‑1 (MQ=255)
           tGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCt                                                                                >  2:296868/1‑91 (MQ=255)
            gCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtg                                                                               >  2:3221185/1‑91 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                             <  2:3262486/91‑1 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                             <  2:748220/91‑1 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                             <  2:2820386/91‑1 (MQ=255)
               aaTATAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTATGc                                                                            >  1:3117981/1‑91 (MQ=255)
               aaTAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGc                                                                            <  1:2385099/91‑1 (MQ=255)
               aaTAAAATGCAATTTATCACTTATAACATTCTTCAAAAAATAATAGCAAAAAAAAAAAAAAAATATCTTGTCATGGCAAGATTGGCTGTGc                                                                            <  2:1287770/91‑1 (MQ=255)
                aTAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCt                                                                           <  1:892312/91‑1 (MQ=255)
                 tAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCta                                                                          <  1:2338603/91‑1 (MQ=255)
                  aaaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtat                                                                         >  2:2782699/1‑91 (MQ=255)
                  aaaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtat                                                                         >  1:296797/1‑91 (MQ=255)
                   aaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtata                                                                        <  1:3069316/91‑1 (MQ=255)
                     aTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAat                                                                      >  1:1277475/1‑91 (MQ=255)
                      tGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatc                                                                     >  1:839166/1‑91 (MQ=255)
                      tGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatc                                                                     <  1:2174978/91‑1 (MQ=255)
                      tGCAATTTTTCACTTATAACATTATTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatc                                                                     <  2:1553242/91‑1 (MQ=255)
                       gCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatct                                                                    <  2:1064867/91‑1 (MQ=255)
                       gCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatct                                                                    <  1:1137322/91‑1 (MQ=255)
                       gCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatct                                                                    <  1:1758632/91‑1 (MQ=255)
                          agttttCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctatc                                                                 <  2:1022255/89‑1 (MQ=255)
                          aTTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctatc                                                                 <  2:2920857/91‑1 (MQ=255)
                          aTTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctatc                                                                 >  2:1501091/1‑91 (MQ=255)
                           tttttCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctatct                                                                <  2:1742550/91‑1 (MQ=255)
                            ttttCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCtt                                                               <  2:220157/91‑1 (MQ=255)
                             tttCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTtg                                                              >  1:3703333/1‑91 (MQ=255)
                               tCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTtgtg                                                            <  1:1065962/91‑1 (MQ=255)
                               tCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTtgtg                                                            >  1:1634093/1‑91 (MQ=255)
                                cACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGc                                                           <  2:2867802/91‑1 (MQ=255)
                                 aCTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCt                                                          <  2:3257913/91‑1 (MQ=255)
                                     ataACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAg                                                      <  1:2860057/91‑1 (MQ=255)
                                     ataACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAg                                                      <  1:2858831/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAACATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                     >  1:2365237/1‑91 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                     <  2:1356053/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                     <  1:1745378/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                     <  2:247061/91‑1 (MQ=255)
                                       aaCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGaa                                                    >  2:676536/1‑91 (MQ=255)
                                        aCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAAc                                                   <  1:547269/91‑1 (MQ=255)
                                        aCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAAc                                                   <  1:327890/91‑1 (MQ=255)
                                            tcttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCt                                               >  2:1664278/1‑91 (MQ=255)
                                             cttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTc                                              <  1:2084727/91‑1 (MQ=255)
                                             cttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTc                                              <  2:609418/91‑1 (MQ=255)
                                                 aaaaaaTAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTg                                          >  1:3891283/1‑91 (MQ=255)
                                                  aaaaaTAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTgg                                         <  2:1714704/91‑1 (MQ=255)
                                                     aataatAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCa                                      <  1:1542195/91‑1 (MQ=255)
                                                             cAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTa                              <  2:1614410/91‑1 (MQ=255)
                                                              aaaaTTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTaa                             <  2:3504400/91‑1 (MQ=255)
                                                              aaaaTTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTaa                             <  2:2584534/91‑1 (MQ=255)
                                                              aaaaTTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTaa                             >  1:741765/1‑91 (MQ=255)
                                                               aaaTTATGTAAAGAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAg                            >  1:3031387/1‑91 (MQ=255)
                                                                   tATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGacac                        >  1:3642074/1‑91 (MQ=255)
                                                                    aTGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACAcc                       <  1:3775404/91‑1 (MQ=255)
                                                                      gTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGc                     >  2:272405/1‑91 (MQ=255)
                                                                        aaaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccc                   >  1:1666227/1‑91 (MQ=255)
                                                                         aaagaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCt                  <  1:1328672/91‑1 (MQ=255)
                                                                         aaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCt                  <  1:80882/91‑1 (MQ=255)
                                                                           aaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccgtga                <  1:1090086/91‑5 (MQ=255)
                                                                           aaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccgtga                <  1:2483551/91‑5 (MQ=255)
                                                                             atatCTTGTCATGGCAAGATTGGCTGTGCTAGAATCGAGCTGGTGCTTAAGAACGGCTCCTGGGTCAAGCGGTTAAGACACCGCCCTTTCa              <  2:3850643/91‑1 (MQ=37)
                                                                               atCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcg            <  1:1367544/91‑1 (MQ=255)
                                                                               atCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcg            <  2:3568434/91‑1 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCCAGCGGTTAAGACACCGCCCTTTCAcgg           >  1:385168/1‑91 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcgg           <  2:3621035/91‑1 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcgg           >  1:854875/1‑91 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcgg           <  1:1809672/91‑1 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcgg           <  1:714292/91‑1 (MQ=255)
                                                                                 cTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggc          <  1:310170/91‑1 (MQ=255)
                                                                                 cTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggc          >  1:1012326/1‑91 (MQ=255)
                                                                                 cTTGCCATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggc          >  2:204216/1‑91 (MQ=255)
                                                                                  ttgccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcg         >  1:1111256/5‑91 (MQ=255)
                                                                                  ttGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcg         >  1:394899/1‑91 (MQ=255)
                                                                                  ttGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGATTAAGAACGGCTCCTTGGTCAAGCGGGTAAGACACCGCCCTTTCAcggcg         >  2:3680857/1‑91 (MQ=255)
                                                                                   tgccATGACAAGATTGGCTGTGCTATCATCTATCTTGTGCTTAAGAACGGCACCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcgg        >  2:3142773/4‑91 (MQ=25)
                                                                                   tgccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGGCAAGCGGTTAAGACACCGCCCTTTCAcggcgg        >  1:3378876/4‑91 (MQ=37)
                                                                                   tGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcgg        <  1:1803649/91‑1 (MQ=255)
                                                                                    gccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGt       >  2:1090484/3‑91 (MQ=255)
                                                                                    gTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGt       <  2:3063427/91‑1 (MQ=255)
                                                                                    gTCATGGCAAGATTGGCTGTGCGATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCTCCTTCCcggcggg       >  2:3148189/1‑90 (MQ=37)
                                                                                     ccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTa      >  2:319493/2‑91 (MQ=255)
                                                                                     ccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCATTCACGGCGGTa      >  1:2020289/2‑91 (MQ=255)
                                                                                      cATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTaa     <  1:366577/91‑1 (MQ=255)
                                                                                        tGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAaca   <  2:648466/91‑1 (MQ=255)
                                                                                        tGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAaca   >  1:1409346/1‑91 (MQ=255)
                                                                                         ggCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAacac  <  1:1030809/91‑1 (MQ=255)
                                                                                          |                                                                                         
AATGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAACAC  >  USA300TCH1516_ALE/24047‑24226

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: