Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,681 T→A 12.2% intergenic (+27/‑302) USA300HOU_RS04805 → / → USA300HOU_RS04810 Ureidoglycolate lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE970,6810TA12.2% 709.1 / 74.0 237intergenic (+27/‑302)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (134/74);  new base A (13/16);  total (147/90)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

Only 100 of 245 total aligned reads displayed.

AGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAA  >  USA300TCH1516_ALE/970591‑970771
                                                                                          |                                                                                          
 ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                           <  1:1251629/91‑1 (MQ=255)
  gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                          <  2:1131850/91‑1 (MQ=255)
     atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       >  1:1535969/1‑91 (MQ=255)
     atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       >  2:2571490/1‑91 (MQ=255)
      tgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatata                                                                                      <  2:3291213/91‑1 (MQ=255)
       gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                     >  2:884170/1‑91 (MQ=255)
        aaGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGa                                                                                    >  2:1469156/1‑91 (MQ=255)
        aaGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGa                                                                                    >  1:3808111/1‑91 (MQ=255)
        aaGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGa                                                                                    <  2:626552/91‑1 (MQ=255)
         aGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAt                                                                                   >  2:3143837/1‑91 (MQ=255)
           ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTg                                                                                 <  1:3281342/91‑1 (MQ=255)
           ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGCCTTTATATAGATTg                                                                                 >  2:1934093/1‑91 (MQ=255)
           ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGactg                                                                                 <  1:2743547/91‑4 (MQ=38)
            tAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTgg                                                                                <  2:1874109/91‑1 (MQ=255)
             aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                               >  1:2123411/1‑91 (MQ=255)
                gTGACTATTGATAATATTGGAACGCTGACGACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                            <  1:3108633/91‑1 (MQ=37)
                gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            >  1:1444166/1‑91 (MQ=255)
                gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            >  1:3025836/1‑91 (MQ=255)
                 tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                           >  2:3784894/1‑91 (MQ=255)
                   aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                         >  1:652692/1‑91 (MQ=255)
                   aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                         >  1:386657/1‑91 (MQ=255)
                   aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                         >  2:106265/1‑91 (MQ=255)
                   aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                         >  1:3733319/1‑91 (MQ=255)
                   aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAgctt                                                                         <  2:243057/91‑4 (MQ=37)
                      aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                      >  2:77612/1‑91 (MQ=255)
                      aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                      <  1:985428/91‑1 (MQ=21)
                       ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGCTTTTTc                                                                     >  1:2184096/1‑91 (MQ=255)
                       ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                     <  1:1179463/91‑2 (MQ=21)
                       ttGATAATATTGGAACGATGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                     <  2:3721721/91‑2 (MQ=18)
                        tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                    <  1:202060/91‑1 (MQ=255)
                        tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                    <  2:3005202/91‑1 (MQ=255)
                         gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                   >  1:297002/1‑91 (MQ=255)
                          ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGctttttaaat                                                                  <  2:802440/91‑5 (MQ=18)
                              tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                              >  1:280788/1‑91 (MQ=255)
                              tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                              <  1:2295049/91‑1 (MQ=255)
                               atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtttt                                                             <  1:3663259/91‑1 (MQ=255)
                                ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCCTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                            >  2:1455753/1‑91 (MQ=255)
                                 tGGAACGGTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                           <  1:1941616/91‑1 (MQ=255)
                                 tGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                           <  2:2189743/91‑1 (MQ=255)
                                  ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                          >  2:3642451/1‑91 (MQ=255)
                                   gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                         <  1:3362658/91‑1 (MQ=255)
                                   gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                         >  2:703231/1‑91 (MQ=255)
                                   gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                         >  1:1486410/1‑91 (MQ=255)
                                   gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                         <  2:1816776/91‑1 (MQ=255)
                                   gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGGTAGTTTTTTCTTGCTTTTCTac                                                         >  1:3554187/1‑90 (MQ=255)
                                    aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                        >  1:3859748/1‑91 (MQ=255)
                                    aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                        >  2:520250/1‑91 (MQ=255)
                                    aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                        <  1:2480715/91‑1 (MQ=255)
                                      cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                      <  2:2256602/91‑1 (MQ=255)
                                      cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                      >  2:2933632/1‑91 (MQ=255)
                                       gCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAg                                                     >  2:2631995/1‑91 (MQ=255)
                                        cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                    >  1:178537/1‑91 (MQ=255)
                                        cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                    >  2:31727/1‑91 (MQ=255)
                                        cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                    >  1:109677/1‑91 (MQ=255)
                                         tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                   >  2:1372759/1‑91 (MQ=255)
                                            cAACTTATATTGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                >  2:29012/1‑91 (MQ=255)
                                            cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                >  1:918605/1‑91 (MQ=255)
                                             aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                               >  1:3853669/1‑91 (MQ=255)
                                             aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                               <  2:3350216/91‑1 (MQ=255)
                                             aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                               <  1:57898/91‑1 (MQ=255)
                                                 tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                           <  2:2834698/91‑1 (MQ=255)
                                                 tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                           <  2:85696/91‑1 (MQ=255)
                                                 tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                           <  1:881176/91‑1 (MQ=255)
                                                 tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                           <  2:3193112/91‑1 (MQ=255)
                                                   tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATa                                         <  1:1481809/91‑1 (MQ=255)
                                                    atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                        <  1:1053105/91‑1 (MQ=255)
                                                    atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                        >  2:2729238/1‑91 (MQ=255)
                                                    atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                        >  1:756470/1‑91 (MQ=255)
                                                     tCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaaa                                       <  1:1608896/91‑1 (MQ=255)
                                                      cGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAat                                      >  1:2692798/1‑91 (MQ=255)
                                                        cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattc                                    >  2:419723/1‑90 (MQ=255)
                                                        cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                    >  1:3440585/1‑91 (MQ=255)
                                                         tAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattat                                   >  1:1334919/1‑91 (MQ=255)
                                                           aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                                 >  2:1326949/1‑91 (MQ=255)
                                                               atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                             >  1:2534444/1‑91 (MQ=255)
                                                                 tATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATg                           >  1:1681785/1‑91 (MQ=255)
                                                                  aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                          <  2:1647757/91‑1 (MQ=255)
                                                                   tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                         <  1:446152/91‑1 (MQ=255)
                                                                    cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                        >  2:1670573/1‑91 (MQ=255)
                                                                    cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                        >  2:3048295/1‑91 (MQ=255)
                                                                     aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                       <  2:2084984/91‑1 (MQ=255)
                                                                       ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                     >  1:3761983/1‑91 (MQ=255)
                                                                        tAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAt                    >  1:1004662/1‑91 (MQ=255)
                                                                           aaaacTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 >  2:2799338/5‑91 (MQ=25)
                                                                           aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 <  1:454661/91‑1 (MQ=255)
                                                                            aaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTg                >  1:3032254/1‑91 (MQ=255)
                                                                              gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag              >  1:2214328/1‑91 (MQ=255)
                                                                              acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag              >  2:2271340/2‑91 (MQ=37)
                                                                              acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag              >  1:1550334/2‑91 (MQ=37)
                                                                               cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga             >  1:2578692/1‑91 (MQ=37)
                                                                               cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga             >  2:1764338/1‑91 (MQ=37)
                                                                                     aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg       >  1:1889218/1‑91 (MQ=255)
                                                                                     aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg       <  2:226718/91‑1 (MQ=255)
                                                                                      gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt      <  1:764672/91‑1 (MQ=255)
                                                                                       tCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg     >  2:1424943/1‑91 (MQ=255)
                                                                                       tCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg     >  2:1424650/1‑91 (MQ=255)
                                                                                        cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa    >  1:1362674/1‑91 (MQ=255)
                                                                                          ttATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGaaa  >  2:2741502/1‑91 (MQ=255)
                                                                                          ttATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGaaa  >  2:2407907/1‑91 (MQ=255)
                                                                                          ttATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGaaa  >  2:1566725/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAA  >  USA300TCH1516_ALE/970591‑970771

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: