Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,362,249 A→T 9.2% intergenic (‑181/+37) 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,2490AT 9.2% 444.9 / 28.3 131intergenic (‑181/+37)USA300TCH1516_01271/lysP_1hypothetical protein/Lysine‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (60/59);  new base T (2/10);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 3.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 149 total aligned reads displayed.

ATATAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAC  >  USA300TCH1516_ALE/1362160‑1362339
                                                                                         |                                                                                          
atatAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCtat                                                                                           <  1:2756221/91‑1 (MQ=255)
atatAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCtat                                                                                           >  1:3983861/1‑91 (MQ=255)
atatAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCtat                                                                                           >  1:2152785/1‑91 (MQ=255)
atatAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCtat                                                                                           >  1:1964986/1‑91 (MQ=255)
  ataAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAg                                                                                         <  1:3152787/91‑1 (MQ=255)
  ataAACAGAATACTAGTCAATACACTTTATACAAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAg                                                                                         <  1:3553229/91‑1 (MQ=255)
   taAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGa                                                                                        <  1:2791490/91‑1 (MQ=255)
     aaCAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCAGTAAAAAGCGGCAAATGCTATAGATa                                                                                      >  1:3852792/1‑91 (MQ=255)
       cAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGAttat                                                                                    <  1:1694763/91‑4 (MQ=255)
       cAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGAttat                                                                                    <  1:1694746/91‑4 (MQ=255)
       cAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAAt                                                                                    <  1:2410124/91‑1 (MQ=255)
             aCTAGTTAATACACTTTATACTAGAAAAAGATATTAGAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAg                                                                              >  1:3062155/1‑91 (MQ=255)
             aCTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTGTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAg                                                                              <  1:3585469/91‑1 (MQ=255)
              cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                             <  1:3789183/91‑1 (MQ=255)
              cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                             >  1:3275922/1‑91 (MQ=255)
              cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                             >  1:677483/1‑91 (MQ=255)
              cTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGc                                                                             >  1:732616/1‑91 (MQ=255)
               tAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCa                                                                            <  1:1482061/91‑1 (MQ=255)
               tAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCa                                                                            >  1:1514240/1‑91 (MQ=255)
                aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCAt                                                                           <  1:3169975/91‑1 (MQ=37)
                aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAATGCAt                                                                           >  1:3461809/1‑91 (MQ=255)
                aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                           >  1:3461631/1‑91 (MQ=255)
                aGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCAt                                                                           <  1:3911212/91‑1 (MQ=255)
                  ttAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCt                                                                         <  1:2929473/91‑1 (MQ=255)
                   tAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTg                                                                        >  1:1389657/1‑91 (MQ=255)
                    aaTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGc                                                                       <  1:4105932/91‑1 (MQ=255)
                     aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                      <  1:698011/91‑1 (MQ=25)
                     aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                      <  1:697950/91‑1 (MQ=25)
                     aTACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAGATGCTTTAGATTATCTATAGCATTTGcc                                                                      <  1:1378486/91‑1 (MQ=25)
                        cacTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTa                                                                   >  1:3252012/1‑91 (MQ=255)
                         acTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTAt                                                                  <  1:3112457/91‑1 (MQ=255)
                          cTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATg                                                                 >  1:2873219/1‑91 (MQ=255)
                          cTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATg                                                                 >  1:710611/1‑91 (MQ=255)
                              ataCTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGaa                                                             >  1:3143510/1‑91 (MQ=255)
                                 cTAGAAAAAGATATTAATGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                          <  1:299874/91‑1 (MQ=255)
                                 cTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAattt                                                          <  1:1146286/91‑1 (MQ=255)
                                   aGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttat                                                        <  1:1886963/91‑1 (MQ=255)
                                    gAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatt                                                       >  1:1123804/1‑91 (MQ=255)
                                      aaaaGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttattta                                                     <  1:1152296/91‑1 (MQ=255)
                                         aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttatttattt                                                  >  1:1506238/1‑91 (MQ=255)
                                         aGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAatttagttattt                                                  <  1:3225464/91‑1 (MQ=255)
                                          gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                 >  1:3717719/1‑91 (MQ=255)
                                          gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                 >  1:255746/1‑91 (MQ=255)
                                          gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                 >  1:1697423/1‑91 (MQ=255)
                                          gATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAtttt                                                 <  1:2451364/91‑1 (MQ=255)
                                           atatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTAttttt                                                >  1:3359013/1‑91 (MQ=255)
                                            tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGTAAATGCTATAGATGATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                               >  1:3855242/1‑91 (MQ=255)
                                            tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                               <  1:2335479/91‑1 (MQ=255)
                                            tatTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTg                                               <  1:2337361/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:2880274/1‑91 (MQ=255)
                                                aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                           >  1:1022567/1‑91 (MQ=255)
                                                aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                           <  1:2670208/91‑1 (MQ=255)
                                                aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                           >  1:2133569/1‑91 (MQ=255)
                                                aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                           >  1:528871/1‑91 (MQ=255)
                                                aaaGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGAGAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTg                                           >  1:3292612/1‑91 (MQ=255)
                                                      tatCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATaa                                     <  1:1558187/91‑1 (MQ=255)
                                                      tatCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATaa                                     <  1:2465589/91‑1 (MQ=255)
                                                        tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                   >  1:308727/1‑91 (MQ=255)
                                                        tCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATc                                   >  1:1874018/1‑91 (MQ=255)
                                                         cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                  <  1:1261353/91‑1 (MQ=255)
                                                         cTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCg                                  >  1:314658/1‑91 (MQ=255)
                                                           gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                <  1:178789/91‑1 (MQ=255)
                                                           gCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCgtg                                <  1:2854467/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)
                                                                acacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtt                           <  1:3190475/91‑1 (MQ=255)
                                                                 cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                          <  1:3674848/91‑1 (MQ=255)
                                                                 cacCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGttt                          >  1:2470728/1‑91 (MQ=255)
                                                                  acCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGtttt                         <  1:1389171/91‑1 (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:3897220/91‑1 (MQ=255)
                                                                       tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                    <  1:3561658/91‑1 (MQ=255)
                                                                       tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                    >  1:4107269/1‑91 (MQ=255)
                                                                       tAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaatt                    <  1:2524636/91‑1 (MQ=255)
                                                                         aaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaa                  >  1:852508/1‑91 (MQ=255)
                                                                         aaaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaa                  >  1:587802/1‑91 (MQ=255)
                                                                          aaaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaat                 <  1:924729/91‑1 (MQ=255)
                                                                           aaGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAaattaatt                <  1:2803066/91‑1 (MQ=255)
                                                                             gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt              >  1:3810867/1‑91 (MQ=255)
                                                                             gCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGt              >  1:1245515/1‑91 (MQ=255)
                                                                               ggCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGttt            >  1:3421450/1‑91 (MQ=255)
                                                                                gCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTc           >  1:3889437/1‑91 (MQ=255)
                                                                                  aaaTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTa         >  1:2807317/1‑91 (MQ=255)
                                                                                    aTGCTTTAGATTATCTATAGCATTTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc       >  1:1367147/1‑91 (MQ=37)
                                                                                    aTGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGc       >  1:2594066/1‑91 (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)
                                                                                         ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc  >  1:2558864/1‑91 (MQ=255)
                                                                                         ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc  >  1:3389453/1‑91 (MQ=255)
                                                                                         ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc  <  1:3859469/91‑1 (MQ=255)
                                                                                         ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc  <  1:3584751/91‑1 (MQ=255)
                                                                                         ataGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAc  >  1:3784520/1‑91 (MQ=255)
                                                                                         |                                                                                          
ATATAAACAGAATACTAGTTAATACACTTTATACTAGAAAAAGATATTAAAGTATATCTGCTTTACACCATTAAAAAGCGGCAAATGCTATAGATAATCTAAAGCATCTGCCGTATGAGAATTTATTTATTTTTGATTGTCATAATCGTGTGGTTTTAAATTAATTGTTTCTAGCTTTAC  >  USA300TCH1516_ALE/1362160‑1362339

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