Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,825 A→T 15.0% intergenic (+106/+154) USA300TCH1516_02485 → / ← USA300HOU_RS13450 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,608,8250AT15.0% 654.1 / 100.8 253intergenic (+106/+154)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (91/123);  new base T (23/15);  total (114/139)
Fisher's exact test for biased strand distribution p-value = 5.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

Only 100 of 265 total aligned reads displayed.

ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGG  >  USA300TCH1516_ALE/2608735‑2608915
                                                                                          |                                                                                          
atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAAtata                                                                                            <  2:326317/91‑1 (MQ=255)
 tataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATACTATAAGTATGACTAAAGCGACATCCAATATAg                                                                                           >  2:1328339/1‑91 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         <  1:2726622/91‑1 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         <  1:1391492/91‑1 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTGTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         <  2:3055195/91‑1 (MQ=255)
    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGat                                                                                        <  1:2309442/91‑2 (MQ=255)
    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                        <  2:1038697/91‑1 (MQ=255)
    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                        <  1:1933696/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGatg                                                                                       <  2:3490132/91‑3 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGatg                                                                                       <  1:1668127/91‑3 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATCTAGGACg                                                                                       >  2:3515376/1‑91 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  2:1866451/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  2:1778793/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAAGCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  2:2630435/91‑1 (MQ=255)
      cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGt                                                                                      <  1:2968151/91‑1 (MQ=255)
      cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGGTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGatgt                                                                                      <  2:2944358/91‑4 (MQ=38)
      aaaaCATGATAGAATATTATGTGAAGATCATGAATGTAATCGTTTTAAAAATAATAGGAGTGTGTGAAAAGCCACATCCAAGGTAGGACGt                                                                                      <  2:2495132/88‑1 (MQ=255)
       aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTg                                                                                     >  2:1418714/1‑91 (MQ=255)
        gACATAATAGAATTTTATATGTAAATCTTGTAGGTAACCATTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGtcg                                                                                    <  2:2941229/91‑3 (MQ=255)
          cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGATGTGGCt                                                                                  <  1:2781669/91‑1 (MQ=38)
          cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                                  <  1:2048898/91‑1 (MQ=255)
          cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                                  <  2:443489/91‑1 (MQ=255)
          cATGATAGAATTTTATATGTAAATCTTGTAGGGAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                                  >  2:3885081/1‑91 (MQ=255)
           aTGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtt                                                                                 >  2:1231130/1‑91 (MQ=255)
            tGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGATGTGGCttt                                                                                <  1:2551959/91‑1 (MQ=37)
             gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                               >  2:1647214/1‑91 (MQ=255)
              aTAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                              >  2:1050306/1‑91 (MQ=255)
                aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGt                                                                            >  2:770689/1‑91 (MQ=255)
                 gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                           >  1:467747/1‑91 (MQ=255)
                  aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTAt                                                                          <  1:795097/91‑1 (MQ=255)
                    ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCg                                                                        >  1:883735/1‑91 (MQ=255)
                     tttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                                       <  2:1521895/91‑1 (MQ=255)
                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                      >  1:2859499/1‑91 (MQ=255)
                        atatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTcc                                                                    <  1:3115578/91‑1 (MQ=255)
                          atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                  >  1:2187750/1‑91 (MQ=255)
                            gTAAATCTTGTAGGTAATCGTTTTAAAAAAAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAgct                                                                >  2:3448326/1‑89 (MQ=255)
                             tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAtt                                                               >  1:3698553/1‑91 (MQ=255)
                             tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAtt                                                               <  2:837123/91‑1 (MQ=255)
                                 tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                           <  1:3539164/91‑1 (MQ=255)
                                 tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                           >  2:3937831/1‑91 (MQ=255)
                                  cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                          <  2:514075/91‑1 (MQ=255)
                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                        >  2:2000307/1‑91 (MQ=255)
                                     gTAGGTAAGCGTTGTAAAAATCATGTAAGTATGCCTAAAGCCACCTCCTCTCTAGGACGTGGCTTTTAGTAGCGTTGCCGATCTGTggggg                                                       >  2:3359969/1‑90 (MQ=255)
                                      tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                      >  2:3724414/1‑91 (MQ=255)
                                        ggTAATCGTTTTTAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgct                                                    >  1:3048397/1‑89 (MQ=255)
                                        ggTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgct                                                    >  2:3657890/1‑89 (MQ=255)
                                          tAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAg                                                  <  1:1486737/91‑1 (MQ=255)
                                           aaTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAgg                                                 >  1:34550/1‑91 (MQ=255)
                                             tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                               <  1:3027899/91‑1 (MQ=255)
                                              cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                              <  2:1744687/91‑1 (MQ=255)
                                              cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                              >  2:1158845/1‑91 (MQ=255)
                                               gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                             <  1:1775344/91‑1 (MQ=255)
                                                 tttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTAc                                           <  2:30937/91‑1 (MQ=255)
                                                 tttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTAc                                           <  2:1540396/91‑1 (MQ=255)
                                                  ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTGAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                          <  1:3334337/91‑1 (MQ=255)
                                                   tAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                         <  1:1071118/91‑1 (MQ=255)
                                                    aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                        >  1:535331/1‑91 (MQ=255)
                                                    aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                        <  2:2286081/91‑1 (MQ=255)
                                                      aaaTAATATAAGTATGGCTAAAGCCACATCGAATATAGGACGTGGATTTTAGTATCGTTCCAGATTTTTGGGGAGTTAGGGGTTACttaat                                      <  2:279213/91‑1 (MQ=255)
                                                      aaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                      <  1:2617466/91‑1 (MQ=255)
                                                       aataatATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatt                                     >  1:122343/1‑91 (MQ=255)
                                                        ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                                    <  1:71420/91‑1 (MQ=255)
                                                        ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                                    <  1:51995/91‑1 (MQ=255)
                                                         taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                                   >  1:1805281/1‑91 (MQ=255)
                                                         taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                                   >  1:894715/1‑91 (MQ=255)
                                                         taataTAAGTATGACCAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCCGATTTGTGGGGAGTTAGGGTTTACttaattaa                                   >  1:3370414/1‑91 (MQ=255)
                                                          aataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTaaa                                  >  2:556959/1‑91 (MQ=255)
                                                             ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                               <  2:911804/91‑1 (MQ=255)
                                                             ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                               <  1:622815/91‑1 (MQ=255)
                                                             ataAGTATGACTAAAGCCACATCCAAGATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                               <  2:643741/91‑1 (MQ=255)
                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                              <  1:3094596/91‑1 (MQ=255)
                                                               aaGTATGACTAAAGCCACATCCATTATAGGACGTGGCTTTTACTATCGTTCCAGATTTGGGGGGAGATAGGGTATACGAAATAAAAGGGGt                             >  2:1377564/1‑91 (MQ=255)
                                                                   aTGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                         <  1:2348577/91‑1 (MQ=255)
                                                                     gACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGccc                       <  2:2834476/91‑1 (MQ=255)
                                                                     gACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGccc                       <  1:2081230/91‑1 (MQ=255)
                                                                      aCTAAAGCCACATCCAATCTAGGACGTGGCCTGTAGTAACGTCCCAGACCTGTGGGGCGGTAGGGTTTACGTACTTAAAGGGGACCGCCCt                      >  2:1961628/1‑91 (MQ=255)
                                                                       cTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTc                     <  1:3919711/91‑1 (MQ=255)
                                                                         aaaGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCaa                   <  1:3798492/91‑1 (MQ=255)
                                                                            gCCACATGCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                <  2:2976990/91‑1 (MQ=255)
                                                                            gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                <  1:2478328/91‑1 (MQ=255)
                                                                             ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               >  2:3828404/1‑91 (MQ=38)
                                                                             ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               >  1:1322603/1‑91 (MQ=38)
                                                                             ccACATCCAATATAGGGCGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               <  2:1600680/91‑1 (MQ=255)
                                                                             ccACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               >  1:2839783/1‑91 (MQ=255)
                                                                                 aTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc           <  2:1827807/91‑1 (MQ=255)
                                                                                  tCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAcc          >  2:2559733/1‑91 (MQ=255)
                                                                                   ccAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCa         <  1:3187287/91‑1 (MQ=255)
                                                                                    ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa        >  2:511435/3‑91 (MQ=255)
                                                                                    ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa        >  1:1210234/3‑91 (MQ=255)
                                                                                     tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTACTTAAAGGGGTCTGCCCTCAATTGCTCACCAAc       >  2:3340595/2‑91 (MQ=255)
                                                                                     tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTTACCAAc       >  2:40945/2‑91 (MQ=255)
                                                                                     tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc       >  1:1028250/2‑91 (MQ=255)
                                                                                      atatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg      >  1:2515009/1‑91 (MQ=255)
                                                                                      atatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg      >  1:366627/1‑91 (MQ=255)
                                                                                       tgtaGGCCGTGTCTTTTCGTATCGTTCCAGATTTGTGGGTACATAGGGTTTACTTAATTAAAGTGGTCTGCCCTCAATTGCTCGCCAACGa     <  2:731421/89‑1 (MQ=255)
                                                                                       tataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa     <  2:616407/91‑1 (MQ=255)
                                                                                        ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg    <  2:1676055/91‑1 (MQ=255)
                                                                                        ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg    <  1:937564/91‑1 (MQ=255)
                                                                                         taGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAgg   <  1:3891085/91‑1 (MQ=255)
                                                                                         taGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAgg   <  2:1266935/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGG  >  USA300TCH1516_ALE/2608735‑2608915

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