Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 890,542 T→C 9.5% intergenic (+17/‑47) 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,5420TC 9.5% 381.6 / 24.3 126intergenic (+17/‑47)gcvH/USA300TCH1516_00820Glycine cleavage system H protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (60/54);  new base C (5/7);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 129 total aligned reads displayed.

AGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTA  >  USA300TCH1516_ALE/890453‑890631
                                                                                         |                                                                                         
aGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTa                                                                                          >  1:1369231/1‑91 (MQ=255)
aGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTa                                                                                          >  1:145138/1‑91 (MQ=255)
 gAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCCTaa                                                                                         >  1:3346815/1‑91 (MQ=255)
  aaaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAAt                                                                                        >  1:2030468/1‑91 (MQ=255)
  aaaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAAt                                                                                        >  1:3570191/1‑91 (MQ=255)
   aaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATc                                                                                       >  1:474597/1‑91 (MQ=255)
       aGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAg                                                                                   >  1:1699221/1‑91 (MQ=255)
       aGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCttag                                                                                   <  1:1510226/91‑4 (MQ=255)
        gTGATGAAAGTCAGATTGAAGCTTTATTATCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGa                                                                                  <  1:3334776/91‑1 (MQ=255)
        gTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTCAgc                                                                                  >  1:3953627/1‑90 (MQ=255)
        gTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGa                                                                                  >  1:3879548/1‑91 (MQ=255)
        gTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGa                                                                                  <  1:2908230/91‑1 (MQ=255)
          gatgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGAtt                                                                                <  1:2597690/91‑1 (MQ=255)
           atgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTg                                                                               >  1:2756134/1‑91 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGATAAATATTCAGAAATGATTGGTGAATCATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                              <  1:3683610/91‑3 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                              <  1:3890453/91‑3 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                              <  1:1970264/91‑3 (MQ=255)
             gaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgag                                                                             >  1:3110310/1‑91 (MQ=255)
             gaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaag                                                                             <  1:1572200/91‑4 (MQ=255)
               aaGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgagga                                                                           <  1:1835163/91‑1 (MQ=255)
                    aGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTa                                                                      >  1:1958455/1‑91 (MQ=255)
                     gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                     <  1:500085/91‑1 (MQ=255)
                     gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                     <  1:500162/91‑1 (MQ=255)
                      aTTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAtt                                                                    >  1:2938127/1‑91 (MQ=255)
                        tGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAtttt                                                                  <  1:1159148/91‑1 (MQ=255)
                        tGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAtttt                                                                  <  1:3109405/91‑1 (MQ=255)
                         gAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttttt                                                                 <  1:3485233/91‑1 (MQ=255)
                         gAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttttt                                                                 >  1:578223/1‑91 (MQ=255)
                            gCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGt                                                              >  1:1956486/1‑91 (MQ=255)
                            gCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGt                                                              <  1:1980689/91‑1 (MQ=255)
                             cTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTc                                                             <  1:1568079/91‑1 (MQ=255)
                              tttATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCt                                                            >  1:3323694/1‑91 (MQ=255)
                               ttattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTg                                                           <  1:3064645/91‑1 (MQ=255)
                               ttattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTg                                                           <  1:3064626/91‑1 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                          >  1:3844104/1‑91 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                          <  1:4080441/91‑1 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                          >  1:2627614/1‑91 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAACATTGAGGAGTTTATTTTTAGTCTga                                                          >  1:3844171/1‑91 (MQ=255)
                                 attaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                         <  1:1199115/91‑1 (MQ=255)
                                  ttaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgaga                                                        >  1:1015489/1‑91 (MQ=255)
                                   taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaa                                                       >  1:1501756/1‑91 (MQ=255)
                                   taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaa                                                       <  1:2791670/91‑1 (MQ=255)
                                     aCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaaa                                                     >  1:4007626/1‑91 (MQ=255)
                                      cAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaat                                                    >  1:3636436/1‑91 (MQ=255)
                                       aGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaata                                                   >  1:2531742/1‑91 (MQ=255)
                                       aGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaata                                                   <  1:3508255/91‑1 (MQ=255)
                                           gAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATg                                               >  1:850673/1‑91 (MQ=255)
                                           gAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATg                                               >  1:4072291/1‑91 (MQ=255)
                                           gAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATg                                               >  1:2889889/1‑91 (MQ=255)
                                            aaaaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGa                                              <  1:2266900/91‑1 (MQ=255)
                                              aaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGata                                            <  1:2466306/91‑1 (MQ=255)
                                              aaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGata                                            >  1:3745794/1‑91 (MQ=255)
                                              aaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGata                                            <  1:2061032/91‑1 (MQ=255)
                                                  atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                        <  1:185404/91‑1 (MQ=255)
                                                  atTCAGAAATGATTGGTGAATAATCACCGAGGAACTCCTTAATCGAAGATTGAGGAGGTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                        <  1:3316667/91‑1 (MQ=255)
                                                   ttCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaa                                       >  1:2432314/1‑91 (MQ=255)
                                                   ttCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGACTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaa                                       >  1:1617564/1‑91 (MQ=255)
                                                    tCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaag                                      <  1:1869676/91‑1 (MQ=255)
                                                     cAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaaga                                     <  1:667177/91‑1 (MQ=255)
                                                     cAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaaga                                     >  1:2519759/1‑91 (MQ=255)
                                                      aGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaa                                    <  1:329748/91‑1 (MQ=255)
                                                       gAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaaa                                   <  1:2088804/91‑1 (MQ=255)
                                                            gATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAatta                              >  1:117583/1‑91 (MQ=255)
                                                             aTTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattat                             >  1:1651948/1‑91 (MQ=255)
                                                             aTTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGACTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattat                             >  1:2208370/1‑91 (MQ=255)
                                                                ggTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAg                          <  1:3916644/91‑1 (MQ=255)
                                                                  tGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGt                        <  1:3133492/91‑1 (MQ=255)
                                                                   gAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTa                       <  1:446816/91‑1 (MQ=255)
                                                                    aataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTAc                      <  1:3898140/91‑1 (MQ=255)
                                                                    aataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTAc                      <  1:749223/91‑1 (MQ=255)
                                                                    aataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTAc                      >  1:338123/1‑91 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                     <  1:2704775/91‑1 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                     <  1:4048560/91‑1 (MQ=255)
                                                                      tattCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACta                    >  1:165227/4‑91 (MQ=255)
                                                                      taatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACta                    <  1:1278493/91‑1 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                >  1:3083170/1‑91 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                >  1:1125294/1‑91 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                >  1:173739/1‑91 (MQ=255)
                                                                           aCCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTattatt               >  1:2420281/1‑91 (MQ=255)
                                                                           aCCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTattatt               >  1:1012662/1‑91 (MQ=255)
                                                                           aCCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTattatt               <  1:105976/91‑1 (MQ=255)
                                                                            ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt              <  1:3964700/91‑1 (MQ=255)
                                                                              aataaaCTCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTgg            >  1:1285461/5‑91 (MQ=37)
                                                                               tgtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGa           >  1:453789/1‑91 (MQ=255)
                                                                               tgtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGa           >  1:1423337/1‑91 (MQ=255)
                                                                                gtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAg          >  1:3820962/1‑91 (MQ=255)
                                                                                 tAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGt         <  1:1515418/91‑1 (MQ=255)
                                                                                   aCTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt       <  1:1350976/91‑1 (MQ=255)
                                                                                   aCTCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt       >  1:3896269/1‑91 (MQ=39)
                                                                                   aCTCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt       >  1:3156010/1‑91 (MQ=38)
                                                                                   aCTCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt       >  1:5169/1‑91 (MQ=38)
                                                                                    ctcctcAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAtt      >  1:95837/1‑91 (MQ=255)
                                                                                    cTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAtt      >  1:1718003/1‑91 (MQ=255)
                                                                                    cTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAtt      >  1:1718189/1‑91 (MQ=255)
                                                                                     tCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg     >  1:2367573/1‑91 (MQ=255)
                                                                                     tCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg     >  1:3086134/1‑91 (MQ=255)
                                                                                     tCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg     >  1:959785/1‑91 (MQ=39)
                                                                                      ccTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGc    >  1:2308334/1‑91 (MQ=255)
                                                                                      ccTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGc    >  1:3070366/1‑91 (MQ=255)
                                                                                       cTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCt   <  1:674986/91‑1 (MQ=255)
                                                                                         |                                                                                         
AGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTA  >  USA300TCH1516_ALE/890453‑890631

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