Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,135,944 A→T 5.4% K146N (AAA→AAT ktrB_2 → Ktr system potassium uptake protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,135,9440AT 5.4% 703.9 / 21.1 204K146N (AAA→AATktrB_2Ktr system potassium uptake protein B
Reads supporting (aligned to +/- strand):  ref base A (91/102);  new base T (4/7);  total (95/109)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 221 total aligned reads displayed.

TGTTATTAAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAACTT  >  USA300TCH1516_ALE/2135854‑2136033
                                                                                          |                                                                                         
tgtTATTAAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCaaaa                                                                                           >  2:1224566/1‑91 (MQ=255)
 gtTATTAAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAAt                                                                                          <  1:2662243/91‑1 (MQ=255)
  ttattaAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAAtt                                                                                         <  2:75488/91‑1 (MQ=255)
  ttattaAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAAtt                                                                                         <  1:2164396/91‑1 (MQ=255)
   tattaAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAAttt                                                                                        >  1:1427248/1‑91 (MQ=255)
   tattaAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAAttt                                                                                        >  2:2043665/1‑91 (MQ=255)
       aaGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGta                                                                                    <  1:3691457/91‑1 (MQ=255)
        aGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGtat                                                                                   >  2:2508080/1‑91 (MQ=255)
           tAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAATTTTGGtataaa                                                                                <  1:3271275/91‑3 (MQ=255)
           tAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAATTTTGGtataaa                                                                                <  1:3113979/91‑3 (MQ=255)
           tAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAgg                                                                                <  1:467563/91‑1 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              >  1:75958/1‑91 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              <  2:1327488/91‑1 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              >  2:2788817/1‑91 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              >  2:2789241/1‑91 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              <  1:2136483/91‑1 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              >  2:1687348/1‑91 (MQ=255)
             aTTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCa                                                                              >  1:2668675/1‑91 (MQ=255)
              ttAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCaa                                                                             >  1:1974772/1‑91 (MQ=255)
               tAAACACTAGGCTATTTATAGTTTAGTCACTGAATTAATGGGTAAGATTTGAATGTGTTTATCTTTTATACCAAAATTTGGTATAGGCaaa                                                                            >  2:3086021/1‑91 (MQ=255)
                aaaCACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGTATAGGCAAAgg                                                                          >  2:881116/1‑91 (MQ=255)
                   cacTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGtt                                                                        >  1:1747883/1‑91 (MQ=255)
                   cacTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGtt                                                                        >  2:1405907/1‑91 (MQ=255)
                     cTTTGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTa                                                                      <  1:3406794/91‑1 (MQ=255)
                     cTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTa                                                                      >  1:2089629/1‑91 (MQ=255)
                      ttGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAt                                                                     <  1:3741922/91‑1 (MQ=255)
                       tGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAtt                                                                    <  1:553091/91‑1 (MQ=255)
                          cTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAttttt                                                                 >  2:3106709/1‑91 (MQ=255)
                          cTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAttttt                                                                 >  2:2443379/1‑91 (MQ=255)
                          cTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAttttt                                                                 >  2:1742485/1‑91 (MQ=255)
                          cTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTAttttt                                                                 >  1:2887700/1‑91 (MQ=255)
                           tATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTa                                                                <  2:2247901/91‑1 (MQ=255)
                               tataGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCt                                                            >  1:2962255/1‑91 (MQ=255)
                               tataGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCt                                                            >  2:3883647/1‑91 (MQ=255)
                               tataGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCt                                                            <  1:168054/91‑1 (MQ=255)
                               tataGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCt                                                            >  1:2724893/1‑91 (MQ=255)
                               tataGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCt                                                            <  1:3548329/91‑1 (MQ=255)
                                 taGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTa                                                          <  2:2932959/91‑1 (MQ=255)
                                    tttAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTc                                                       <  1:891908/91‑1 (MQ=255)
                                     ttAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTca                                                      <  1:3377679/91‑1 (MQ=255)
                                      tAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTcac                                                     <  1:995229/91‑1 (MQ=255)
                                       aGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTAAACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTcaca                                                    <  2:188072/91‑1 (MQ=255)
                                            cTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGt                                               <  1:1654999/91‑1 (MQ=255)
                                            cTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGt                                               <  1:2907340/91‑1 (MQ=255)
                                            cTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGt                                               <  1:1687348/91‑1 (MQ=255)
                                             tGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTa                                              >  1:1173236/1‑91 (MQ=255)
                                             tGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTa                                              <  1:2778180/91‑1 (MQ=255)
                                              gAATTAATTGGTATGATTTGTGTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTAt                                             <  2:300775/91‑1 (MQ=255)
                                               aattaattGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATc                                            <  2:2360216/91‑1 (MQ=255)
                                                 ttaattGTTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAg                                          <  1:1937701/91‑1 (MQ=255)
                                                  taattGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGc                                         >  2:2734515/1‑91 (MQ=255)
                                                    attGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCtt                                       >  2:666590/1‑91 (MQ=255)
                                                      tGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCtttt                                     >  2:1643529/1‑91 (MQ=255)
                                                       ggTATGATTTGTTTGTGTTTATCTTCTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTta                                    <  2:2303830/91‑1 (MQ=255)
                                                        gTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTtaa                                   >  1:3470995/1‑91 (MQ=255)
                                                        gTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTtaa                                   >  2:3303582/1‑91 (MQ=255)
                                                         tATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTtaat                                  >  1:847662/1‑91 (MQ=255)
                                                          aTGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTtaata                                 >  2:1531519/1‑91 (MQ=255)
                                                            gATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTaataat                               >  1:2246053/1‑91 (MQ=255)
                                                             aTTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATg                              >  2:3339351/1‑91 (MQ=255)
                                                             aTTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATg                              <  2:1730782/91‑1 (MQ=255)
                                                              tttgtttgTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGc                             >  2:452988/1‑91 (MQ=255)
                                                               ttgtttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCt                            <  1:1127455/91‑1 (MQ=255)
                                                               ttgttcgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCt                            <  2:3581715/91‑1 (MQ=255)
                                                                tgtttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGc                             <  1:614688/89‑1 (MQ=255)
                                                                 gtttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTgg                          >  1:915379/1‑91 (MQ=255)
                                                                 gtttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTgg                          >  2:571700/1‑91 (MQ=255)
                                                                 gtttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTgg                          <  2:51212/91‑1 (MQ=255)
                                                                  tttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGa                         >  1:3733864/1‑91 (MQ=255)
                                                                  tttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGa                         <  2:1380850/91‑1 (MQ=255)
                                                                   ttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAt                        <  1:1020668/91‑1 (MQ=255)
                                                                   ttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAt                        <  2:3546807/91‑1 (MQ=255)
                                                                   ttgtGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTGTTAAGCTTATTCACATCAGTATCAGCTTTCAATAATGCTGGAt                        >  1:3880038/1‑91 (MQ=255)
                                                                    tgtgTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAtt                       >  2:857576/1‑91 (MQ=255)
                                                                    tgtgTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAtt                       <  2:3489904/91‑1 (MQ=255)
                                                                    tgtgTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAtt                       <  2:1148158/91‑1 (MQ=255)
                                                                    tgtgTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAtt                       <  2:2490581/91‑1 (MQ=255)
                                                                     gtgtTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGAttt                      >  1:1929805/1‑91 (MQ=255)
                                                                      tgtTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTg                     <  1:744196/91‑1 (MQ=255)
                                                                       gtTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGc                    <  1:196848/91‑1 (MQ=255)
                                                                        tttATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGcc                   <  1:2484460/91‑1 (MQ=255)
                                                                        tttATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGcc                   <  1:1405907/91‑1 (MQ=255)
                                                                             cTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCtttt              <  1:283692/91‑1 (MQ=255)
                                                                              ttttATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCttttt             <  1:317305/91‑1 (MQ=255)
                                                                               tttATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTa            <  1:976982/91‑1 (MQ=255)
                                                                               tttATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTa            >  1:2137453/1‑91 (MQ=255)
                                                                               tttATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAACTTTTAATAATGCTGGATTTGCCCTTTTTa            >  1:3047457/1‑91 (MQ=255)
                                                                                 tataCCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAg          >  2:1186822/1‑91 (MQ=255)
                                                                                   taCCAAATTTTGGTGTAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGaa        >  2:2164332/1‑91 (MQ=255)
                                                                                   taCCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGaa        <  2:2761438/91‑1 (MQ=255)
                                                                                   taCCAAAATTTGGTATAGGCAAAGGTTTATTTATAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGaa        >  2:3802678/1‑91 (MQ=255)
                                                                                    aCCAAATTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAAt       >  1:1816462/1‑91 (MQ=255)
                                                                                    aCCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAAt       >  1:1855571/1‑91 (MQ=255)
                                                                                     ccAAATTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATa      >  1:3850896/1‑91 (MQ=255)
                                                                                       aaaaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAAc    >  1:3284107/1‑91 (MQ=255)
                                                                                       aaaaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAAc    <  2:1427248/91‑1 (MQ=255)
                                                                                       aaaaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAAc    <  2:3033402/91‑1 (MQ=255)
                                                                                        aaaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAACt   >  1:3425621/1‑91 (MQ=255)
                                                                                         aaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAACtt  >  1:3260632/1‑91 (MQ=255)
                                                                                         aaTTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAACtt  <  2:2785369/91‑1 (MQ=255)
                                                                                          |                                                                                         
TGTTATTAAGCTAATTAAACACTTGGCTATTTATAGTTTAGTCACTGAATTAATTGGTATGATTTGTTTGTGTTTATCTTTTATACCAAAATTTGGTATAGGCAAAGGTTTATTTTTAAGCTTATTCACATCAGTATCAGCTTTTAATAATGCTGGATTTGCCCTTTTTAAGAATAACTT  >  USA300TCH1516_ALE/2135854‑2136033

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: