Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,032,588 G→T 5.1% intergenic (+61/+118) yfkN_2 → / ← USA300TCH1516_00957 Trifunctional nucleotide phosphoesterase protein YfkN/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,032,5880GT 5.1% 806.0 / 20.1 232intergenic (+61/+118)yfkN_2/USA300TCH1516_00957Trifunctional nucleotide phosphoesterase protein YfkN/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (105/115);  new base T (8/4);  total (113/119)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 252 total aligned reads displayed.

GTGATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACAAAAAAGCT  >  USA300TCH1516_ALE/1032499‑1032677
                                                                                         |                                                                                         
gTGATAATTTTATAGTTAAACACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGagaaga                                                                                          >  1:88480/1‑91 (MQ=255)
gTGATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGagaaga                                                                                          <  2:695475/91‑1 (MQ=255)
gTGATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGagaaga                                                                                          <  2:2705154/91‑1 (MQ=255)
 tGATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGAt                                                                                         >  1:3675062/1‑91 (MQ=255)
  gATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGAtt                                                                                        <  1:1771725/91‑1 (MQ=255)
    tAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTAt                                                                                      <  1:3562984/91‑1 (MQ=255)
      aTTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCt                                                                                    <  1:2699490/91‑1 (MQ=255)
       ttttATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAACAACCTTTCTCAGAGAATATTATCTa                                                                                   >  2:3365041/1‑91 (MQ=255)
       ttttATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTa                                                                                   <  1:133297/91‑1 (MQ=255)
        tttATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAg                                                                                  <  2:1818363/91‑1 (MQ=255)
        tttATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAg                                                                                  >  1:1285440/1‑91 (MQ=255)
         ttATAGTTACAAACAGGTAAATAGTGCTATGTAATAAATGAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAgg                                                                                 <  2:403173/91‑1 (MQ=255)
         ttATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAgg                                                                                 <  2:465750/91‑1 (MQ=255)
         ttAGAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATGATCTAgg                                                                                 <  2:402407/91‑1 (MQ=255)
          tataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGa                                                                                <  2:3428488/91‑1 (MQ=255)
          tataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGa                                                                                >  2:676463/1‑91 (MQ=255)
          tataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGa                                                                                <  2:772913/91‑1 (MQ=255)
           ataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGaa                                                                               >  2:1760732/1‑91 (MQ=255)
           ataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGaa                                                                               <  2:3145300/91‑1 (MQ=255)
           ataGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGaa                                                                               <  1:1745384/91‑1 (MQ=255)
            taGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGaaa                                                                              >  1:1449544/1‑91 (MQ=255)
             aGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAg                                                                             >  1:2004092/1‑91 (MQ=255)
              gTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAgg                                                                            <  1:2775384/91‑1 (MQ=255)
              gTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAgg                                                                            >  1:3229002/1‑91 (MQ=255)
              gTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAgg                                                                            >  1:2418733/1‑91 (MQ=255)
               ttAAAAACATGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGt                                                                           <  2:2329576/91‑1 (MQ=255)
               ttAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCCTAGATAATCTTCTCTAAGAAAGGt                                                                           <  1:31638/91‑1 (MQ=21)
                tAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGtt                                                                          <  1:2680973/91‑1 (MQ=255)
                 aaaaaCAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGttt                                                                         <  2:3857515/91‑1 (MQ=255)
                 aaaaaCAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCCTAGATAATCTTCTCTAAGAAAGGttt                                                                         <  1:2866299/91‑1 (MQ=21)
                  aaaaCAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGtttt                                                                        <  2:1710364/91‑1 (MQ=255)
                   aaaCAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGttttt                                                                       <  1:1863893/91‑1 (MQ=255)
                     aCAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAg                                                                     <  2:1923777/91‑1 (MQ=255)
                       aGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGaa                                                                   <  1:428934/91‑1 (MQ=255)
                       aGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGaa                                                                   <  1:1834303/91‑1 (MQ=255)
                        ggTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAAt                                                                  >  1:773056/1‑91 (MQ=255)
                          tAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAAttt                                                                <  1:461588/91‑1 (MQ=255)
                              tatTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTaaa                                                            >  2:124640/1‑91 (MQ=255)
                                 tGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAAtt                                                         >  1:1708722/1‑91 (MQ=255)
                                 tGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAAtt                                                         <  2:2224594/91‑1 (MQ=255)
                                 tGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAAtt                                                         <  1:1577394/91‑1 (MQ=255)
                                 tGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAAtt                                                         <  2:2863598/91‑1 (MQ=255)
                                 tGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAAtt                                                         <  2:3557612/91‑1 (MQ=255)
                                  gCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTg                                                        >  2:260310/1‑91 (MQ=255)
                                  gCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTg                                                        >  2:3004984/1‑91 (MQ=255)
                                     aTGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGata                                                     >  2:3079849/1‑91 (MQ=255)
                                      tGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGatat                                                    >  2:2404623/1‑91 (MQ=255)
                                       gTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGAtata                                                   <  1:3580688/91‑1 (MQ=255)
                                        taataaATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATaa                                                  >  1:2484809/1‑91 (MQ=255)
                                        taataaATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATaa                                                  <  1:2963745/91‑1 (MQ=255)
                                         aataaATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATaaa                                                 <  2:2819335/91‑1 (MQ=255)
                                            aaaTTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGa                                              >  2:2513941/1‑91 (MQ=255)
                                             aattaattATAAGTTCAACAAGTATAAAAACCTTTCTTAGTGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGaa                                             >  2:3762726/1‑91 (MQ=255)
                                               ttaattaTAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACa                                           <  2:2679778/91‑1 (MQ=255)
                                               ttaattaTAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACa                                           <  2:504315/91‑1 (MQ=255)
                                                 aattaTAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGt                                         >  2:1463701/1‑91 (MQ=255)
                                                   ttaTAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGa                                       >  1:94334/1‑91 (MQ=255)
                                                    tataGGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGAGaa                                      >  1:3062335/1‑91 (MQ=255)
                                                    tataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaa                                      <  2:3182617/91‑1 (MQ=255)
                                                    tataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaa                                      >  2:97089/1‑91 (MQ=255)
                                                    tataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaa                                      <  1:1402663/91‑1 (MQ=255)
                                                    tataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaa                                      <  1:334353/91‑1 (MQ=255)
                                                    tataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaa                                      <  1:3703157/91‑1 (MQ=255)
                                                     ataAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGaaa                                     <  2:2417465/91‑1 (MQ=255)
                                                      taAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAAt                                    <  1:1778430/91‑1 (MQ=255)
                                                      taAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTAACTAGGAAAGGTTTTTAGAATTTTAATCTTGATATAAAAGAACAGTGAAAt                                    >  2:3608733/1‑91 (MQ=255)
                                                       aaGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAAtt                                   <  2:788376/91‑1 (MQ=255)
                                                       aaGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAAtt                                   <  2:2960414/91‑1 (MQ=255)
                                                       aaGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGATTTTAGAATCTTAAAATTGATATAAAAGAACAGTGAAAtt                                   >  2:419396/1‑91 (MQ=255)
                                                        aGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATtc                                  >  1:3467905/1‑91 (MQ=255)
                                                        aGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATtc                                  >  1:3671452/1‑91 (MQ=255)
                                                          ttCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATtctc                                <  1:1028897/91‑1 (MQ=255)
                                                           tCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtca                               <  1:2188402/91‑1 (MQ=255)
                                                             aacaaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcatt                             <  1:114130/91‑1 (MQ=255)
                                                              acaaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattc                            >  2:227276/1‑91 (MQ=255)
                                                              acaaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattc                            <  1:637011/91‑1 (MQ=255)
                                                              acaaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattc                            <  2:110341/91‑1 (MQ=255)
                                                               caaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattca                           >  2:2462294/1‑91 (MQ=255)
                                                               caaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattca                           >  1:1991785/1‑91 (MQ=255)
                                                                aaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattcat                          <  2:1627288/91‑1 (MQ=255)
                                                                aaGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCtcattcat                          <  2:3675062/91‑1 (MQ=255)
                                                                  gTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATta                        <  2:773056/91‑1 (MQ=255)
                                                                   tataAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATtat                       >  2:1588470/1‑91 (MQ=255)
                                                                    ataAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATtata                      >  1:548842/1‑91 (MQ=255)
                                                                        aaaCCTTTCTTAGAGAAGATTATCTAGGAAAGGTATTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCCCATTCATtatatata                  >  2:2041714/1‑91 (MQ=255)
                                                                         aaCCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATGCATTATATATaa                 >  2:3229259/1‑91 (MQ=255)
                                                                         aaCCTTTCCTAGATAATCTTCTCTAAGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATaa                 >  1:1358864/1‑91 (MQ=21)
                                                                          aCCTTTCCTAGATAATCTTCTCTAAGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAg                >  1:2919104/1‑91 (MQ=25)
                                                                          aCCTTTCCTAGATAATCTTCTCTAAGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAg                >  2:3488863/1‑91 (MQ=25)
                                                                           ccTTTCCTAGATAATCTTCTCTAAGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGt               >  1:1095950/1‑91 (MQ=25)
                                                                             tttCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATCGATATAAAAGAACAGTGAAATTCTCATTCATtatatat                   <  2:1803498/85‑1 (MQ=255)
                                                                                cTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACa          >  2:896007/1‑91 (MQ=255)
                                                                                 ttAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaa         <  1:2231246/91‑1 (MQ=255)
                                                                                  tAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaaa        >  2:3729449/1‑91 (MQ=255)
                                                                                  tAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaaa        >  2:3202211/1‑91 (MQ=255)
                                                                                     agaagaTTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaaaaaa     <  1:1608239/91‑1 (MQ=255)
                                                                                     agaagaTTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaaaaaa     <  1:742442/91‑1 (MQ=255)
                                                                                     agaagaTTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACaaaaaa     <  1:3853354/91‑1 (MQ=255)
                                                                                      gaagaTTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAGGAACAGTGAAATTCTCATTCATTATATATAAGTAAACAAAAAAg    <  2:3321840/91‑1 (MQ=255)
                                                                                       aagaTTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACAAAAAAGc   <  2:1013296/91‑1 (MQ=255)
                                                                                         |                                                                                         
GTGATAATTTTATAGTTAAAAACAGGTAAATATTGCTATGTAATAAATTAATTATAAGTTCAACAAGTATAAAAACCTTTCTTAGAGAAGATTATCTAGGAAAGGTTTTTAGAATTTTAAAATTGATATAAAAGAACAGTGAAATTCTCATTCATTATATATAAGTAAACAAAAAAGCT  >  USA300TCH1516_ALE/1032499‑1032677

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: