Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 690,098 A→T 10.8% intergenic (‑27/+8) USA300HOU_RS03320 ← / ← USA300HOU_RS03325 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE690,0980AT10.8% 863.9 / 69.3 282intergenic (‑27/+8)USA300HOU_RS03320/USA300HOU_RS03325hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (123/128);  new base T (11/20);  total (134/148)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 288 total aligned reads displayed.

CTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAA  >  USA300TCH1516_ALE/690008‑690188
                                                                                          |                                                                                          
  tCGACAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCTCAAt                                                                                          <  1:3122353/91‑1 (MQ=255)
   cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                         >  2:3524391/1‑91 (MQ=255)
   cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAATAACCCCAATAg                                                                                        >  2:2858391/1‑91 (MQ=255)
    gCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAg                                                                                        <  2:2330602/91‑1 (MQ=255)
     ccAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAgg                                                                                       >  1:3428995/1‑91 (MQ=255)
      cAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAggg                                                                                      <  2:563460/91‑1 (MQ=255)
      cAGTTTTCGTTAACAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAggg                                                                                      >  1:469067/1‑91 (MQ=255)
       aGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAgggg                                                                                     <  1:3855885/91‑1 (MQ=255)
       aGTTTTCGTTAATAACAAATGATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAgggg                                                                                     <  2:3246617/91‑1 (MQ=255)
        gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCGTACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                    <  2:622849/91‑1 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                   <  1:3500223/91‑3 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                   <  1:41007/91‑3 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                   <  1:2781796/91‑3 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCt                                                                                   >  1:72198/1‑91 (MQ=255)
         ttttCGTCAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                   <  1:1974929/91‑3 (MQ=255)
          tttCGTTAAGAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                  <  1:3795205/91‑4 (MQ=255)
             cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                               >  2:3002815/1‑91 (MQ=255)
               ttAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTaa                                                                             >  2:1570419/1‑91 (MQ=255)
               ttAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTaa                                                                             >  2:1916993/1‑91 (MQ=255)
                 aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                           <  2:1617992/91‑1 (MQ=255)
                  ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                          >  2:226036/1‑91 (MQ=255)
                  ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                          <  1:163103/91‑1 (MQ=255)
                   taaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGt                                                                         >  2:637359/1‑91 (MQ=255)
                    aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                        <  1:2733254/91‑1 (MQ=255)
                     aCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTa                                                                       <  2:3855705/91‑1 (MQ=255)
                       aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                     <  2:2003757/91‑1 (MQ=255)
                        aaTTATCTCTTAGACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttt                                                                    >  1:3036462/1‑91 (MQ=255)
                        aaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttt                                                                    <  2:1258834/91‑1 (MQ=255)
                         aTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttt                                                                   <  1:2188858/91‑1 (MQ=255)
                             tctcTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTa                                                               >  1:1016683/1‑91 (MQ=255)
                              ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                              <  1:950700/91‑1 (MQ=255)
                              ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                              <  2:1386291/91‑1 (MQ=255)
                              ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                              <  1:114743/91‑1 (MQ=255)
                               tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                             <  1:3930449/91‑1 (MQ=255)
                                 ttAAACCTGCATACGGTGCTTTTAAATTCATGTCATTCCTCCAGTAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATcgt                                                           <  2:2446301/91‑3 (MQ=255)
                                 ttAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCAt                                                           >  2:1831387/1‑91 (MQ=255)
                                  tAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                          >  2:3493027/1‑91 (MQ=255)
                                  tAAACCTGCATACGGTGATTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                          <  1:3389271/91‑1 (MQ=255)
                                   aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                         >  1:1572194/1‑91 (MQ=255)
                                      ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                      <  1:1479183/91‑1 (MQ=255)
                                      ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                      <  1:810225/91‑1 (MQ=255)
                                      ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                      <  2:3704544/91‑1 (MQ=255)
                                       cTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATAtt                                                     <  1:293840/91‑1 (MQ=255)
                                        tGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTc                                                    >  1:565041/1‑91 (MQ=255)
                                        tGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTc                                                    >  1:2171662/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                   >  2:1077707/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                   >  1:1321591/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                   >  1:1900173/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                   >  1:268227/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                   >  1:2638324/1‑91 (MQ=255)
                                          cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                  <  2:2620501/91‑1 (MQ=255)
                                          cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                  >  2:380048/1‑91 (MQ=255)
                                           aTACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAc                                                 >  1:3012633/1‑91 (MQ=255)
                                            tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                >  1:3339405/1‑91 (MQ=255)
                                              cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                              <  2:2535915/91‑1 (MQ=255)
                                                 tGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTGTTTATCATGATATTCTACttcttc                                           <  2:3361706/91‑1 (MQ=255)
                                                  gCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtca                                          >  1:3652624/1‑91 (MQ=255)
                                                  gCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtca                                          >  2:3879841/1‑91 (MQ=255)
                                                   cTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtcat                                         <  2:1363572/91‑1 (MQ=255)
                                                    ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtcatc                                        <  1:993816/91‑1 (MQ=255)
                                                      ttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTcatcat                                      >  1:2053666/1‑91 (MQ=255)
                                                       tAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATCCTACTTCTTCATCATa                                     >  2:771505/1‑91 (MQ=255)
                                                       tAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCGATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATa                                     <  1:3138944/91‑1 (MQ=255)
                                                        aaaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATaa                                    <  2:403610/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                   <  1:921566/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                   >  2:3068015/1‑91 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                   >  2:2605469/1‑91 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                   <  2:3222435/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                   <  1:1751850/91‑1 (MQ=255)
                                                          aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                  <  1:2910849/91‑1 (MQ=255)
                                                             cATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTaa                               >  2:2761205/1‑91 (MQ=255)
                                                              aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                              <  1:764027/91‑1 (MQ=255)
                                                              aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACGTAAt                              >  1:3285843/1‑91 (MQ=255)
                                                                  cATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTaa                          >  1:3499199/1‑91 (MQ=255)
                                                                   aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                         >  1:1030914/1‑91 (MQ=255)
                                                                   aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                         >  2:678080/1‑91 (MQ=255)
                                                                     cccTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtt                       <  2:2195245/91‑1 (MQ=255)
                                                                     cccTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtt                       <  2:2084000/91‑1 (MQ=255)
                                                                       cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                     <  2:1745539/91‑1 (MQ=255)
                                                                        tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                    >  1:2548057/1‑91 (MQ=255)
                                                                        tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                    >  2:3316027/1‑91 (MQ=255)
                                                                            atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTCTTCTTaa                >  2:232773/5‑91 (MQ=37)
                                                                            aTAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                <  1:1376515/91‑1 (MQ=255)
                                                                             tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc               <  1:3764952/91‑1 (MQ=255)
                                                                               aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg             <  2:270991/91‑1 (MQ=255)
                                                                               aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg             >  1:2475489/1‑91 (MQ=255)
                                                                               aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg             <  1:2881457/91‑1 (MQ=255)
                                                                               aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg             >  1:520011/1‑91 (MQ=255)
                                                                                aaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGc            >  2:1789461/1‑91 (MQ=38)
                                                                                aaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGc            >  2:3160829/1‑91 (MQ=39)
                                                                                  aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcac          >  2:1481910/1‑91 (MQ=255)
                                                                                    agccccTATTGGGGTTATTTAAATGTTATTTTTTATCATGATATTCTACTTTTTCATTCTAACTTAATTTTAACTTTCTTCACTGCACAtt       >  2:2595483/3‑91 (MQ=255)
                                                                                    aaCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt        >  1:325249/1‑91 (MQ=255)
                                                                                     gccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt       >  1:2092180/2‑91 (MQ=255)
                                                                                     gccccTATTGGGGTTATTTAAATGTTATTTTTTTATAATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt       >  1:1294516/2‑91 (MQ=255)
                                                                                      ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa      >  1:1413441/1‑91 (MQ=255)
                                                                                        ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa    <  1:362681/91‑1 (MQ=255)
                                                                                         cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa   <  1:3358715/91‑1 (MQ=255)
                                                                                          aaTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa  <  1:168126/91‑1 (MQ=255)
                                                                                          aaTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa  >  2:1045223/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAA  >  USA300TCH1516_ALE/690008‑690188

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