Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 895,942 C→T 6.9% intergenic (+141/+2) metQ_2 → / ← Int‑Tn_1 Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE895,9420CT 6.9% 429.5 / 14.4 131intergenic (+141/+2)metQ_2/Int‑Tn_1Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916
Reads supporting (aligned to +/- strand):  ref base C (44/77);  new base T (6/3);  total (51/80)
Fisher's exact test for biased strand distribution p-value = 8.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 131 total aligned reads displayed.

TGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTA  >  USA300TCH1516_ALE/895853‑896028
                                                                                         |                                                                                      
 gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                      <  1:3089410/91‑1 (MQ=255)
  tgGTATTGAAAATAAAAAAGGGCACAAANGGGGNAGTATATGCAANTTATCGANTTTTAANCAAAAAAATGCGCTCAAATCTGAGCGTAtt                                                                                     <  1:1174561/91‑1 (MQ=37)
  tgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAGCGTAtt                                                                                     <  1:2704912/91‑1 (MQ=255)
   ggTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTa                                                                                    <  1:1576542/91‑1 (MQ=255)
    gTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATGCGCTCAAATCTGAGCGTAtttt                                                                                   <  1:120019/91‑3 (MQ=38)
     tATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATc                                                                                  <  1:2714647/91‑1 (MQ=255)
      aTTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                                 <  1:2040849/91‑1 (MQ=255)
      aTTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                                 <  1:2041212/91‑1 (MQ=255)
       ttGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTAGCGAATTTTGAACGAAAAAATACGCTCAGATTTGAGCGCATTATCCa                                                                                <  1:3306739/91‑1 (MQ=255)
            aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGttt                                                                           >  1:6115/1‑91 (MQ=255)
             aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                          >  1:2507860/1‑91 (MQ=255)
                aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                       <  1:201997/91‑1 (MQ=255)
                  aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                     <  1:3676537/91‑1 (MQ=255)
                  aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                     <  1:3241460/91‑1 (MQ=255)
                  aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                     <  1:3935695/91‑1 (MQ=255)
                  aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTGTTACCTg                                                                     >  1:2076002/1‑91 (MQ=255)
                  aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTCTTACCTg                                                                     >  1:2076780/1‑91 (MQ=255)
                   aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                    >  1:2540568/1‑91 (MQ=255)
                   aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                    >  1:864920/1‑91 (MQ=255)
                     gggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTc                                                                  <  1:621645/91‑1 (MQ=255)
                     gggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTc                                                                  <  1:1102230/91‑1 (MQ=255)
                       gCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCaa                                                                <  1:3803639/91‑1 (MQ=255)
                        cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                               >  1:1846263/1‑91 (MQ=255)
                        cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                               >  1:1843798/1‑91 (MQ=255)
                        cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                               >  1:735590/1‑91 (MQ=255)
                        cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                               >  1:735523/1‑91 (MQ=255)
                          caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                             >  1:368565/1‑91 (MQ=255)
                             aaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGt                                                          >  1:1905326/1‑91 (MQ=255)
                             aaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGt                                                          <  1:2923328/91‑1 (MQ=255)
                              aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTGAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                         <  1:496907/91‑1 (MQ=255)
                               ggggCAGTATATGCAGATTATCTAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                        <  1:2434880/91‑1 (MQ=255)
                               ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                        <  1:2433944/91‑1 (MQ=255)
                               ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                        <  1:2242311/91‑1 (MQ=255)
                               ggggCAGTATAGGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                        <  1:3765351/91‑1 (MQ=255)
                               ggggCAGTATAGGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                        <  1:3764762/91‑1 (MQ=255)
                                 ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                      <  1:410825/91‑1 (MQ=255)
                                 ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                      >  1:13014/1‑91 (MQ=255)
                                   cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                    <  1:4118130/91‑1 (MQ=255)
                                   cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                    <  1:4117805/91‑1 (MQ=255)
                                   cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                    <  1:2027788/91‑1 (MQ=255)
                                    aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                                   <  1:3630626/91‑1 (MQ=255)
                                     gTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCata                                                  <  1:3616843/91‑1 (MQ=255)
                                       atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                                <  1:1725645/91‑1 (MQ=255)
                                         atGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCt                                              <  1:2244378/91‑1 (MQ=255)
                                         atGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCt                                              >  1:1572104/1‑91 (MQ=255)
                                         atGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCt                                              <  1:1311240/91‑1 (MQ=255)
                                           gCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTa                                            >  1:2779/1‑91 (MQ=255)
                                            cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCGTCATATCCTTAt                                           <  1:3190126/91‑1 (MQ=255)
                                            cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                           <  1:3688346/91‑1 (MQ=255)
                                            cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                           <  1:2434971/91‑1 (MQ=255)
                                            cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                           <  1:1609433/91‑1 (MQ=255)
                                             aaaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATc                                          <  1:2284188/91‑1 (MQ=255)
                                             aaaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATc                                          <  1:611623/91‑1 (MQ=255)
                                              aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                         >  1:3946425/1‑91 (MQ=255)
                                              aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                         <  1:74686/91‑1 (MQ=255)
                                               aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                        >  1:423419/1‑91 (MQ=255)
                                               aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                        >  1:729343/1‑91 (MQ=255)
                                                ttATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCAt                                       <  1:3843203/91‑1 (MQ=255)
                                                ttATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCAt                                       <  1:3843181/91‑1 (MQ=255)
                                                  aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCt                                     >  1:457234/1‑91 (MQ=255)
                                                   tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                    <  1:3734584/91‑1 (MQ=255)
                                                   tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                    <  1:2880664/91‑1 (MQ=255)
                                                   tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                    <  1:1010957/91‑1 (MQ=255)
                                                   tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                    <  1:1011062/91‑1 (MQ=255)
                                                    cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                                   >  1:3015335/1‑91 (MQ=255)
                                                    cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                                   >  1:1467440/1‑91 (MQ=255)
                                                    cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                                   >  1:1406838/1‑91 (MQ=255)
                                                         tttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGt                              <  1:87741/91‑1 (MQ=255)
                                                         tttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGt                              <  1:1469128/91‑1 (MQ=255)
                                                            aaaCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAAc                           >  1:4017514/1‑91 (MQ=255)
                                                              aCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                         <  1:1375699/91‑1 (MQ=255)
                                                              aCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                         <  1:2783266/91‑1 (MQ=255)
                                                               cAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                        <  1:3919242/91‑1 (MQ=255)
                                                                  aaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgaga                     >  1:1689565/1‑91 (MQ=255)
                                                                   aaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGaa                    >  1:1485478/1‑91 (MQ=255)
                                                                   aaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGaa                    >  1:1485477/1‑91 (MQ=255)
                                                                    aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                   <  1:1279998/91‑1 (MQ=255)
                                                                    aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                   <  1:2559339/91‑1 (MQ=255)
                                                                     aaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATa                  <  1:381360/91‑1 (MQ=255)
                                                                       tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa                <  1:2310180/91‑1 (MQ=255)
                                                                       tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa                <  1:4007460/91‑1 (MQ=255)
                                                                        gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt               >  1:3929530/2‑91 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              >  1:642296/1‑91 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              <  1:3490857/91‑1 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              >  1:842169/1‑91 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              <  1:3490568/91‑1 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCAGATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              <  1:3490704/91‑1 (MQ=255)
                                                                         cGCTCAGATTTGAGCGCATTATCCAAGTTCTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              >  1:205584/1‑91 (MQ=255)
                                                                         cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg              >  1:3711002/1‑91 (MQ=255)
                                                                          gCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg             >  1:2685397/1‑91 (MQ=255)
                                                                          gCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg             >  1:708688/1‑91 (MQ=255)
                                                                           cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa            >  1:1112307/1‑91 (MQ=255)
                                                                            tCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGaa           >  1:1377068/1‑91 (MQ=255)
                                                                               gATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt        <  1:2318374/91‑1 (MQ=255)
                                                                               gATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt        <  1:253312/91‑1 (MQ=255)
                                                                                 tctGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtgt      >  1:3360401/3‑91 (MQ=255)
                                                                                  ttGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgtt     >  1:178252/1‑91 (MQ=255)
                                                                                   tGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgttg    <  1:2322538/91‑1 (MQ=255)
                                                                                   tGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGtgttg    >  1:1221626/1‑91 (MQ=255)
                                                                                     aGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTa  <  1:439199/91‑1 (MQ=255)
                                                                                         |                                                                                      
TGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTTGTA  >  USA300TCH1516_ALE/895853‑896028

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