Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 890,548 A→T 11.3% intergenic (+23/‑41) gcvH → / → USA300TCH1516_00820 Glycine cleavage system H protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE890,5480AT11.3% 748.2 / 68.2 246intergenic (+23/‑41)gcvH/USA300TCH1516_00820Glycine cleavage system H protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (116/102);  new base T (9/19);  total (125/121)
Fisher's exact test for biased strand distribution p-value = 4.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTACATTTTA  >  USA300TCH1516_ALE/890458‑890638
                                                                                          |                                                                                          
ttAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTa                                                                                            >  1:2529095/1‑91 (MQ=255)
ttAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTa                                                                                            <  2:1052181/91‑1 (MQ=255)
 tAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTTGTGAATAATCACCGTGTAACTCCTTAATCTaa                                                                                           <  2:3815066/91‑1 (MQ=255)
  aGTGATGAAAGTCAGATTGAAGCTTTATCAACAGCTGTGAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTCATCTAAg                                                                                          >  2:2437129/1‑91 (MQ=255)
   gTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGa                                                                                         <  2:1217229/91‑1 (MQ=255)
   gTGATGAAAGCCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGa                                                                                         <  1:2743980/91‑1 (MQ=255)
    tgatgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAt                                                                                        <  1:1063012/91‑1 (MQ=255)
       tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                                     <  2:2642117/91‑3 (MQ=255)
       tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                                     <  1:736099/91‑3 (MQ=255)
       tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                                     <  1:2838299/91‑3 (MQ=255)
        gaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgag                                                                                    <  1:2694842/91‑1 (MQ=255)
        gaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaag                                                                                    <  2:2239878/91‑4 (MQ=255)
         aaaGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGATTaagg                                                                                   <  2:2023949/91‑1 (MQ=255)
          aaGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgagga                                                                                  <  1:2761213/91‑1 (MQ=255)
            gTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGt                                                                                <  2:460371/91‑1 (MQ=255)
              cAGCTTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGttt                                                                              <  1:3497128/91‑1 (MQ=255)
              cAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGttt                                                                              >  1:1620337/1‑91 (MQ=255)
               aGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTa                                                                             >  1:389268/1‑91 (MQ=255)
                gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                            >  2:2190291/1‑91 (MQ=255)
                gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                            <  2:890543/91‑1 (MQ=255)
                 aTTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAtt                                                                           >  1:951430/1‑91 (MQ=255)
                 aTTGAAGCTTTAATAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAtt                                                                           >  1:2990176/1‑91 (MQ=255)
                  ttGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGATTAttt                                                                          >  2:40841/1‑91 (MQ=255)
                    gAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttttt                                                                        <  1:2934184/91‑1 (MQ=255)
                    gAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttttt                                                                        <  1:2937775/91‑1 (MQ=255)
                       gCTTTATTAACAGCTGCAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGt                                                                     <  2:3700561/91‑1 (MQ=255)
                        cTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTc                                                                    <  1:529109/91‑1 (MQ=255)
                        cTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTc                                                                    <  2:288531/91‑1 (MQ=255)
                         tttATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTActct                                                                   >  1:3121858/1‑91 (MQ=255)
                           tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                                 >  1:148390/1‑91 (MQ=255)
                           tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                                 >  1:1759232/1‑91 (MQ=255)
                            attaACAGCTGAAAAATATTCAGATATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                                <  2:3616489/91‑1 (MQ=255)
                            attaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                                <  1:87076/91‑1 (MQ=255)
                              taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCCCCGTGTAACTCCTTAATCTAAGATTCAGGAGTTTATTTTTAGTCTGAGaa                                                              >  2:3600077/1‑91 (MQ=255)
                              taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaa                                                              <  2:3862255/91‑1 (MQ=255)
                               aaCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaa                                                             >  1:3230493/1‑91 (MQ=255)
                               aaCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaa                                                             <  1:700162/91‑1 (MQ=255)
                                aCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaaa                                                            >  1:268031/1‑91 (MQ=255)
                                aCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaaa                                                            >  1:2888719/1‑91 (MQ=255)
                                 cAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaat                                                           >  1:96772/1‑91 (MQ=255)
                                 cAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaat                                                           <  1:2365333/91‑1 (MQ=255)
                                   gCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataa                                                         >  1:3202525/1‑91 (MQ=255)
                                   gCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataa                                                         >  2:330793/1‑91 (MQ=255)
                                    cTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataaa                                                        <  1:597771/91‑1 (MQ=255)
                                        aaaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGat                                                    >  1:910766/1‑91 (MQ=255)
                                          aaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGatat                                                  >  1:1978024/1‑91 (MQ=255)
                                           atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGtt                                                                               <  1:3766871/61‑1 (MQ=255)
                                           atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGtt                                                                               >  2:3766871/1‑61 (MQ=255)
                                           atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATg                                                 >  2:2721454/1‑91 (MQ=255)
                                           atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATg                                                 >  1:1636697/1‑91 (MQ=255)
                                             atTCAGAAATGATTGGTGAATAATCACCGTGTAGCCCCTTAATCTAAGCTTGAGGAGTTTACTTTTAGTCTGAGAAAATAAACGATATgaa                                               >  2:2501652/1‑91 (MQ=255)
                                             atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                               <  1:3701973/91‑1 (MQ=255)
                                             atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                               <  2:1293057/91‑1 (MQ=255)
                                             atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                               >  1:1704949/1‑91 (MQ=255)
                                              ttCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaa                                              >  2:2600530/1‑91 (MQ=255)
                                               tCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaag                                             >  1:3802358/1‑91 (MQ=255)
                                                cAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaaga                                            >  1:2082273/1‑91 (MQ=255)
                                                cAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaaga                                            >  2:3465985/1‑91 (MQ=255)
                                                  gAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaaa                                          <  2:608717/91‑1 (MQ=255)
                                                  gAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaaa                                          <  2:1934249/91‑1 (MQ=255)
                                                      tGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAatt                                      >  1:2096473/1‑91 (MQ=255)
                                                       gATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAatta                                     <  2:1324317/91‑1 (MQ=255)
                                                         ttGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattatt                                   >  1:331487/1‑91 (MQ=255)
                                                         ttGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattatt                                   >  1:2257598/1‑91 (MQ=255)
                                                         ttGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattatt                                   >  1:637342/1‑91 (MQ=255)
                                                           ggTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAat                                       >  2:1558652/1‑85 (MQ=255)
                                                            gTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAgg                                >  1:2486991/1‑91 (MQ=255)
                                                             tGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGt                               <  1:3385533/91‑1 (MQ=255)
                                                             tGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGt                               <  2:2038841/91‑1 (MQ=255)
                                                              gAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTa                              <  1:1344474/91‑1 (MQ=255)
                                                               aataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTAc                             <  2:531020/91‑1 (MQ=255)
                                                                ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                                <  2:730486/55‑1 (MQ=255)
                                                                ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                            <  2:3049969/91‑1 (MQ=255)
                                                                 tattCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACta                           >  1:2194477/4‑91 (MQ=255)
                                                                  attCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                          >  2:845399/3‑91 (MQ=255)
                                                                  aatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                          >  2:3179381/1‑91 (MQ=255)
                                                                  aatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                          >  1:3789740/1‑91 (MQ=255)
                                                                  aatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                          >  2:2129258/1‑91 (MQ=255)
                                                                  aatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                          <  2:1580795/91‑1 (MQ=255)
                                                                     cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                       >  2:2976645/1‑91 (MQ=255)
                                                                     cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                       >  1:1548115/1‑91 (MQ=255)
                                                                      aCCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTattatt                      <  2:637335/91‑1 (MQ=255)
                                                                       ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt                     <  2:2333997/91‑1 (MQ=255)
                                                                       ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt                     >  1:3653484/1‑91 (MQ=255)
                                                                       ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt                     <  1:64144/91‑1 (MQ=255)
                                                                        cGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTg                    >  1:101198/1‑91 (MQ=255)
                                                                        cGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCCGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTg                    >  2:2703597/1‑91 (MQ=255)
                                                                         gtgtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTgg                   <  2:3089305/91‑1 (MQ=255)
                                                                          tgtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGa                  >  1:1899144/1‑91 (MQ=255)
                                                                           gtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAg                 <  2:478283/91‑1 (MQ=255)
                                                                           gtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAg                 >  1:266041/1‑91 (MQ=255)
                                                                           gtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAg                 <  2:1013719/91‑1 (MQ=255)
                                                                                tCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg            >  2:3504194/1‑91 (MQ=255)
                                                                                tCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg            >  2:2532644/1‑91 (MQ=255)
                                                                                tCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTCTTAGGTACTATTATTTGGAGTATTg            >  1:3616839/1‑91 (MQ=38)
                                                                                 ccTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTGCTATTATTTGGAGTATTGc           >  1:1564305/1‑91 (MQ=255)
                                                                                       tCTTAGCTTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATGTGGAGTATTGCTACAtt     >  2:1827136/1‑91 (MQ=255)
                                                                                       tCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTACAtt     >  2:439179/1‑91 (MQ=255)
                                                                                       tCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTACAtt     <  2:1622734/91‑1 (MQ=255)
                                                                                       tCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTACAtt     >  2:1417937/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTACATTTTA  >  USA300TCH1516_ALE/890458‑890638

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: