Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 367,644 G→T 5.6% intergenic (‑96/‑62) nanA ← / → bglK N‑acetylneuraminate lyase/Beta‑glucoside kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE367,6440GT 5.6% 742.2 / 22.1 218intergenic (‑96/‑62)nanA/bglKN‑acetylneuraminate lyase/Beta‑glucoside kinase
Reads supporting (aligned to +/- strand):  ref base G (87/117);  new base T (8/5);  total (96/122)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 222 total aligned reads displayed.

GTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGG  >  USA300TCH1516_ALE/367555‑367734
                                                                                         |                                                                                          
gTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGt                                                                                           <  2:724556/91‑1 (MQ=255)
gTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGt                                                                                           <  1:3225833/91‑1 (MQ=255)
   ccTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCat                                                                                        <  1:791543/91‑1 (MQ=255)
      ccAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATAtt                                                                                     >  2:703541/1‑91 (MQ=255)
       cAAAAATCATTGGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATAttt                                                                                    <  1:3589191/91‑1 (MQ=255)
        aaaaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTa                                                                                   <  1:71898/91‑1 (MQ=255)
        aaaaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTa                                                                                   <  1:2018015/91‑1 (MQ=255)
         aaaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTAt                                                                                  >  1:1152213/1‑91 (MQ=255)
         aaaaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTAt                                                                                  >  2:1289844/1‑91 (MQ=255)
           aaTCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGt                                                                                <  2:3104691/91‑1 (MQ=255)
             tCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGtat                                                                              <  1:835461/91‑1 (MQ=255)
              cATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGtata                                                                             >  2:2845942/1‑91 (MQ=255)
               aTTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaa                                                                            <  2:2877690/91‑1 (MQ=255)
               aTTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaa                                                                            <  1:3730778/91‑1 (MQ=255)
                 ttGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaatt                                                                          <  2:1365318/91‑1 (MQ=255)
                  tGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaatta                                                                         <  1:400018/91‑1 (MQ=255)
                  tGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaatta                                                                         >  1:1512971/1‑91 (MQ=255)
                   gATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaattaa                                                                        <  2:1249623/91‑1 (MQ=255)
                   gATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTAtaattaa                                                                        >  2:1978776/1‑91 (MQ=255)
                      aataTAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGa                                                                     <  2:3856233/91‑1 (MQ=255)
                         ataGTTTACAGCTATAATTGTAAACGCTAGCATAAAAGGTAACAATATCTTTGTGAAAATAGTAGTCATAGCTATGGATAATTAGTGaaaa                                                                  >  2:1888225/1‑91 (MQ=255)
                           aGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATg                                                                <  1:2590043/91‑1 (MQ=255)
                            gTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGt                                                               <  1:1638338/91‑1 (MQ=255)
                             tttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGtt                                                              <  2:3491262/91‑1 (MQ=255)
                             tttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGtt                                                              <  1:1690982/91‑1 (MQ=255)
                             tttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGtt                                                              <  2:107617/91‑1 (MQ=255)
                             tttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGtt                                                              <  2:232611/91‑1 (MQ=255)
                              ttACAGCTATAATTGTAAACGCTATCTTAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGttt                                                             <  1:483942/91‑1 (MQ=255)
                              ttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGttt                                                             <  2:2546672/91‑1 (MQ=255)
                              ttACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGttt                                                             <  1:198928/91‑1 (MQ=255)
                               tACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGtttt                                                            >  2:3309321/1‑91 (MQ=255)
                                aCAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGttttt                                                           >  1:2728385/1‑91 (MQ=255)
                                 cAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAAATAAAAAAAATGTTTTTc                                                          <  2:1386568/91‑1 (MQ=255)
                                  aGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCa                                                         <  2:1639163/91‑1 (MQ=255)
                                   gCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCaa                                                        <  2:3214464/91‑1 (MQ=255)
                                   gCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCaa                                                        <  2:3226124/91‑1 (MQ=255)
                                     tataATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCaaaa                                                      >  2:3950371/1‑91 (MQ=255)
                                       taATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATc                                                    >  2:1395498/1‑91 (MQ=255)
                                       taATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATc                                                    <  2:2336387/91‑1 (MQ=255)
                                       taATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATc                                                    <  1:1672780/91‑1 (MQ=255)
                                        aaTTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCa                                                   <  1:1093961/91‑1 (MQ=255)
                                         aTTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCaa                                                  >  1:550111/1‑91 (MQ=255)
                                         aTTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCaa                                                  >  2:1322960/1‑91 (MQ=255)
                                           tGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTTATGAAAATGTTTTTCAAAATCAATa                                                >  1:3689230/1‑91 (MQ=255)
                                           tGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCtata                                                >  2:1387817/1‑91 (MQ=255)
                                            gTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAg                                               >  1:216551/1‑91 (MQ=255)
                                            gTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAg                                               <  2:112651/91‑1 (MQ=255)
                                            gTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAg                                               <  2:3527225/91‑1 (MQ=255)
                                             tAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGa                                              <  2:110725/91‑1 (MQ=255)
                                              aaaCTCTATCATAGAATGTAACGATGTCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAATTCAATAGaa                                             <  2:1933025/91‑1 (MQ=255)
                                              aaaCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaa                                             <  2:2313196/91‑1 (MQ=255)
                                              aaaCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaa                                             <  2:2477003/91‑1 (MQ=255)
                                              aaaCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaa                                             >  1:1394789/1‑91 (MQ=255)
                                               aaCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaaa                                            >  2:122696/1‑91 (MQ=255)
                                               aaCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGaaa                                            >  2:3053094/1‑91 (MQ=255)
                                                aCGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAAt                                           >  2:3083893/1‑91 (MQ=255)
                                                 cGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATg                                          >  1:2821072/1‑91 (MQ=255)
                                                   cTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGga                                        <  2:808143/91‑1 (MQ=255)
                                                   cTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGga                                        <  1:2345956/91‑1 (MQ=255)
                                                     aTCATAAAATGTAACAATATCTTTTTGAAAATTTTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagt                                      >  2:1392887/1‑91 (MQ=255)
                                                     aTCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagt                                      <  1:1340010/91‑1 (MQ=255)
                                                     aTCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagt                                      >  1:459098/1‑91 (MQ=255)
                                                     aTCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagt                                      <  2:2125566/91‑1 (MQ=255)
                                                       cATAAAATGTAACAATATCTTTTTTTTTTTTTTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagtga                                    >  2:1895192/1‑91 (MQ=255)
                                                        aTAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagtgag                                   <  2:3542413/91‑1 (MQ=255)
                                                         tAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagtgagt                                  <  1:600627/91‑1 (MQ=255)
                                                         tAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGgagtgagt                                  >  2:2185150/1‑91 (MQ=255)
                                                            aaTGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAAAGGAGTGAGTAAg                               <  2:2658877/91‑1 (MQ=255)
                                                             aTGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAgg                              <  2:1003054/91‑1 (MQ=255)
                                                              tGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGgt                             <  2:561546/91‑1 (MQ=255)
                                                               gTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAggag                            >  2:3362213/1‑89 (MQ=255)
                                                               gTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGgtg                            >  1:25079/1‑91 (MQ=255)
                                                               gTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGgtg                            <  2:1654179/91‑1 (MQ=255)
                                                                 aaCAATATCTTTTTTTTTTTTTTTTTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTa                          >  2:628444/1‑91 (MQ=255)
                                                                 aaCAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTa                          >  2:3407047/1‑91 (MQ=255)
                                                                  aCAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTAt                         >  1:2674825/1‑91 (MQ=255)
                                                                  aCAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTAt                         >  2:1141490/1‑91 (MQ=255)
                                                                   cAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTAtt                        >  1:2910706/1‑91 (MQ=255)
                                                                    aaTATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTa                       >  2:916161/1‑91 (MQ=255)
                                                                     atatCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTAc                      >  1:1845150/1‑91 (MQ=255)
                                                                         cTTTTTGAAAATTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCg                  <  1:2073408/91‑1 (MQ=255)
                                                                          tttttGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGc                 <  1:3408477/91‑1 (MQ=255)
                                                                            tttGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCaa               <  2:2980217/91‑1 (MQ=255)
                                                                            tttGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCaa               <  2:846689/91‑1 (MQ=255)
                                                                             ttGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAAt              <  2:861301/91‑1 (MQ=255)
                                                                              tGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATc             <  1:652040/91‑1 (MQ=255)
                                                                               gAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCg            <  2:2511061/91‑1 (MQ=255)
                                                                                aaaaTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGGGTAAGGTGTATTACATCGCAATCGa           <  2:448308/91‑1 (MQ=255)
                                                                                aaaaTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGa           <  2:3723161/91‑1 (MQ=255)
                                                                                 aaaTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGat          <  2:3600722/91‑1 (MQ=255)
                                                                                  aaTTGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGata         <  1:2474410/91‑1 (MQ=255)
                                                                                  aaTTGTAGTCATCTTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGata         >  1:3060540/1‑91 (MQ=255)
                                                                                  aaTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGata         <  2:3637635/91‑1 (MQ=255)
                                                                                   aTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGatat        <  2:1827947/91‑1 (MQ=255)
                                                                                   aTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGatat        >  1:3298254/1‑91 (MQ=255)
                                                                                   aTTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGatat        <  1:2393183/91‑1 (MQ=255)
                                                                                     tGTATTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTg      >  1:3345272/1‑91 (MQ=255)
                                                                                         gTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAgg  >  2:61473/1‑91 (MQ=255)
                                                                                         gTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAgg  >  1:1573254/1‑91 (MQ=255)
                                                                                         gTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAgg  <  1:3231920/91‑1 (MQ=255)
                                                                                         |                                                                                          
GTCCCTCCAAAAATCATTTGATAATATAGTTTACAGCTATAATTGTAAACGCTATCATAAAATGTAACAATATCTTTTTGAAAATTGTAGTCATATTTATGTATAATTAATGAAAATGTTTTTCAAAATCAATAGAAATGGAGTGAGTAAGGTGTATTACATCGCAATCGATATTGGAGG  >  USA300TCH1516_ALE/367555‑367734

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