Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,516,445 A→G 5.5% intergenic (‑96/+5) USA300TCH1516_01381 ← / ← gpsB hypothetical protein/Cell cycle protein GpsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,516,4450AG 5.5% 411.1 / 10.3 126intergenic (‑96/+5)USA300TCH1516_01381/gpsBhypothetical protein/Cell cycle protein GpsB
Reads supporting (aligned to +/- strand):  ref base A (80/39);  new base G (3/4);  total (83/43)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 127 total aligned reads displayed.

TTCCTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTA  >  USA300TCH1516_ALE/1516355‑1516534
                                                                                          |                                                                                         
ttCCTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTtttattta                                                                                           <  1:3223375/91‑1 (MQ=255)
  ccTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGt                                                                                         <  1:2970462/91‑1 (MQ=255)
  ccTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGt                                                                                         >  1:3100256/1‑91 (MQ=255)
  ccTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGt                                                                                         >  1:1906381/1‑91 (MQ=255)
  ccTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGt                                                                                         >  1:1905757/1‑91 (MQ=255)
   ctctATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTa                                                                                        >  1:3676810/1‑91 (MQ=255)
   ctctATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTa                                                                                        >  1:3676866/1‑91 (MQ=255)
   ctctATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTCTTGTATTTAGTa                                                                                        <  1:3201517/91‑1 (MQ=255)
    tctATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTaa                                                                                       >  1:2733594/1‑91 (MQ=255)
    tctATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTaa                                                                                       >  1:3822879/1‑91 (MQ=255)
        tataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattat                                                                                   >  1:2565933/1‑91 (MQ=255)
        tataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattat                                                                                   >  1:1382246/1‑91 (MQ=255)
         ataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattatt                                                                                  <  1:642629/91‑1 (MQ=255)
          taatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTAttt                                                                                 <  1:146813/91‑1 (MQ=255)
            atataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAc                                                                               >  1:76251/1‑91 (MQ=255)
             tataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAcc                                                                              >  1:2471918/1‑91 (MQ=255)
             tataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAcc                                                                              >  1:2471919/1‑91 (MQ=255)
              ataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTActa                                                                             <  1:4023304/91‑3 (MQ=255)
              ataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCa                                                                             <  1:2391383/91‑1 (MQ=255)
                  cGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATa                                                                         >  1:2503685/1‑91 (MQ=255)
                   gATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATAc                                                                        <  1:2873263/91‑1 (MQ=255)
                     ttACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAg                                                                      >  1:686598/1‑91 (MQ=255)
                      tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAgt                                                                     >  1:3094696/1‑90 (MQ=255)
                      tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAgt                                                                     >  1:3094384/1‑90 (MQ=255)
                      tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGc                                                                     <  1:3698705/91‑1 (MQ=255)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTACTAAATACaaaa                                                                    <  1:3926135/91‑4 (MQ=37)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTACTAAATACaaaa                                                                    <  1:2849317/91‑4 (MQ=37)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTACTAAATACaaaa                                                                    <  1:2849207/91‑4 (MQ=37)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                    >  1:2680315/1‑91 (MQ=255)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                    >  1:2627360/1‑91 (MQ=255)
                       aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                    >  1:2627013/1‑91 (MQ=255)
                          cAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCtttt                                                                 <  1:2858557/91‑1 (MQ=255)
                           aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                >  1:540230/1‑91 (MQ=255)
                           aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                >  1:3793414/1‑91 (MQ=255)
                           aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                >  1:3793397/1‑91 (MQ=255)
                           aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                <  1:3793298/91‑1 (MQ=255)
                           aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCctttt                                                                >  1:1399833/1‑91 (MQ=255)
                                 tCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGt                                                          <  1:4066447/91‑1 (MQ=255)
                                 tCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGt                                                          <  1:2878441/91‑1 (MQ=255)
                                  cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                         >  1:1165486/1‑91 (MQ=255)
                                  cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                         >  1:3601256/1‑91 (MQ=255)
                                  cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                         >  1:2329094/1‑91 (MQ=255)
                                  cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                         >  1:1165130/1‑91 (MQ=255)
                                  cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATATACCAAATACAGCTTTTTCTAAGtt                                                         >  1:3853038/1‑91 (MQ=255)
                                    cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                       >  1:2402491/1‑91 (MQ=255)
                                    cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                       >  1:2089925/1‑91 (MQ=255)
                                    cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                       >  1:1105648/1‑91 (MQ=255)
                                    cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGATTg                                                       >  1:3341260/1‑91 (MQ=255)
                                    cTTTTAAAATTAGAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                       >  1:3177375/1‑91 (MQ=255)
                                     ttttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGa                                                      >  1:3099133/1‑91 (MQ=255)
                                     ttttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGa                                                      >  1:3029892/1‑91 (MQ=255)
                                      tttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAATAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaa                                                     <  1:3489215/91‑1 (MQ=255)
                                       ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                    <  1:2955089/91‑1 (MQ=255)
                                       ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                    >  1:3156740/1‑91 (MQ=255)
                                              tataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGttt                                             >  1:3255208/1‑91 (MQ=255)
                                               ataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGtttt                                            <  1:172590/91‑1 (MQ=255)
                                               ataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGtttt                                            >  1:3134177/1‑91 (MQ=255)
                                                taACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTa                                           >  1:1566510/1‑91 (MQ=255)
                                                    aTAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAat                                       >  1:3263589/1‑91 (MQ=255)
                                                       gTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATc                                    <  1:1421165/91‑1 (MQ=255)
                                                          attattAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTAc                                 <  1:3802999/91‑1 (MQ=255)
                                                           ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                <  1:4033937/91‑1 (MQ=255)
                                                           ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                >  1:1902589/1‑91 (MQ=255)
                                                           ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                <  1:4033612/91‑1 (MQ=255)
                                                           ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                <  1:4033564/91‑1 (MQ=255)
                                                             attaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACatt                              <  1:2505824/91‑1 (MQ=255)
                                                               taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                            >  1:3919962/1‑91 (MQ=255)
                                                               taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                            >  1:348656/1‑91 (MQ=255)
                                                               taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                            >  1:1905928/1‑91 (MQ=255)
                                                                 gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                          >  1:2970002/1‑91 (MQ=255)
                                                                 gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                          >  1:3285268/1‑91 (MQ=255)
                                                                 gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                          >  1:2664597/1‑91 (MQ=255)
                                                                  tAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTg                         <  1:1747388/91‑1 (MQ=255)
                                                                   aaGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGa                        >  1:3255635/1‑91 (MQ=255)
                                                                     gACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAg                      >  1:2373826/1‑91 (MQ=255)
                                                                     gACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAg                      >  1:2377737/1‑91 (MQ=255)
                                                                      aCAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGa                     <  1:2606198/91‑1 (MQ=255)
                                                                         gTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGt                  >  1:1956514/1‑91 (MQ=255)
                                                                          ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                 >  1:2643275/1‑91 (MQ=255)
                                                                          ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                 >  1:693210/1‑91 (MQ=255)
                                                                          ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                 <  1:3367851/91‑1 (MQ=255)
                                                                             aaCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAGGATGTAttt              <  1:2082011/91‑1 (MQ=255)
                                                                             aaCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTAttt              <  1:2865823/91‑1 (MQ=255)
                                                                              aCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttg             >  1:1823435/1‑91 (MQ=255)
                                                                              aCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttg             <  1:4091816/91‑1 (MQ=255)
                                                                               cTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgt            <  1:2285515/91‑1 (MQ=255)
                                                                                agctGTATTTGGTAAATAATTACTAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgtt           >  1:3954567/4‑91 (MQ=37)
                                                                                  ttGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgt         >  1:1939255/1‑91 (MQ=255)
                                                                                   tGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgtt        >  1:2156302/1‑91 (MQ=255)
                                                                                   tGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgtt        >  1:2155908/1‑91 (MQ=255)
                                                                                    gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg       >  1:1753293/1‑91 (MQ=255)
                                                                                    gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg       <  1:2776788/91‑1 (MQ=255)
                                                                                    gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg       >  1:1753068/1‑91 (MQ=255)
                                                                                     tATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTtgttgttgt      <  1:1642570/91‑1 (MQ=255)
                                                                                     tATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTtgttgttgt      <  1:2915564/91‑1 (MQ=255)
                                                                                      aTTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta     >  1:1035861/1‑91 (MQ=255)
                                                                                      aTTTAGTAAGTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta     <  1:3976831/91‑1 (MQ=255)
                                                                                      aTTTAGTAAGTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta     <  1:3976705/91‑1 (MQ=255)
                                                                                         tAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGtatta  >  1:332402/1‑91 (MQ=255)
                                                                                         tAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGtatta  >  1:2543468/1‑91 (MQ=255)
                                                                                          |                                                                                         
TTCCTCTATATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTA  >  USA300TCH1516_ALE/1516355‑1516534

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: