Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,690 T→C 11.5% intergenic (+36/‑293) 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,6900TC11.5% 680.9 / 61.6 243intergenic (+36/‑293)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (133/82);  new base C (13/15);  total (146/97)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 249 total aligned reads displayed.

AGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGATA  >  USA300TCH1516_ALE/970600‑970779
                                                                                          |                                                                                         
aGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAt                                                                                           >  2:3143837/1‑91 (MQ=255)
  ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTg                                                                                         <  1:3281342/91‑1 (MQ=255)
     aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                                      >  2:2245337/1‑91 (MQ=255)
       gTGACTATTGATAATATTGGAACGCTGACGACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                                    <  1:3108633/91‑1 (MQ=37)
        tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                                   >  2:3784894/1‑91 (MQ=255)
         gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtt                                                                                  >  1:1416446/1‑91 (MQ=255)
         gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgct                                                                                  <  2:156199/91‑3 (MQ=21)
          aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                                 >  1:3733319/1‑91 (MQ=255)
          aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                                 >  1:652692/1‑91 (MQ=255)
          aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAgctt                                                                                 <  2:243057/91‑4 (MQ=37)
           cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                                >  1:2189873/1‑91 (MQ=255)
           cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                                >  1:1595686/1‑91 (MQ=255)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                              >  2:77612/1‑91 (MQ=255)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                              <  1:985428/91‑1 (MQ=21)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                              <  1:895093/91‑1 (MQ=21)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                              <  1:3815608/91‑1 (MQ=21)
              ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                             <  1:1179463/91‑2 (MQ=21)
               tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                            <  2:3005202/91‑1 (MQ=255)
               tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                            <  1:202060/91‑1 (MQ=255)
               tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGCTTAGTTTTTTCt                                                                            <  1:2990491/91‑1 (MQ=255)
                 ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                          >  2:1412858/1‑91 (MQ=255)
                 ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGctttttaaat                                                                          <  2:803368/91‑5 (MQ=17)
                     tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                                      >  1:280788/1‑91 (MQ=255)
                      atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtttt                                                                     <  1:3663259/91‑1 (MQ=255)
                        tGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                                   <  2:2189743/91‑1 (MQ=255)
                         ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                                  >  2:3642451/1‑91 (MQ=255)
                         ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                                  <  2:1447830/91‑1 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                                 >  2:1927616/1‑91 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                                 <  1:3362658/91‑1 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                                 >  1:1486392/1‑91 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGGTAGTTTTTTCTTGCTTTTCTac                                                                 >  1:3554187/1‑90 (MQ=255)
                           aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                                >  1:3859954/1‑91 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                               >  1:445806/1‑91 (MQ=255)
                             cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGGTTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                              >  2:2784629/1‑91 (MQ=255)
                             cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                              >  2:2933632/1‑91 (MQ=255)
                               cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                            >  2:31727/1‑91 (MQ=255)
                               cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                            >  1:109677/1‑91 (MQ=255)
                                tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                           >  1:555169/1‑91 (MQ=255)
                                tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                           >  2:1373546/1‑91 (MQ=255)
                                tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                           >  2:3270693/1‑91 (MQ=255)
                                 gACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtg                                                          >  2:1172260/1‑91 (MQ=255)
                                   cAACTTATATTGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                        >  2:193471/1‑91 (MQ=255)
                                   cAACTTATATTGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                        >  2:29011/1‑91 (MQ=255)
                                   cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                        >  2:2560851/1‑91 (MQ=255)
                                   cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                        >  2:1340214/1‑91 (MQ=255)
                                       ttATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAg                                                    <  2:772390/91‑1 (MQ=255)
                                        tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                   <  2:2834698/91‑1 (MQ=255)
                                        tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                   <  1:881176/91‑1 (MQ=255)
                                        tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                   <  2:85696/91‑1 (MQ=255)
                                        tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                   <  2:3193112/91‑1 (MQ=255)
                                         atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                                  <  1:3857896/91‑1 (MQ=255)
                                          tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATa                                                 <  1:1481809/91‑1 (MQ=255)
                                           atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                                >  2:2728845/1‑91 (MQ=255)
                                           atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                                >  2:2729238/1‑91 (MQ=255)
                                           atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                                <  1:1053105/91‑1 (MQ=255)
                                           atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                                <  1:2680857/91‑1 (MQ=255)
                                           atCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATaa                                                >  1:756470/1‑91 (MQ=255)
                                              gCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatt                                             >  2:2261169/1‑91 (MQ=255)
                                              gCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatt                                             >  1:44933/1‑91 (MQ=255)
                                               cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                            >  2:529609/1‑91 (MQ=255)
                                               cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                            >  1:3440585/1‑91 (MQ=255)
                                                     aatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTtata                                      >  1:2703347/1‑91 (MQ=255)
                                                       ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                                    >  1:2983600/1‑91 (MQ=255)
                                                        tATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATg                                   >  1:1681785/1‑91 (MQ=255)
                                                         aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                                  <  2:1647757/91‑1 (MQ=255)
                                                         aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                                  <  2:3425070/91‑1 (MQ=255)
                                                         aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                                  >  1:2462629/1‑91 (MQ=255)
                                                         aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                                  <  2:414100/91‑1 (MQ=255)
                                                           cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                                >  2:1670573/1‑91 (MQ=255)
                                                           cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                                >  2:3048295/1‑91 (MQ=255)
                                                           cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                                <  2:1846868/91‑1 (MQ=255)
                                                              ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                             >  1:3761983/1‑91 (MQ=255)
                                                               tAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAt                            >  1:1004662/1‑91 (MQ=255)
                                                                  aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                         >  2:2745114/1‑91 (MQ=255)
                                                                     gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag                      >  1:2214328/1‑91 (MQ=255)
                                                                     acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag                      >  1:1550334/2‑91 (MQ=37)
                                                                     acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag                      >  2:2271340/2‑91 (MQ=37)
                                                                      cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                     >  1:3027561/1‑91 (MQ=255)
                                                                      cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                     >  1:1284153/1‑91 (MQ=255)
                                                                      cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                     <  2:2360522/91‑1 (MQ=255)
                                                                      cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                     >  2:1764338/1‑91 (MQ=37)
                                                                      cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                     >  1:2578692/1‑91 (MQ=37)
                                                                         aCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGa                  >  1:2117359/1‑91 (MQ=255)
                                                                           cAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaaa                >  1:2531531/1‑91 (MQ=255)
                                                                            aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg               >  1:1889263/1‑91 (MQ=255)
                                                                            aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg               >  1:1889218/1‑91 (MQ=255)
                                                                             gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt              <  1:764672/91‑1 (MQ=255)
                                                                             gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt              <  1:1759052/91‑1 (MQ=255)
                                                                              tCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg             >  2:1424943/1‑91 (MQ=255)
                                                                               cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa            >  1:1362674/1‑91 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAGGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAat         >  2:520939/1‑91 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAat         >  2:1040645/1‑91 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTCGAGATGAAAGTGAAat         >  2:3710859/1‑91 (MQ=255)
                                                                                   atataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAata        <  2:2123411/91‑1 (MQ=255)
                                                                                    tataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAatat       <  2:555169/91‑1 (MQ=255)
                                                                                     ataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATAtt      <  2:515047/91‑1 (MQ=255)
                                                                                     ataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATAtt      <  2:1498260/91‑1 (MQ=255)
                                                                                     ataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATAtt      <  1:2309213/91‑1 (MQ=255)
                                                                                       aGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGa    >  1:3163171/1‑91 (MQ=255)
                                                                                       aGATTGGTTAGTTTTTTCTTGCTTTGCTAAAGAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGa    <  2:1058310/91‑1 (MQ=255)
                                                                                          |                                                                                         
AGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGATA  >  USA300TCH1516_ALE/970600‑970779

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: