Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,362,255 A→T 8.6% intergenic (‑187/+31) 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,2550AT 8.6% 817.5 / 46.7 245intergenic (‑187/+31)USA300TCH1516_01271/lysP_1hypothetical protein/Lysine‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (138/86);  new base T (15/6);  total (153/92)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 250 total aligned reads displayed.

AACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTT  >  USA300TCH1516_ALE/1362165‑1362345
                                                                                          |                                                                                          
aaCAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATa                                                                                            >  1:23352/1‑91 (MQ=255)
 aCAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATaa                                                                                           >  2:2877471/1‑91 (MQ=255)
    gAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCt                                                                                        >  2:1518154/1‑91 (MQ=255)
     aaTACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTa                                                                                       >  1:1289295/1‑91 (MQ=255)
     aaTACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTa                                                                                       <  2:2045588/91‑1 (MQ=255)
        aCTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAg                                                                                    >  1:1084641/1‑91 (MQ=255)
         cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                                   >  1:2853872/1‑91 (MQ=255)
         cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                                   >  2:2813961/1‑91 (MQ=255)
           aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                                 >  1:2984378/1‑91 (MQ=255)
           aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                                 <  2:2276934/91‑1 (MQ=255)
            gTTAATACACATTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAACAGCGGCAAATGCTATAGATAATCTAAAGCATc                                                                                >  2:3084387/1‑91 (MQ=255)
             ttAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCAttt                                                                               <  1:3401480/91‑3 (MQ=37)
             ttAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCt                                                                               <  2:471908/91‑1 (MQ=255)
              tAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTg                                                                              <  1:828428/91‑1 (MQ=255)
                 tACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCg                                                                           >  1:693557/1‑91 (MQ=255)
                 tACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCg                                                                           >  1:2373895/1‑91 (MQ=255)
                      tttATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATga                                                                      >  1:144792/1‑91 (MQ=255)
                       ttATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATgag                                                                     <  1:3458671/91‑1 (MQ=255)
                         ataCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGaa                                                                   <  2:404262/91‑1 (MQ=255)
                         ataCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGaa                                                                   <  2:851468/91‑1 (MQ=255)
                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                >  2:2160610/1‑91 (MQ=255)
                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                <  2:1254552/91‑1 (MQ=255)
                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                >  2:1219729/1‑91 (MQ=255)
                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                >  2:2091746/1‑91 (MQ=255)
                            cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                                >  1:778655/1‑91 (MQ=255)
                            cTAGAAAAAGATATCAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGCTAATCTAAAGCATCCGCCGTATGAGAattt                                                                >  2:275009/1‑91 (MQ=255)
                             tAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattta                                                               <  1:761152/91‑1 (MQ=255)
                             tAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattta                                                               >  1:640457/1‑91 (MQ=255)
                             tAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattta                                                               <  2:1783538/91‑1 (MQ=255)
                             tAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattta                                                               <  1:407948/91‑1 (MQ=255)
                               gAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatt                                                             >  2:2901529/1‑91 (MQ=255)
                                  aaaGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttat                                                          <  1:715352/91‑1 (MQ=255)
                                    aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttattt                                                        <  1:2327008/91‑1 (MQ=255)
                                    aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttattt                                                        <  1:3551118/91‑1 (MQ=255)
                                    aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttattt                                                        <  1:3578294/91‑1 (MQ=255)
                                     gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                       >  1:1131445/1‑91 (MQ=255)
                                     gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                       <  2:3567382/91‑1 (MQ=255)
                                      atatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAttttt                                                      >  1:1259360/1‑91 (MQ=255)
                                       tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                                     <  1:100712/91‑1 (MQ=255)
                                        atTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGa                                                    >  2:2735903/1‑91 (MQ=255)
                                        atTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGa                                                    >  2:2314308/1‑91 (MQ=255)
                                         ttAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGAt                                                   >  2:634420/1‑91 (MQ=255)
                                          tAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGAtt                                                  <  1:3647816/91‑1 (MQ=255)
                                          tAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGAtt                                                  >  2:692981/1‑91 (MQ=255)
                                           aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                                 >  1:507373/1‑91 (MQ=255)
                                            aaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGt                                                >  2:1172694/1‑91 (MQ=255)
                                            aaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGt                                                <  2:1701455/91‑1 (MQ=255)
                                              gTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCa                                              >  2:2371410/1‑91 (MQ=255)
                                              gTATATCTGCTGTACACCTGTAGGAAGCGGCAAATGCTATAGATAATCTGAAGCATCTGTCGTATGAGAATTTATTTATTTTTGATTGTCa                                              <  2:2955273/91‑1 (MQ=255)
                                                    cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                        <  1:2776472/91‑1 (MQ=255)
                                                       cTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtgt                                     >  2:1087792/1‑91 (MQ=255)
                                                        tttACACCATTAAAAAGCGGCAAATGCTATAGGTAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGtgtg                                    <  1:3301162/91‑1 (MQ=255)
                                                         ttACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGAGTgg                                   >  1:3237642/1‑91 (MQ=255)
                                                           acacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCAAAATCGTGTGGtt                                 >  2:3232219/1‑91 (MQ=255)
                                                           acacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTATTTATTTTTGATTGTCATAATCGTGTGGttt                                >  1:460026/1‑91 (MQ=255)
                                                            cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                                >  1:3874377/1‑91 (MQ=255)
                                                            cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                                >  2:1490017/1‑91 (MQ=255)
                                                             acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                               >  2:1960185/1‑91 (MQ=255)
                                                             acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                               >  2:100520/1‑91 (MQ=255)
                                                               cATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaa                             >  1:1713296/1‑91 (MQ=255)
                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            >  1:1619695/1‑91 (MQ=255)
                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            >  2:669099/1‑91 (MQ=255)
                                                                aTTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTaaa                            >  1:2152604/1‑91 (MQ=255)
                                                                 ttAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaat                           <  2:1289295/91‑1 (MQ=255)
                                                                 ttAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaat                           >  1:1902114/1‑91 (MQ=255)
                                                                 ttAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaat                           >  2:2415618/1‑91 (MQ=255)
                                                                  tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                          <  2:3742485/91‑1 (MQ=255)
                                                                   aaaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatta                         >  1:2042706/1‑91 (MQ=255)
                                                                   aaaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatta                         <  1:2588763/91‑1 (MQ=255)
                                                                     aaaGCGGCATATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                       <  1:3722150/91‑1 (MQ=255)
                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      >  2:2662818/1‑91 (MQ=255)
                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      >  2:3702802/1‑91 (MQ=255)
                                                                      aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                      >  1:916617/1‑91 (MQ=255)
                                                                        gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt                    >  1:3051429/1‑91 (MQ=255)
                                                                         cGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGtt                   >  1:3190502/1‑91 (MQ=255)
                                                                          ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt                  <  1:206731/91‑1 (MQ=255)
                                                                            cAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCt                >  2:3658336/1‑91 (MQ=255)
                                                                            cAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCt                >  1:268453/1‑91 (MQ=255)
                                                                             agaTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa               >  1:466537/3‑91 (MQ=37)
                                                                             agaTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa               >  1:2458652/3‑91 (MQ=37)
                                                                             aaaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa               >  2:2517174/1‑91 (MQ=255)
                                                                              gaTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTTATTGTTTCTAg              >  2:3193275/2‑91 (MQ=37)
                                                                              aaTGCTATAGATAATCTAAAGCATCTGCCGTGTGGGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAg              >  1:1741072/1‑91 (MQ=255)
                                                                               aTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc             >  2:1320183/1‑91 (MQ=38)
                                                                                tGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCt            >  1:214251/1‑91 (MQ=255)
                                                                                tGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCt            <  1:1514521/91‑1 (MQ=255)
                                                                                 gCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCtt           <  1:116494/91‑1 (MQ=255)
                                                                                  ctttaGATTATCTATAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt          >  2:3629570/4‑91 (MQ=255)
                                                                                  cTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCttt          <  2:570004/91‑1 (MQ=255)
                                                                                   tttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc             <  2:1839552/85‑1 (MQ=38)
                                                                                   tttaGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         >  2:1417134/3‑91 (MQ=38)
                                                                                   tataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTa         <  2:886480/91‑1 (MQ=255)
                                                                                     taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa       <  2:3365262/91‑1 (MQ=255)
                                                                                     taGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACa       <  1:3203181/91‑1 (MQ=255)
                                                                                      aGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaa      >  2:229290/1‑91 (MQ=255)
                                                                                       gATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACaaa     >  1:110108/1‑91 (MQ=255)
                                                                                        aTAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAt    >  2:1670149/1‑91 (MQ=255)
                                                                                         tAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAtt   <  2:3587693/91‑1 (MQ=255)
                                                                                         tAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAtt   <  1:1779250/91‑1 (MQ=255)
                                                                                          aaTCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAAttt  <  2:2093749/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTACAAATTT  >  USA300TCH1516_ALE/1362165‑1362345

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: