Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 462,826 T→A 7.4% intergenic (+85/+18) USA300HOU_RS02185 → / ← USA300HOU_RS02190 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE462,8260TA 7.4% 731.7 / 34.4 216intergenic (+85/+18)USA300HOU_RS02185/USA300HOU_RS02190hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (95/105);  new base A (10/6);  total (105/111)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 220 total aligned reads displayed.

AATAATTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAA  >  USA300TCH1516_ALE/462737‑462916
                                                                                         |                                                                                          
aataatTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTa                                                                                           <  1:446449/91‑1 (MQ=255)
  taatTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTAACa                                                                                         <  2:1598072/91‑1 (MQ=255)
  taatTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTAACa                                                                                         <  1:2254046/91‑1 (MQ=255)
    atTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACaaa                                                                                       >  1:1223367/1‑91 (MQ=255)
    atTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACaaa                                                                                       >  1:3358787/1‑91 (MQ=255)
    atTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACaaa                                                                                       >  1:3358627/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTTCTTTGTTACAAAGc                                                                                     >  2:531678/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  2:2634513/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  2:1345803/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:2056653/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  2:3947590/1‑91 (MQ=255)
      taaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATCGAGCTGTTTTTTGTGCTTTGTTACAAAGc                                                                                     >  1:2630149/1‑91 (MQ=255)
       aaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCa                                                                                    <  1:2624707/91‑1 (MQ=255)
       aaataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCa                                                                                    >  2:2883289/1‑91 (MQ=255)
         ataaaACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCatt                                                                                  <  1:304137/91‑1 (MQ=255)
           aaaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattat                                                                                <  1:2827926/91‑1 (MQ=255)
           aaaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattat                                                                                <  2:1582935/91‑1 (MQ=255)
           aaaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattat                                                                                <  1:2652289/91‑1 (MQ=255)
            aaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattatt                                                                               <  1:2933375/91‑1 (MQ=255)
            aaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattatt                                                                               <  1:2175356/91‑1 (MQ=255)
            aaaCAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCattatt                                                                               <  1:2175486/91‑1 (MQ=255)
                aGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAat                                                                           <  2:3893764/91‑1 (MQ=255)
                aGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAat                                                                           <  1:2485959/91‑1 (MQ=255)
                 gCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatt                                                                          >  2:2539089/1‑91 (MQ=255)
                 gCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatt                                                                          >  2:719942/1‑91 (MQ=255)
                  cTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAattt                                                                         >  2:2182977/1‑91 (MQ=255)
                    tAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttat                                                                       <  2:195702/91‑1 (MQ=255)
                     aaTGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttatt                                                                      >  1:144277/1‑91 (MQ=255)
                     aaTGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttatt                                                                      <  1:627720/91‑1 (MQ=255)
                      atgatgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttattt                                                                     <  1:1893680/91‑1 (MQ=255)
                      atgatgTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAatttattt                                                                     <  1:797504/91‑1 (MQ=255)
                           gTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACgt                                                                >  1:377369/1‑91 (MQ=255)
                           gTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACgt                                                                >  1:3562239/1‑91 (MQ=255)
                                  aaTTAGTGATTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCata                                                         >  1:1980260/1‑91 (MQ=255)
                                  aaTTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCata                                                         >  1:1994133/1‑91 (MQ=255)
                                   aTTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTCTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCatat                                                        >  2:3051493/1‑91 (MQ=255)
                                      aGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTTTGCTTTTTTTCAAAGCATTATTGAATTTATTTTACGTGTTCATAtttt                                                     >  2:671852/1‑91 (MQ=255)
                                      aGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATAtttt                                                     >  2:2352363/1‑91 (MQ=255)
                                       gTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAATCATTATTGAATTTATTTTACGTGTTCATATTTTg                                                    >  2:1411002/1‑91 (MQ=255)
                                       gTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTg                                                    >  1:1346050/1‑91 (MQ=255)
                                       gTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTg                                                    >  1:2654025/1‑91 (MQ=255)
                                        tGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGa                                                   <  1:1181773/91‑1 (MQ=255)
                                        tGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGa                                                   <  1:1128977/91‑1 (MQ=255)
                                           gTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAAc                                                <  2:3727037/91‑1 (MQ=255)
                                             tAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTATTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACAt                                              >  2:1246565/1‑91 (MQ=255)
                                                tCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCaa                                           <  2:122851/91‑1 (MQ=255)
                                                tCACTGACTTCTACGTCATTGAGCTGTTTTTTGGGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCaa                                           <  1:3504366/91‑1 (MQ=255)
                                                 cACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCaaa                                          >  1:2966210/1‑91 (MQ=255)
                                                  aCTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAg                                         <  2:2056653/91‑1 (MQ=255)
                                                  aCTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTGGAAACATCAAAg                                         <  1:3027458/91‑1 (MQ=255)
                                                   cTGACTTCTGCGTTATTGAGCTGGTTTTGGTGCTTTGTTACAAAGCAGTATTGAATTGTTTCTAGGTGTTCACATTCTGAAACAGCAAAGc                                        <  2:3352354/91‑1 (MQ=255)
                                                   cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  2:1726548/91‑1 (MQ=255)
                                                   cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  1:727578/91‑1 (MQ=255)
                                                   cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  2:2014643/91‑1 (MQ=255)
                                                   cTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGc                                        <  1:996985/91‑1 (MQ=255)
                                                     gACTTCTACGTCATTGAGCTGTTTTTTTTTCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCg                                      >  2:322786/1‑91 (MQ=255)
                                                     gACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCg                                      >  2:3444782/1‑91 (MQ=255)
                                                        ttCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCt                                   >  1:3188570/1‑91 (MQ=255)
                                                          cTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTg                                 <  1:3286400/91‑1 (MQ=255)
                                                              gTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTa                             >  1:348856/1‑91 (MQ=255)
                                                               tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            >  1:50839/1‑91 (MQ=255)
                                                               tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            <  1:1353260/91‑1 (MQ=255)
                                                               tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            >  1:3887227/1‑91 (MQ=255)
                                                               tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            <  1:3580797/91‑1 (MQ=255)
                                                               tCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAg                            >  1:1818421/1‑91 (MQ=255)
                                                                cATTAAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGc                           <  2:3299895/91‑1 (MQ=255)
                                                                  ttGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTCGAAACATCAAAGCCGTCTTGTTTAGCtt                         >  1:903318/1‑91 (MQ=255)
                                                                   tGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCttt                        <  2:1885233/91‑1 (MQ=255)
                                                                   tGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCttt                        <  2:3796723/91‑1 (MQ=255)
                                                                    gAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTttg                       <  1:2749957/91‑1 (MQ=255)
                                                                     aGCTGTTTTTTGTGCTTTTTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTCGCTttgt                      >  2:3847731/1‑91 (MQ=255)
                                                                      gCTGTTTTTTTTTCTTTTTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTttgtt                     >  2:3129135/1‑91 (MQ=255)
                                                                         gTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGAt                  <  2:3867940/91‑1 (MQ=255)
                                                                          ttttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATa                 >  1:204935/1‑91 (MQ=255)
                                                                           tttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATaa                <  2:863717/91‑1 (MQ=255)
                                                                           tttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATaa                <  1:3731666/91‑1 (MQ=255)
                                                                           tttttGTGCTTTGTTACAAAGCATTATTGAATTTATTTGACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATaa                <  2:3845764/91‑1 (MQ=255)
                                                                             tttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATg              <  1:2328651/91‑1 (MQ=255)
                                                                             tttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATg              <  1:2526901/91‑1 (MQ=255)
                                                                             tttGTGCTTTGTTACAAAGCATGATTGAATTTATTTTACGTGTTCATATCTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATg              <  2:3007059/91‑1 (MQ=255)
                                                                              ttGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGt             >  1:1281822/1‑91 (MQ=255)
                                                                               tgtgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTc            >  1:2672186/1‑91 (MQ=255)
                                                                                gtgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCt           <  1:1928098/91‑1 (MQ=255)
                                                                                gtgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCt           <  1:2597903/91‑1 (MQ=255)
                                                                                gtgCTTTGTTACAAAGCATGCTTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCt           >  1:3531442/1‑91 (MQ=255)
                                                                                 tgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCtt          >  2:503265/1‑91 (MQ=255)
                                                                                 tgCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCtt          <  2:3033213/91‑1 (MQ=255)
                                                                                 tgCTTTGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTTTTGTTTAGCTTTGTTGATAATGTCtt          >  1:1242382/1‑91 (MQ=255)
                                                                                  gCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCttt         <  1:3655003/91‑1 (MQ=255)
                                                                                  gCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCttt         >  2:2441280/1‑91 (MQ=255)
                                                                                   cTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCAGATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttg        >  1:1029428/1‑91 (MQ=255)
                                                                                    tttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttga       >  1:1081305/1‑91 (MQ=255)
                                                                                    tttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttga       >  1:921016/1‑91 (MQ=255)
                                                                                    tttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttga       >  2:1355882/1‑91 (MQ=255)
                                                                                     ttGTAACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgat      >  1:61519/1‑91 (MQ=255)
                                                                                      tGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgatt     <  2:1609861/91‑1 (MQ=255)
                                                                                       gTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTttgattg    >  1:3779113/1‑91 (MQ=255)
                                                                                         tACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGaa  <  2:3531506/91‑1 (MQ=255)
                                                                                         tACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGaa  <  1:57045/91‑1 (MQ=255)
                                                                                         tACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGaa  <  1:648949/91‑1 (MQ=255)
                                                                                         |                                                                                          
AATAATTAAATAAAACAGCTTAATGATGTAATGAAATTAGTGAGTTAATCACTGACTTCTACGTCATTGAGCTGTTTTTTGTGCTTTGTTACAAAGCATTATTGAATTTATTTTACGTGTTCATATTTTGAAACATCAAAGCCGTCTTGTTTAGCTTTGTTGATAATGTCTTTGATTGAA  >  USA300TCH1516_ALE/462737‑462916

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