Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,829 C→T 11.1% intergenic (+110/+150) 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,8290CT11.1% 726.5 / 57.7 251intergenic (+110/+150)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (94/128);  new base T (22/6);  total (117/134)
Fisher's exact test for biased strand distribution p-value = 4.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 260 total aligned reads displayed.

AGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAG  >  USA300TCH1516_ALE/2608739‑2608919
                                                                                          |                                                                                          
aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTACAGCCACATCCAATATAGGAc                                                                                            >  1:91723/1‑91 (MQ=255)
aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                            >  1:3055309/1‑91 (MQ=255)
aGCAGACATGATAGAATTCTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATTACTAAAGTCACATCCAATATTGGAc                                                                                            <  2:1038713/91‑1 (MQ=255)
   aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTg                                                                                         >  2:1792465/1‑91 (MQ=255)
   aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTg                                                                                         >  2:1418714/1‑91 (MQ=255)
    gACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                                        <  1:1765639/91‑1 (MQ=255)
      cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGATGTGGCt                                                                                      <  1:2781549/91‑1 (MQ=38)
      cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                                      <  2:3042872/91‑1 (MQ=255)
      cATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCt                                                                                      <  1:2048898/91‑1 (MQ=255)
       atgatgGAATATGATATGTAGATCATGTAGGTAATCGTGTTGAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGATGTGGCtt                                                                                     <  2:1169753/91‑1 (MQ=19)
        tGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGATGTGGCttt                                                                                    <  1:2551959/91‑1 (MQ=37)
        tGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCttt                                                                                    <  2:45100/91‑1 (MQ=255)
        tGATAGAAGTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAGTATAAGTATGACTGAAGCCACATCCAATATAGGACGTGGCttt                                                                                    <  2:216217/91‑1 (MQ=255)
           tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                                 <  2:332518/91‑1 (MQ=255)
            aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGt                                                                                >  2:770613/1‑91 (MQ=255)
             gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                               >  2:8469/1‑91 (MQ=255)
              aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCCATATAGGACGTGGCTTTTAGTAt                                                                              <  2:1288947/91‑1 (MQ=255)
              aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTAt                                                                              <  1:490373/91‑1 (MQ=255)
                 tttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                                           <  2:3839868/91‑1 (MQ=255)
                  ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                          >  1:2859473/1‑91 (MQ=255)
                  ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                          >  1:2860024/1‑91 (MQ=255)
                   tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                         <  1:2481751/91‑1 (MQ=255)
                   tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                         <  2:43500/91‑1 (MQ=255)
                    atatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTcc                                                                        <  1:3115578/91‑1 (MQ=255)
                     tatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCa                                                                       >  2:2541935/1‑91 (MQ=255)
                      atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                      >  1:2187750/1‑91 (MQ=255)
                        gTAAATCTTGTAGGTAATCGTTTTAAAAAAAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAgct                                                                    >  2:3448326/1‑89 (MQ=255)
                         tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAtt                                                                   <  2:837123/91‑1 (MQ=255)
                         tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAtt                                                                   >  1:773776/1‑91 (MQ=255)
                            aTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgt                                                                >  2:2845311/1‑91 (MQ=255)
                            aTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgt                                                                <  1:3894824/91‑1 (MQ=255)
                             tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                               <  1:1847177/91‑1 (MQ=255)
                             tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                               <  1:3539164/91‑1 (MQ=255)
                              cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                              <  2:513948/91‑1 (MQ=255)
                               ttGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTggg                                                             >  2:1524064/1‑91 (MQ=255)
                                tGTAGGTAATCGTTTTAAAAATTATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                            >  2:1007747/1‑91 (MQ=255)
                                tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                            <  2:269830/91‑1 (MQ=255)
                                cgTAGGTAATCGTTTTAAAAATACTATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCCGCCTTGTgggg                                                            >  2:3599655/2‑91 (MQ=255)
                                 gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGGCGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                           <  2:395458/91‑1 (MQ=255)
                                 gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                           >  1:2517812/1‑91 (MQ=255)
                                 gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                           >  2:1775257/1‑91 (MQ=255)
                                   aGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                         <  2:1884318/91‑1 (MQ=255)
                                    ggTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgct                                                        >  2:3657890/1‑89 (MQ=255)
                                    ggTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                        >  2:1993992/1‑91 (MQ=255)
                                     gTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                                       <  1:3840883/91‑1 (MQ=255)
                                          cGTTTTAAAAATAGTATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                                  <  2:165970/91‑1 (MQ=255)
                                          cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                                  <  2:1744598/91‑1 (MQ=255)
                                           gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                                 <  1:1775707/91‑1 (MQ=255)
                                           gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                                 >  2:2943224/1‑91 (MQ=255)
                                              ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                              >  1:2447485/1‑91 (MQ=255)
                                               tAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                             <  1:1071118/91‑1 (MQ=255)
                                               tAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                             <  2:548734/91‑1 (MQ=255)
                                                aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                            <  2:32589/91‑1 (MQ=255)
                                                aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtta                                            >  1:535199/1‑91 (MQ=255)
                                                 aaaaTAATATAAGTATGACTAAAGCCACATCCAAAATTGGTAGTGGTTTTTAGTATCGTTTCAGTTTTGTGGATAGTTAGGGTTTACttga                                           >  2:603635/1‑89 (MQ=255)
                                                  aaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                          <  1:2617466/91‑1 (MQ=255)
                                                   aataatATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatt                                         >  1:122343/1‑91 (MQ=255)
                                                    ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                                        <  1:71420/91‑1 (MQ=255)
                                                    ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                                        <  2:2597127/91‑1 (MQ=255)
                                                     taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTGTCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                                       >  1:893966/1‑91 (MQ=255)
                                                     taataTAAGTATGACCAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCCGATTTGTGGGGAGTTAGGGTTTACttaattaa                                       >  1:3370414/1‑91 (MQ=255)
                                                        tataATTATGACTAAAGCCACATCCAATGTAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgg                                    <  2:280653/91‑1 (MQ=255)
                                                         ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                                   >  1:2094826/1‑91 (MQ=255)
                                                         ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                                   <  1:622815/91‑1 (MQ=255)
                                                           aaGTATGACTAAAGCCACATCCATTATAGGACGTGGCTTTTACTATCGTTCCAGATTTGGGGGGAGATAGGGTATACGAAATAAAAGGGGt                                 >  2:1377564/1‑91 (MQ=255)
                                                           aaGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                                 <  2:1459361/91‑1 (MQ=255)
                                                           aaGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTGAGGGTTTACTTAATTAAAggggc                                 >  2:639959/1‑90 (MQ=255)
                                                           aaGTATGACTAAAGCCACATACAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                                 >  1:2340286/1‑91 (MQ=255)
                                                            aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTc                                >  1:310991/1‑91 (MQ=255)
                                                             gTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                               <  2:3349127/91‑1 (MQ=255)
                                                               aTGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                             <  1:3472584/91‑1 (MQ=255)
                                                                  aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                          >  2:154691/1‑91 (MQ=255)
                                                                      aaGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAAt                      >  1:3930021/1‑91 (MQ=255)
                                                                       aGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAAtc                     >  2:3923159/1‑90 (MQ=255)
                                                                        gCCACATGCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                    <  2:2976990/91‑1 (MQ=255)
                                                                        gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                    >  1:3038043/1‑91 (MQ=255)
                                                                         ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc                   >  2:3828404/1‑91 (MQ=38)
                                                                         ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc                   >  1:1322446/1‑91 (MQ=38)
                                                                         ccACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATCTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc                   >  2:965737/1‑91 (MQ=255)
                                                                          cacaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCt                  <  1:2985605/91‑1 (MQ=255)
                                                                            caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                >  2:228863/1‑91 (MQ=255)
                                                                            caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATGTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                <  2:3710184/91‑1 (MQ=255)
                                                                               ccAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCa             >  1:3756753/1‑91 (MQ=255)
                                                                                ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa            >  1:1210234/3‑91 (MQ=255)
                                                                                ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa            >  1:148470/3‑91 (MQ=255)
                                                                                ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa            >  1:3205202/3‑91 (MQ=255)
                                                                                 tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc           >  1:1028250/2‑91 (MQ=255)
                                                                                 tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc           >  1:2822522/2‑91 (MQ=255)
                                                                                 tatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGCGGTCTGCCCTCAATTGCTCACCAAc           >  1:3687727/2‑91 (MQ=255)
                                                                                 aaTATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc           <  2:1864845/91‑1 (MQ=255)
                                                                                  atatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg          >  1:3365143/1‑91 (MQ=255)
                                                                                  atatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg          >  2:1641971/1‑91 (MQ=255)
                                                                                  atatAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg          >  1:878272/1‑91 (MQ=255)
                                                                                   tgtaGGCCGTGTCTTTTCGTATCGTTCCAGATTTGTGGGTACATAGGGTTTACTTAATTAAAGTGGTCTGCCCTCAATTGCTCGCCAACGa         <  2:731421/89‑1 (MQ=255)
                                                                                    ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTTACCAACGAg        <  1:1970600/91‑1 (MQ=255)
                                                                                    ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg        >  1:548545/1‑91 (MQ=255)
                                                                                    ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg        >  2:1150360/1‑91 (MQ=255)
                                                                                     taGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAgg       <  2:1266935/91‑1 (MQ=255)
                                                                                       ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGc     <  1:1522368/91‑1 (MQ=255)
                                                                                       ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGc     <  2:3506433/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAG  >  USA300TCH1516_ALE/2608739‑2608919

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: