Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,117,393 C→G 6.9% intergenic (+106/+48) suhB_1 → / ← USA300TCH1516_01040 Inositol‑1‑monophosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,117,3930CG 6.9% 406.0 / 15.9 116intergenic (+106/+48)suhB_1/USA300TCH1516_01040Inositol‑1‑monophosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/67);  new base G (2/6);  total (43/73)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 136 total aligned reads displayed.

TATGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTGTAGTAA  >  USA300TCH1516_ALE/1117303‑1117481
                                                                                          |                                                                                        
TATGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTtcatc                                                                                          <  1:3254376‑M1/91‑1 (MQ=255)
TATGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTtcatc                                                                                          <  1:3470220‑M1/91‑1 (MQ=255)
TATGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAACGTtcatc                                                                                          <  1:2070733‑M1/91‑1 (MQ=255)
  tgttacctaattgttaggactacatattacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTcatcat                                                                                        >  1:3558223‑M2/37‑91 (MQ=255)
  TGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTcatcat                                                                                        >  1:2519496‑M1/1‑91 (MQ=255)
   GTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATg                                                                                       >  1:69175‑M1/1‑91 (MQ=255)
     aGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGaa                                                                                     <  1:2680417/91‑1 (MQ=255)
     aGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGaa                                                                                     <  1:2019125/91‑1 (MQ=255)
     aGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGaa                                                                                     <  1:2096266/91‑1 (MQ=255)
     aGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGaa                                                                                     <  1:1483770/91‑1 (MQ=255)
     aGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGaa                                                                                     <  1:2791869/91‑1 (MQ=255)
      gTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAAc                                                                                    <  1:759560/91‑1 (MQ=255)
      acctaattgttaggactacatattacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAAc                                                                                    >  1:3939427‑M2/33‑91 (MQ=255)
        ccTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAt                                                                                  >  1:1243152/1‑91 (MQ=255)
         cTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtg                                                                                 <  1:2842471/91‑1 (MQ=255)
          tAACAGTTAGGTAACATATTACGCTTTTATACCTTAAATGATAAAAACACGCGGCAATATAAAACCAGTTAAATGTTCATCATGAACAtgt                                                                                <  1:3856578/91‑1 (MQ=255)
           aaCAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtg                                                                               <  1:3231897/91‑1 (MQ=255)
           aaCAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtg                                                                               <  1:3232519/91‑1 (MQ=255)
            aCAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATGATGAACAtttaa                                                                              <  1:162538/91‑5 (MQ=255)
            aCAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtgt                                                                              <  1:2791160/91‑1 (MQ=255)
            aCAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtgt                                                                              <  1:259507/91‑1 (MQ=255)
             cAATTAGGTAACATATTACGCTTTTATACCTTGAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtgtg                                                                             <  1:3537544/91‑1 (MQ=255)
             cAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtgtg                                                                             >  1:3828732/1‑91 (MQ=255)
             cAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACAtgtgtg                                                                             >  1:553730/1‑91 (MQ=255)
                 tAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCt                                                                         <  1:394833/91‑1 (MQ=255)
                     tAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTaaa                                                                     >  1:207626/1‑91 (MQ=255)
                     tAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTaaa                                                                     <  1:914899/91‑1 (MQ=255)
                     tAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTaaa                                                                     <  1:3564524/91‑1 (MQ=255)
                     tAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTaaa                                                                     <  1:3201672/91‑1 (MQ=255)
                       acatattacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATc                                                                   >  1:240475‑M2/16‑91 (MQ=255)
                        catattacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCa                                                                  <  1:3899963‑M2/77‑1 (MQ=255)
                          tattacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAc                                                                <  1:1055597‑M2/79‑1 (MQ=255)
                           attacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcg                                                               <  1:3649285‑M2/80‑1 (MQ=255)
                            ttacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgc                                                              <  1:1710098‑M2/81‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:3928133‑M2/82‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:3927940‑M2/82‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:2439682‑M2/82‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:2880712‑M2/82‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:1974980‑M2/82‑1 (MQ=255)
                             tacgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAAcgcg                                                             <  1:3606475‑M2/82‑1 (MQ=255)
                              acgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCgt                                                            >  1:2903676‑M2/9‑91 (MQ=255)
                               cgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCgtg                                                           >  1:3643939‑M2/8‑91 (MQ=255)
                               cgcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCgtg                                                           >  1:3643983‑M2/8‑91 (MQ=255)
                                gcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCCGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGcgtgg                                                          >  1:2314986‑M2/7‑90 (MQ=255)
                                gcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCgtgt                                                          >  1:1054742‑M2/7‑91 (MQ=255)
                                gcttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCgtgt                                                          <  1:346961‑M2/85‑1 (MQ=255)
                                 cttttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGtt                                                         >  1:3377872‑M2/6‑91 (MQ=255)
                                   tttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTAt                                                       <  1:1629583‑M2/88‑1 (MQ=255)
                                   tttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTAt                                                       <  1:1629262‑M2/88‑1 (MQ=255)
                                    ttATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTAtt                                                      <  1:1465716‑M2/89‑1 (MQ=255)
                                       taCCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTc                                                   <  1:1770974/91‑1 (MQ=255)
                                         ccTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGa                                                 >  1:513014/1‑91 (MQ=255)
                                           ttAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGata                                               >  1:1384238/1‑91 (MQ=255)
                                             aaaTAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGAtata                                             >  1:3089350/1‑91 (MQ=255)
                                             aaaTAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGAtata                                             <  1:2979778/91‑1 (MQ=255)
                                                  ataaAAACACGCGTCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTta                                        <  1:4145334/91‑1 (MQ=255)
                                                  ataaAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTta                                        >  1:3392367/1‑91 (MQ=255)
                                                  ataaAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTta                                        <  1:3585167/91‑1 (MQ=255)
                                                  ataaAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTta                                        <  1:1935672/91‑1 (MQ=255)
                                                  ataaAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTta                                        >  1:239120/1‑91 (MQ=255)
                                                   taaAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGCCTTAAATCAACGCGTGTTATTGTCGATATAACTtat                                       >  1:1477391/1‑91 (MQ=255)
                                                     aaaaCACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAg                                     <  1:173224/91‑1 (MQ=255)
                                                          acGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAAt                                <  1:1038866/91‑1 (MQ=255)
                                                          acGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAAt                                <  1:481326/91‑1 (MQ=255)
                                                            gcgcCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCg                              >  1:2625089/1‑91 (MQ=255)
                                                             cgcCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGt                             >  1:1827824/1‑91 (MQ=255)
                                                             cgcCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGt                             >  1:3141587/1‑91 (MQ=255)
                                                             cgcCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGt                             >  1:3141586/1‑91 (MQ=255)
                                                                 aaTATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTtc                         <  1:839605/91‑1 (MQ=255)
                                                                 aaTATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTtc                         <  1:3057825/91‑1 (MQ=255)
                                                                 aaTATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTtc                         >  1:451384/1‑91 (MQ=255)
                                                                  atatAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTtct                        <  1:2531895/91‑1 (MQ=255)
                                                                   tataAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTtctc                       >  1:934762/1‑91 (MQ=255)
                                                                    ataAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCg                      <  1:1211379/91‑1 (MQ=255)
                                                                     taAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGa                     <  1:3747724/91‑1 (MQ=255)
                                                                     taAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGa                     <  1:2502920/91‑1 (MQ=255)
                                                                       aaaCCAGTTAAATGTTCATCATGAACATGTGTGGGCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGata                   <  1:3012262/91‑1 (MQ=255)
                                                                        aaCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGatat                  >  1:1944983/1‑91 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:231791/91‑1 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:1878649/91‑1 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:63875/91‑1 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:231213/91‑1 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:3959440/91‑1 (MQ=255)
                                                                         aCCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtt                 <  1:2776639/91‑1 (MQ=255)
                                                                           cAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAtttt               <  1:1740955/91‑1 (MQ=255)
                                                                              ttAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTc            >  1:3287666/1‑91 (MQ=255)
                                                                               tAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATAttttttat           >  1:4061114/1‑89 (MQ=255)
                                                                               tAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCt           >  1:1970275/1‑91 (MQ=255)
                                                                               tAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCt           >  1:813025/1‑91 (MQ=255)
                                                                                 aaTGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCttt         >  1:3259928/1‑91 (MQ=255)
                                                                                  aTGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTg        <  1:1865624/91‑1 (MQ=255)
                                                                                  aTGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGAGATTTTTTCTTTg        <  1:3175905/91‑1 (MQ=255)
                                                                                      tcatgaTGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagt    >  1:2449626/1‑91 (MQ=255)
                                                                                      tcatgaTGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagt    >  1:2544822/1‑91 (MQ=255)
                                                                                      tcatgaTGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagt    >  1:834009/1‑91 (MQ=255)
                                                                                      tcatgaTGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagt    >  1:4098027/1‑91 (MQ=255)
                                                                                      tcatcaTGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagt    <  1:1716929/91‑1 (MQ=255)
                                                                                       catcatGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTgtagta   >  1:2839057/1‑91 (MQ=255)
                                                                                        atcatGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTGTAGTaa  <  1:3948075/91‑1 (MQ=255)
                                                                                        atcatGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTGTAGTaa  >  1:4127962/1‑91 (MQ=255)
                                                                                          |                                                                                        
TATGTAGTCCTAACAATTAGGTAACATATTACGCTTTTATACCTTAAATAATAAAAACACGCGCCAATATAAAACCAGTTAAATGTTCATCATGAACATGTGTGGTCTTAAATCAACGCGTGTTATTGTCGATATAACTTATAGCCAATCGTTTTCTCGATATTTTTTCTTTGTAGTAA  >  USA300TCH1516_ALE/1117303‑1117481

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: