Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,677 G→A 12.4% intergenic (+23/‑306) 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,6770GA12.4% 681.7 / 60.8 234intergenic (+23/‑306)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (132/73);  new base A (13/16);  total (145/89)
Fisher's exact test for biased strand distribution p-value = 6.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 239 total aligned reads displayed.

AACCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGT  >  USA300TCH1516_ALE/970587‑970767
                                                                                          |                                                                                          
     ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                       <  1:1251629/91‑1 (MQ=255)
      gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                      <  2:1131850/91‑1 (MQ=255)
       tgatgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtat                                                                                     >  2:1801400/1‑91 (MQ=255)
         atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                   >  1:1534853/1‑91 (MQ=255)
          tgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatata                                                                                  <  2:3291311/91‑1 (MQ=255)
           gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                 >  2:884170/1‑91 (MQ=255)
           gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                 >  1:1294844/1‑91 (MQ=255)
           gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                 >  1:2834698/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)
              gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                              >  1:879459/1‑91 (MQ=255)
               ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGactg                                                                             <  1:2743547/91‑4 (MQ=38)
                tAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTgg                                                                            <  2:3757641/91‑1 (MQ=37)
                 aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                           >  1:2123411/1‑91 (MQ=255)
                  aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                          <  2:221675/91‑1 (MQ=255)
                    gTGACTATTGATAATATTGGAACGCTGACGACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                        <  1:3108633/91‑1 (MQ=37)
                    gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                        >  2:1666603/1‑91 (MQ=255)
                      gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtt                                                                      >  1:1416446/1‑91 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     >  1:3733319/1‑91 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     >  1:738761/1‑91 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     >  1:652692/1‑91 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     >  1:386657/1‑91 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAgctt                                                                     <  2:243057/91‑4 (MQ=37)
                        cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtatt                                                                    >  2:2820920/1‑91 (MQ=255)
                         tATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttttt                                                                   >  2:3045249/1‑91 (MQ=255)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                  >  2:77612/1‑91 (MQ=255)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                  >  2:1759052/1‑91 (MQ=255)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                  <  1:985428/91‑1 (MQ=21)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                  <  1:895093/91‑1 (MQ=21)
                           ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                 >  2:1548879/1‑91 (MQ=255)
                           ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGCTTTTTc                                                                 >  1:2184096/1‑91 (MQ=255)
                           ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                 <  1:1179463/91‑2 (MQ=21)
                            tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                <  1:202060/91‑1 (MQ=255)
                            tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGCTTAGTTTTTTCt                                                                <  1:2990491/91‑1 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                              >  2:1412858/1‑91 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGctttttaaat                                                              <  2:802440/91‑5 (MQ=18)
                                  tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                          <  1:2295049/91‑1 (MQ=255)
                                     tGGAACGGTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                       <  1:1941616/91‑1 (MQ=255)
                                      ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                      <  2:1447830/91‑1 (MQ=255)
                                       gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                     >  2:1927616/1‑91 (MQ=255)
                                       gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                     >  1:1486392/1‑91 (MQ=255)
                                        aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                    >  1:3859748/1‑91 (MQ=255)
                                        aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                    <  1:2480715/91‑1 (MQ=255)
                                          cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                  <  2:2256602/91‑1 (MQ=255)
                                            cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                >  1:178537/1‑91 (MQ=255)
                                            cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                >  1:109677/1‑91 (MQ=255)
                                             tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                               >  1:555169/1‑91 (MQ=255)
                                             tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                               >  2:3270693/1‑91 (MQ=255)
                                              gACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtg                                              <  2:1047755/91‑1 (MQ=255)
                                              gACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtg                                              >  2:1172260/1‑91 (MQ=255)
                                               aCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGCCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtgt                                             >  1:3849294/1‑91 (MQ=255)
                                                cAACTTATATTGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                            >  2:193471/1‑91 (MQ=255)
                                                cAACTTATATTGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                            >  2:29012/1‑91 (MQ=255)
                                                cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                            >  2:2560851/1‑91 (MQ=255)
                                                cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                            <  1:1563450/91‑1 (MQ=255)
                                                cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                            >  2:1340214/1‑91 (MQ=255)
                                                    ttATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAg                                        <  2:772390/91‑1 (MQ=255)
                                                     tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                       <  2:85696/91‑1 (MQ=255)
                                                     tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                       <  2:3193112/91‑1 (MQ=255)
                                                      atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                      >  1:2211437/1‑91 (MQ=255)
                                                       tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATa                                     <  1:1481809/91‑1 (MQ=255)
                                                        atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATTTAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                    >  2:2729188/1‑91 (MQ=255)
                                                        atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                    <  1:1053105/91‑1 (MQ=255)
                                                        atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                    <  1:2680857/91‑1 (MQ=255)
                                                         tCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaaa                                   <  1:1608896/91‑1 (MQ=255)
                                                            cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattc                                >  2:419723/1‑90 (MQ=255)
                                                            cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                >  2:529609/1‑91 (MQ=255)
                                                            cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                <  2:3156949/91‑1 (MQ=255)
                                                              aaaTAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattatt                              <  1:3250993/91‑1 (MQ=255)
                                                               aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                             <  2:1595686/91‑1 (MQ=255)
                                                                 taatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttatttat                           <  1:3696299/91‑1 (MQ=255)
                                                                       tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                     <  1:446152/91‑1 (MQ=255)
                                                                       tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                     <  1:1367692/91‑1 (MQ=255)
                                                                        cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                    <  2:1846868/91‑1 (MQ=255)
                                                                        cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                    <  1:3294362/91‑1 (MQ=255)
                                                                        cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                    >  2:3048295/1‑91 (MQ=255)
                                                                        cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTCTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                    >  2:2998910/1‑91 (MQ=255)
                                                                         aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                   <  1:3436438/91‑1 (MQ=255)
                                                                           ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                 >  1:3761983/1‑91 (MQ=255)
                                                                            tAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAt                >  1:1004662/1‑91 (MQ=255)
                                                                            tAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTAACCAt                >  2:3569084/1‑91 (MQ=255)
                                                                               aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt             >  2:2745114/1‑91 (MQ=255)
                                                                               aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt             >  1:2766886/1‑91 (MQ=255)
                                                                               aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt             <  1:454661/91‑1 (MQ=255)
                                                                                 aacTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTga           >  2:3879070/3‑91 (MQ=25)
                                                                                  gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          >  2:2864629/1‑91 (MQ=255)
                                                                                  gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          >  2:2126719/1‑91 (MQ=255)
                                                                                  acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          >  2:2271326/2‑91 (MQ=37)
                                                                                  acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          >  2:2271706/2‑91 (MQ=37)
                                                                                   cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga         >  1:3027561/1‑91 (MQ=255)
                                                                                   cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga         >  1:2578692/1‑91 (MQ=37)
                                                                                    tAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGAt        >  1:3460495/1‑91 (MQ=37)
                                                                                     aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg       >  2:2588026/1‑91 (MQ=255)
                                                                                     aaCCAATCTATATAAAGACTGGTTAGCTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg       <  2:3647000/91‑1 (MQ=37)
                                                                                     aaCCAATCTATATAAAGACTGGTTAGCTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg       >  1:3647000/1‑91 (MQ=37)
                                                                                        cAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaaa    >  1:2531531/1‑91 (MQ=255)
                                                                                         aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg   >  1:1889218/1‑91 (MQ=255)
                                                                                          gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt  <  1:764672/91‑1 (MQ=255)
                                                                                          gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt  <  1:1758783/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGT  >  USA300TCH1516_ALE/970587‑970767

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: