Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 890,542 T→C 10.7% 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,5420TC10.7% 716.3 / 61.5 253intergenic (+17/‑47)gcvH/USA300TCH1516_00820Glycine cleavage system H protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (121/105);  new base C (8/19);  total (129/124)
Fisher's exact test for biased strand distribution p-value = 2.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 256 total aligned reads displayed.

AGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTAC  >  USA300TCH1516_ALE/890453‑890632
                                                                                         |                                                                                          
aGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTa                                                                                           >  1:288531/1‑91 (MQ=255)
aGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTa                                                                                           >  1:2236347/1‑91 (MQ=255)
 gAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTaa                                                                                          <  1:764740/91‑1 (MQ=255)
 gAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTaa                                                                                          <  1:558258/91‑1 (MQ=255)
  aaaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAAt                                                                                         <  2:913923/91‑1 (MQ=255)
  aaaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAAt                                                                                         >  2:1937807/1‑91 (MQ=255)
   aaTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATc                                                                                        >  1:769707/1‑91 (MQ=255)
    aTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCt                                                                                       <  1:1960100/91‑1 (MQ=255)
    aTTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCt                                                                                       <  2:1732775/91‑1 (MQ=255)
         tgatgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAt                                                                                  <  1:1063012/91‑1 (MQ=255)
          gatgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGAtt                                                                                 >  2:1599650/1‑91 (MQ=255)
          gatgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAtt                                                                                 <  1:2669094/91‑1 (MQ=255)
           atgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTg                                                                                >  1:3795585/1‑91 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                               <  2:2642117/91‑3 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                               <  1:2838299/91‑3 (MQ=255)
            tgaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGAttaa                                                                               <  1:1515188/91‑3 (MQ=255)
             gaAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgag                                                                              <  2:2469574/91‑1 (MQ=255)
              aaaGTCAGCTTGAAGCATTATTATCAGCTGAAAAATACGCAGAAATGATTGGTGAATAATCACCGTATGACTCCTCAATCTTAGATTaagg                                                                             <  2:1216604/91‑1 (MQ=255)
              aaaGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgagg                                                                             <  2:452866/91‑1 (MQ=255)
              aaaGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTCAATCTTAGATTaagg                                                                             <  2:2023949/91‑1 (MQ=255)
                aGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTgaggag                                                                           <  1:98658/91‑1 (MQ=255)
                  tCAGATTGAAGCTTTATTAACAGCTGAAAAATCTTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGtt                                                                         >  2:3811904/1‑91 (MQ=255)
                   cAGCTTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGttt                                                                        <  1:3497128/91‑1 (MQ=255)
                    aGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTa                                                                       >  1:389809/1‑91 (MQ=255)
                     gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                      <  2:479105/91‑1 (MQ=255)
                     gATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGAGTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAt                                                                      <  2:477987/91‑1 (MQ=255)
                       ttGAAGCTNTATNANCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttt                                                                    <  1:899474/91‑1 (MQ=255)
                         gAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTAttttt                                                                  <  1:2936860/91‑1 (MQ=255)
                          aaGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTa                                                                 <  2:404166/91‑1 (MQ=255)
                             cTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTc                                                              <  2:288531/91‑1 (MQ=255)
                             cTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTc                                                              <  1:2441628/91‑1 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                           <  2:636837/91‑1 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                           <  2:213117/91‑1 (MQ=255)
                                tattaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTga                                                           >  1:1759232/1‑91 (MQ=255)
                                 attaACAGCTGAAAAATATTCAGATATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                          <  2:3616489/91‑1 (MQ=255)
                                 attaACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                          <  1:1528359/91‑1 (MQ=255)
                                   taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaa                                                        <  2:3862255/91‑1 (MQ=255)
                                   taACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaa                                                        >  1:3876257/1‑91 (MQ=255)
                                    aaCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTGTGAGaaa                                                       >  2:3456848/1‑91 (MQ=255)
                                    aaCAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGaaa                                                       <  1:3236074/91‑1 (MQ=255)
                                      cAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaat                                                     >  1:2912266/1‑91 (MQ=255)
                                       aGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaata                                                    >  2:1565463/1‑91 (MQ=255)
                                        gCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataa                                                   >  1:3202525/1‑91 (MQ=255)
                                         cTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataaa                                                  <  1:597771/91‑1 (MQ=255)
                                          tGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAaaataaat                                                 <  1:1599650/91‑1 (MQ=255)
                                           gAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATg                                                >  1:2352878/1‑91 (MQ=255)
                                             aaaaTATTCAGAAATGATTGGTGACTAATCACCGTGTAACTCCTCAATCTAAGATTGAGGAGTTTATTTTTAGTCTGTGTAAATCAATGat                                              >  2:2303639/1‑91 (MQ=255)
                                             aaaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGat                                              >  1:3160362/1‑91 (MQ=255)
                                             aaaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGat                                              >  1:3625119/1‑91 (MQ=255)
                                              aaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGata                                             >  2:3012469/1‑91 (MQ=255)
                                              aaaTATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGata                                             >  2:1803508/1‑91 (MQ=255)
                                                atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGtt                                                                         <  1:3766871/61‑1 (MQ=255)
                                                atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGtt                                                                         >  2:3766871/1‑61 (MQ=255)
                                                atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATg                                           >  1:216327/1‑91 (MQ=255)
                                                atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATg                                           >  2:2846891/1‑91 (MQ=255)
                                                atatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATg                                           >  1:1636697/1‑91 (MQ=255)
                                                 tatTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATga                                          >  1:640500/1‑91 (MQ=255)
                                                  atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                         >  1:1704949/1‑91 (MQ=255)
                                                  atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                         <  2:1293057/91‑1 (MQ=255)
                                                  atTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaa                                         <  1:3701973/91‑1 (MQ=255)
                                                    tCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaag                                       <  2:2431369/91‑1 (MQ=255)
                                                    tCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaag                                       <  1:275943/91‑1 (MQ=255)
                                                     cAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaaga                                      >  1:851298/1‑91 (MQ=255)
                                                      aGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaa                                     <  2:3173995/91‑1 (MQ=255)
                                                       gAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaaa                                    <  2:608717/91‑1 (MQ=255)
                                                       gAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATgaaagaaa                                    <  1:2588852/91‑1 (MQ=255)
                                                            gATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAatta                               <  2:1324317/91‑1 (MQ=255)
                                                              ttGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAattatt                             >  1:331487/1‑91 (MQ=255)
                                                                ggTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAg                           >  2:1017993/1‑91 (MQ=255)
                                                                 gTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAgg                          >  1:2487004/1‑91 (MQ=255)
                                                                  tGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGt                         <  2:2038841/91‑1 (MQ=255)
                                                                  tGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGt                         <  1:3385474/91‑1 (MQ=255)
                                                                   gAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTa                        <  1:1344474/91‑1 (MQ=255)
                                                                    aataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTAc                       <  2:2598460/91‑1 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                          >  1:730486/1‑55 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTgag                                                          <  2:730486/55‑1 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                      <  2:250274/91‑1 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                      <  2:3049969/91‑1 (MQ=255)
                                                                     ataatCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                      <  2:2225639/91‑1 (MQ=255)
                                                                       attCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtat                    >  2:845399/3‑91 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                 >  1:2684090/1‑91 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACtattat                 >  1:956679/1‑91 (MQ=255)
                                                                          cACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACt                      <  2:3222949/86‑1 (MQ=255)
                                                                           aCCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTattatt                >  1:945596/1‑91 (MQ=255)
                                                                            ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt               <  1:64186/91‑1 (MQ=255)
                                                                            ccGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTAttt               >  1:3436047/1‑91 (MQ=255)
                                                                             cGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTg              >  1:101198/1‑91 (MQ=255)
                                                                             cGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCCGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTg              >  2:2703597/1‑91 (MQ=255)
                                                                               tgtAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGa            >  1:1899144/1‑91 (MQ=255)
                                                                                taaaCTCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAg           >  1:1253524/3‑91 (MQ=38)
                                                                                   aCTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt        >  2:3664025/1‑91 (MQ=255)
                                                                                   aCTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt        >  1:2083896/1‑91 (MQ=255)
                                                                                   aCTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAt        >  1:2023135/1‑91 (MQ=255)
                                                                                    cTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAtt       >  2:359230/1‑91 (MQ=255)
                                                                                    cTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTAtt       <  1:2699655/91‑1 (MQ=255)
                                                                                     tCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTg      >  2:3504194/1‑91 (MQ=255)
                                                                                     tCCTCAATCTTAGATTAAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTCTTAGGTACTATTATTTGGAGTATTg      >  1:3616839/1‑91 (MQ=38)
                                                                                      ccTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTGCTATTATTTGGAGTATTGc     >  1:1564305/1‑91 (MQ=255)
                                                                                      ccTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGc     <  1:3032723/91‑1 (MQ=255)
                                                                                         tAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTAGAGTATTGCTAc  <  1:2976645/91‑1 (MQ=255)
                                                                                         |                                                                                          
AGAAATTAGTGATGAAAGTCAGATTGAAGCTTTATTAACAGCTGAAAAATATTCAGAAATGATTGGTGAATAATCACCGTGTAACTCCTTAATCTAAGATTGAGGAGTTTATTTTTAGTCTGAGAAAATAAATGATATGAAAGAAAAATTATTAGGTACTATTATTTGGAGTATTGCTAC  >  USA300TCH1516_ALE/890453‑890632

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