Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,909,134 T→A 5.4% intergenic (‑200/‑60) tal ← / → USA300HOU_RS09460 Transaldolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,909,1340TA 5.4% 924.3 / 27.2 270intergenic (‑200/‑60)tal/USA300HOU_RS09460Transaldolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (110/145);  new base A (2/13);  total (112/158)
Fisher's exact test for biased strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

ATATCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTAT  >  USA300TCH1516_ALE/1909044‑1909224
                                                                                          |                                                                                          
atatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACATGTTAATAAAACt                                                                                            <  1:3049756/91‑1 (MQ=255)
atatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACt                                                                                            <  1:3592965/91‑1 (MQ=255)
 tatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAACAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTa                                                                                           >  1:493877/1‑91 (MQ=255)
 tatCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTa                                                                                           <  2:1449665/91‑1 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                       <  2:2511014/91‑1 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                       <  1:820850/91‑1 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                       >  1:219109/1‑91 (MQ=255)
     aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAAtagt                                                                                       <  2:1749468/91‑4 (MQ=255)
     aCGCTTTTATAAGATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                       <  2:3286315/91‑1 (MQ=255)
      cGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtt                                                                                      >  2:2322453/1‑91 (MQ=255)
      cGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtt                                                                                      >  2:1867273/1‑91 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATGAAACAATAGttt                                                                                     <  1:3047716/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                     >  2:3037964/1‑91 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                     <  1:3375505/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                     <  2:1328038/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                     <  1:1721581/91‑1 (MQ=255)
       gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                     <  1:1516672/91‑1 (MQ=255)
        cTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                    <  1:3749599/91‑1 (MQ=255)
        cTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                    <  2:2424607/91‑1 (MQ=255)
         ttttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTa                                                                                   <  2:2126083/91‑1 (MQ=255)
          tttATAATATATAGCTAATTTTTTGCCTTTTTTTTTTAAAAAAGAGTGTATCAAACATTGTGTTTACAATTTAATAAAACTATTTTTTTAt                                                                                  >  2:1349983/1‑91 (MQ=255)
          tttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTAt                                                                                  <  1:1438661/91‑1 (MQ=255)
           ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                 <  1:1381103/91‑1 (MQ=255)
           ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                 >  2:2051813/1‑91 (MQ=255)
           ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                 <  2:1003992/91‑1 (MQ=255)
           ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                 >  1:935912/1‑91 (MQ=255)
            tataATATATAGCTAATTTTTTTCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGa                                                                                >  2:3591647/1‑91 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATTTTTTTTAAAAATTACTTTAGCACACATTGTGTTTACAAGTTAATATTACTATTGTTTTATGAt                                                                               >  2:58861/1‑91 (MQ=255)
             ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                               <  2:887131/91‑1 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAATATGACTTTAGCAAACATTGTGTCTACCAGTTAATAAAACTATTGTTTTTTGAtt                                                                              >  2:2716805/1‑91 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                              >  1:3473552/1‑91 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                              <  2:2317503/91‑1 (MQ=255)
              taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                              <  2:632410/91‑1 (MQ=255)
               aataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTg                                                                             >  1:3427486/1‑91 (MQ=255)
                 tatataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGtt                                                                           <  2:109662/91‑1 (MQ=255)
                 tatataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGtt                                                                           <  1:1674756/91‑1 (MQ=255)
                  atataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGttt                                                                          <  1:460684/91‑1 (MQ=255)
                  atataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGttt                                                                          <  1:202613/91‑1 (MQ=255)
                   tataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTa                                                                         <  1:2348245/91‑1 (MQ=255)
                      aGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGa                                                                      <  2:289598/91‑1 (MQ=255)
                        cTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATa                                                                    >  1:2063354/1‑91 (MQ=255)
                         tAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATAc                                                                   >  2:2790340/1‑91 (MQ=255)
                           aTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGa                                                                 <  1:3888495/91‑1 (MQ=255)
                              ttttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGatat                                                              >  2:48567/1‑91 (MQ=255)
                              ttatGCCTCTTTTTTGTAAAACGTAGTTTAGAGAACATTTTGTATACAAGTTAAAACAACTATTGTTGTAGGATTGTTTATGATAAGatat                                                              <  2:2573218/88‑1 (MQ=255)
                               tttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATg                                                             <  1:2950455/91‑1 (MQ=255)
                                ttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGt                                                            <  2:2667263/91‑1 (MQ=255)
                                ttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGt                                                            >  2:1314399/1‑91 (MQ=255)
                                ttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGCCTACAAGTTAATAAAACTATTGTTTGATGATTGTTTATGATACGATATGt                                                            >  2:3667519/1‑91 (MQ=255)
                                  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaa                                                          <  2:3217903/91‑1 (MQ=255)
                                  gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaa                                                          >  1:2216735/1‑91 (MQ=255)
                                   ccTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaaa                                                         <  1:741467/91‑1 (MQ=255)
                                   ccTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaaa                                                         >  2:1585352/1‑91 (MQ=255)
                                    cTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAAt                                                        >  2:2695030/1‑91 (MQ=255)
                                     tataTTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATg                                                       <  2:3147924/91‑1 (MQ=255)
                                       tatTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGAt                                                     >  2:2979049/1‑91 (MQ=255)
                                       tatTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGAt                                                     <  1:2531606/91‑1 (MQ=255)
                                         ttttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATaa                                                   <  1:3277326/91‑1 (MQ=255)
                                         ttttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATaa                                                   >  1:129259/1‑91 (MQ=255)
                                           ttttAAAAATGACTTTAGCAAACTTTGTGTCTACCAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatc                                                 >  2:1939096/1‑91 (MQ=255)
                                           ttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatc                                                 <  2:1534544/91‑1 (MQ=255)
                                             ttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcat                                               >  1:1087188/1‑91 (MQ=255)
                                             ttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAAATATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcat                                               <  1:3134866/91‑1 (MQ=255)
                                              tAAAATTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                              >  2:1508130/1‑91 (MQ=255)
                                              tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                              <  2:3376607/91‑1 (MQ=255)
                                               aaaaaTGACTTTAGCAGACATGGTGGCGACAAGTTAATAAAACTATTGTTGTATGATTGTTTATGAGACGATATGTAAATGATAATCATcc                                             <  2:1199429/91‑1 (MQ=255)
                                                aaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCa                                            <  2:3213733/91‑1 (MQ=255)
                                                aaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGAGAATCATCCa                                            <  1:1969784/91‑1 (MQ=255)
                                                  aaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCagg                                          >  1:90048/1‑91 (MQ=255)
                                                    tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                        <  2:1210102/91‑1 (MQ=255)
                                                    tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                        <  1:1002621/91‑1 (MQ=255)
                                                    tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                        <  2:39144/91‑1 (MQ=255)
                                                     gACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggagg                                       <  2:1407720/91‑1 (MQ=255)
                                                      aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                      <  1:3643949/91‑1 (MQ=255)
                                                      aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                      <  1:3643921/91‑1 (MQ=255)
                                                      aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                      <  2:3089219/91‑1 (MQ=255)
                                                         ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGATTATGATACGATATGTAAATGATAATCATCTAGGAGGTCTa                                   >  2:3469431/1‑91 (MQ=255)
                                                          tAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAg                                  >  2:2628619/1‑91 (MQ=255)
                                                           aGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGc                                 <  1:2483444/91‑1 (MQ=255)
                                                           aGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGc                                 >  1:2841824/1‑91 (MQ=255)
                                                            gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                                <  1:3657466/91‑1 (MQ=255)
                                                            gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                                >  2:1027798/1‑91 (MQ=255)
                                                              aaaCATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGACTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATg                              >  1:3420125/1‑91 (MQ=255)
                                                                  aTTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtctc                          <  1:1944821/91‑1 (MQ=255)
                                                                    tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                        <  1:1844336/91‑1 (MQ=255)
                                                                     gtgtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGtt                       >  2:1649395/1‑91 (MQ=255)
                                                                      tgtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTc                      >  1:1860187/1‑91 (MQ=255)
                                                                       gtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCa                     >  2:94591/1‑91 (MQ=255)
                                                                       gtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCa                     >  2:1990150/1‑91 (MQ=255)
                                                                        tCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaa                    >  1:163393/1‑91 (MQ=255)
                                                                          tACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaa                  <  1:54272/91‑1 (MQ=255)
                                                                           aCAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaa                 >  2:2358346/1‑91 (MQ=255)
                                                                            cAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaaa                <  1:2603128/91‑1 (MQ=255)
                                                                                    taaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt        <  1:189036/91‑1 (MQ=255)
                                                                                    taaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt        >  2:695355/1‑91 (MQ=255)
                                                                                    taaAACAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATACAGGAGGTCTAGCATGTCTCGTTCAAAACAATACttct        >  2:3590880/1‑89 (MQ=255)
                                                                                        aCTATTGTTTTATTTTTTTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACtt    >  2:1056975/1‑91 (MQ=255)
                                                                                        aCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACtt    >  1:881217/1‑91 (MQ=255)
                                                                                        aCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACtt    >  1:1492465/1‑91 (MQ=255)
                                                                                          tATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTAt  <  1:3304863/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATATCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTAT  >  USA300TCH1516_ALE/1909044‑1909224

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