Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,362,267 C→T 7.3% intergenic (‑199/+19) USA300TCH1516_01271 ← / ← lysP_1 hypothetical protein/Lysine‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,362,2670CT 7.3% 489.6 / 20.0 151intergenic (‑199/+19)USA300TCH1516_01271/lysP_1hypothetical protein/Lysine‑specific permease
Reads supporting (aligned to +/- strand):  ref base C (69/71);  new base T (6/5);  total (75/76)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 153 total aligned reads displayed.

TTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGAAT  >  USA300TCH1516_ALE/1362178‑1362357
                                                                                         |                                                                                          
ttAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCt                                                                                           <  1:2929473/91‑1 (MQ=255)
 tAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTg                                                                                          >  1:1389657/1‑91 (MQ=255)
  aaTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGc                                                                                         <  1:4105932/91‑1 (MQ=255)
   aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                                        <  1:697936/91‑1 (MQ=25)
   aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                                        <  1:698011/91‑1 (MQ=25)
   aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                                        <  1:698024/91‑1 (MQ=25)
     acacTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGt                                                                                      <  1:459197/91‑1 (MQ=255)
       acTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTAt                                                                                    <  1:3112457/91‑1 (MQ=255)
        cTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATg                                                                                   >  1:2873219/1‑91 (MQ=255)
        cTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATg                                                                                   <  1:1695878/91‑1 (MQ=255)
          ttATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATgag                                                                                 <  1:2587289/91‑1 (MQ=255)
            ataCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGaa                                                                               >  1:3467090/1‑91 (MQ=255)
               cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                            <  1:1146286/91‑1 (MQ=255)
                 aGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttat                                                                          <  1:1886963/91‑1 (MQ=255)
                    aaaaGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttattta                                                                       <  1:1152296/91‑1 (MQ=255)
                       aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttagttattt                                                                    <  1:3225464/91‑1 (MQ=255)
                        gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                                   >  1:255746/1‑91 (MQ=255)
                        gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                                   <  1:2451364/91‑1 (MQ=255)
                        gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                                   >  1:3720295/1‑91 (MQ=255)
                          tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                                                 <  1:2335479/91‑1 (MQ=255)
                          tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGAAAATGCTATAGATCATCTAAAGCATCTGCCGTATGAGAATTTATTTACTTTTg                                                                 >  1:3857455/1‑91 (MQ=255)
                           atTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGa                                                                <  1:3037784/91‑1 (MQ=255)
                             tAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGAtt                                                              <  1:3915879/91‑1 (MQ=255)
                              aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                             >  1:4104132/1‑91 (MQ=255)
                              aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                             >  1:1022567/1‑91 (MQ=255)
                              aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGAGAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                             >  1:3292612/1‑91 (MQ=255)
                               aaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCCGCCGTATGAGAATTTATTTATTTTTGATTGt                                                            <  1:2970533/91‑1 (MQ=255)
                                  tataTCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCAt                                                         <  1:3647788/91‑1 (MQ=255)
                                      tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                                     >  1:308727/1‑91 (MQ=255)
                                      tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                                     >  1:1874018/1‑91 (MQ=255)
                                      tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                                     <  1:1552157/91‑1 (MQ=255)
                                       cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                                    >  1:314723/1‑91 (MQ=255)
                                       cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                                    >  1:314658/1‑91 (MQ=255)
                                       cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                                    <  1:1261353/91‑1 (MQ=255)
                                         gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                                  <  1:178789/91‑1 (MQ=255)
                                          cTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtgt                                                 <  1:3246334/91‑1 (MQ=255)
                                           tttACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGtgtg                                                <  1:3936408/91‑1 (MQ=255)
                                            ttACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTgg                                               >  1:2498239/1‑91 (MQ=255)
                                             tACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGt                                              >  1:1603081/1‑91 (MQ=255)
                                               cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                                            >  1:2470728/1‑91 (MQ=255)
                                                acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                                           <  1:1389171/91‑1 (MQ=255)
                                                acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                                           <  1:2844424/91‑1 (MQ=255)
                                                 ccATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTa                                          >  1:3894696/1‑91 (MQ=255)
                                                 ccATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTa                                          >  1:400118/1‑91 (MQ=255)
                                                  cATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaa                                         >  1:1088112/1‑91 (MQ=255)
                                                   aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                                        <  1:2174503/91‑1 (MQ=255)
                                                   aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                                        <  1:73690/91‑1 (MQ=255)
                                                     tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                                      <  1:2524636/91‑1 (MQ=255)
                                                       aaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaa                                    >  1:587802/1‑91 (MQ=255)
                                                        catacGGCAGATGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                                   >  1:2687099/5‑91 (MQ=21)
                                                        aaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                                   <  1:924729/91‑1 (MQ=255)
                                                         aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                                  >  1:2024162/1‑91 (MQ=255)
                                                         aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                                  <  1:2803065/91‑1 (MQ=255)
                                                           gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt                                >  1:1245515/1‑91 (MQ=255)
                                                           gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt                                >  1:3810867/1‑91 (MQ=255)
                                                            cGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGtt                               >  1:1363194/1‑91 (MQ=255)
                                                             ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                              <  1:2460102/91‑1 (MQ=255)
                                                             ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                              >  1:3421450/1‑91 (MQ=255)
                                                              gCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTc                             >  1:1487307/1‑91 (MQ=255)
                                                              gCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTc                             >  1:3889437/1‑91 (MQ=255)
                                                                aaaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa                           >  1:2807317/1‑91 (MQ=255)
                                                                 aaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAg                          >  1:3802361/1‑91 (MQ=255)
                                                                  aTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc                         >  1:1367147/1‑91 (MQ=37)
                                                                  aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTGTTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc                         <  1:3288553/91‑1 (MQ=255)
                                                                  aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc                         >  1:2594066/1‑91 (MQ=255)
                                                                  aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc                         <  1:2181886/91‑1 (MQ=255)
                                                                   tGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTCTTTCTAGCt                        <  1:2092406/91‑1 (MQ=255)
                                                                     cTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt                      <  1:648488/91‑1 (MQ=255)
                                                                     cTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt                      <  1:1981511/91‑1 (MQ=255)
                                                                      tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa                     >  1:1486388/1‑91 (MQ=255)
                                                                      tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa                     >  1:2876222/1‑91 (MQ=255)
                                                                       ttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    >  1:1731862/2‑91 (MQ=39)
                                                                       ttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    >  1:3063126/2‑91 (MQ=39)
                                                                       ttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    >  1:3063147/2‑91 (MQ=39)
                                                                       ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    <  1:3584751/91‑1 (MQ=255)
                                                                       ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    >  1:3389453/1‑91 (MQ=255)
                                                                       ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    >  1:2558864/1‑91 (MQ=255)
                                                                       ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc                    <  1:3858503/91‑1 (MQ=255)
                                                                        taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa                   <  1:2682773/91‑1 (MQ=255)
                                                                        taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa                   >  1:3563691/1‑91 (MQ=255)
                                                                         cgATAAGCTAACGCATCTGCCGTATGAGAATGTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa                  <  1:1509776/90‑1 (MQ=255)
                                                                         aGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa                  <  1:3023419/91‑1 (MQ=255)
                                                                            tAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAtt               <  1:426216/91‑1 (MQ=255)
                                                                             aaTCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttt              >  1:536260/1‑91 (MQ=255)
                                                                              aTCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAtttt             >  1:2317023/1‑91 (MQ=255)
                                                                              aTCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTAACAAAtttt             >  1:3724118/1‑91 (MQ=255)
                                                                                cTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgt           <  1:3375850/91‑1 (MQ=255)
                                                                                cTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgt           <  1:3015427/91‑1 (MQ=255)
                                                                                   aaGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgtttt        >  1:3676276/1‑91 (MQ=255)
                                                                                   aaGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgtttt        <  1:3576041/91‑1 (MQ=255)
                                                                                    aGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttttgttttg       <  1:1999220/91‑1 (MQ=255)
                                                                                      cATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTtgtg     <  1:369658/91‑1 (MQ=255)
                                                                                      cATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTtgtg     >  1:885652/1‑91 (MQ=255)
                                                                                       aTCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGa    >  1:721984/1‑91 (MQ=255)
                                                                                        tCTGCCGTATGAGAATTTATTTATTTTTTATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGaa   <  1:2131904/91‑1 (MQ=255)
                                                                                        tCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGaa   <  1:2131223/91‑1 (MQ=255)
                                                                                        tCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGaa   <  1:3371513/91‑1 (MQ=255)
                                                                                        tCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGaa   <  1:862452/91‑1 (MQ=255)
                                                                                         cTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGAAt  >  1:1590455/1‑91 (MQ=255)
                                                                                         cTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGAAt  >  1:2780765/1‑91 (MQ=255)
                                                                                         |                                                                                          
TTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTTTGTTTTGTGAAT  >  USA300TCH1516_ALE/1362178‑1362357

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: