Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 848,638 C→T intergenic (+643/‑237) trxB → / → USA300HOU_RS04145 Thioredoxin reductase/Nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE848,6380CT99.5% 1380.1 / ‑2.6 369intergenic (+643/‑237)trxB/USA300HOU_RS04145Thioredoxin reductase/Nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (179/187);  total (179/190)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 371 total aligned reads displayed.

AATTTCTTTTCGAAATTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCCCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACATT  >  USA300TCH1516_ALE/848548‑848724
                                                                                          |                                                                                      
aaTTTCTTTTCGAAATTCTCTGTGTTGGGGCGTACACCCCAACTCGCATTACCTTTTGATTATCTTTACGAAACTCTTTGTGTTGGGGccc                                                                                        >  2:10383/1‑91 (MQ=17)
   ttCTTTTCGAAATTCTCTGTGTTGGGGCCCACACCTCAACTCGCATTGCCTGTATAGTTTCTTTTCGTACTTCCCTGTCATGGGGCctctt                                                                                     >  2:4499267/1‑90 (MQ=17)
    tCTTTTCGAAATTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGa                                                                                    >  1:1870534/1‑91 (MQ=255)
         tCGAAATTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTaga                                                                               <  1:5522850/91‑1 (MQ=255)
            aaaTTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGtt                                                                            >  1:6887264/1‑91 (MQ=255)
             aaTTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGCATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGttt                                                                           >  2:6778347/1‑90 (MQ=255)
              aTTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGa                                                                          <  1:2814612/91‑1 (MQ=255)
              aTTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGa                                                                          <  2:2674637/91‑1 (MQ=255)
                 ctctGTGTTGGGGCCCACACCCCACCTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaa                                                                       <  1:3848925/91‑1 (MQ=255)
                 ctctGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaa                                                                       <  2:3844101/91‑1 (MQ=255)
                 ctctGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaa                                                                       <  2:1866433/91‑1 (MQ=255)
                  tctGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaaa                                                                      >  2:5498190/1‑91 (MQ=255)
                  tctGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaaa                                                                      <  1:3141858/91‑1 (MQ=255)
                   ctGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGaaaaaa                                                                     >  1:2070317/1‑91 (MQ=255)
                    tgtgTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAg                                                                    >  2:2278058/1‑91 (MQ=255)
                     gtgtTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGc                                                                   >  1:1408113/1‑91 (MQ=255)
                     gtgtTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGc                                                                   <  1:6345876/91‑1 (MQ=255)
                      tgtTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCt                                                                  >  2:413154/1‑91 (MQ=255)
                         tGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCttgt                                                               <  2:6027200/91‑1 (MQ=255)
                          ggggCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCttgtt                                                              >  1:4469219/1‑91 (MQ=255)
                          ggggCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCttgtt                                                              >  2:1699047/1‑91 (MQ=255)
                            ggCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGc                                                            >  1:1030054/1‑91 (MQ=255)
                            ggCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTATCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGc                                                            <  1:5813996/91‑1 (MQ=255)
                             gCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCa                                                           >  1:2538581/1‑91 (MQ=255)
                             gCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCa                                                           <  2:5375141/91‑1 (MQ=255)
                              cccACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCaa                                                          <  1:6611223/91‑1 (MQ=255)
                              cccACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCaa                                                          >  1:5924837/1‑91 (MQ=255)
                               ccccacCCCACCTCGCTTTGCCTGTAGATTTTCTTTTCGAATTTCTCTGGGTTGGGCCCTCTGCCTAGAGTTAAAAAAAGCTTGTTGCAAg                                                         <  2:5064357/88‑1 (MQ=17)
                               ccACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAg                                                         <  1:5270663/91‑1 (MQ=255)
                               ccACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAg                                                         >  1:1376310/1‑91 (MQ=255)
                               ccACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAg                                                         <  1:3071446/91‑1 (MQ=255)
                                cacaCCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCCCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgc                                                        <  1:1645541/91‑1 (MQ=38)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       >  2:5207683/1‑91 (MQ=37)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       >  2:4291957/1‑91 (MQ=255)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       >  2:2432174/1‑91 (MQ=255)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       <  1:4235030/91‑1 (MQ=255)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       >  2:566595/1‑91 (MQ=255)
                                 acacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcg                                                       >  1:4059890/1‑91 (MQ=255)
                                  cacCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAgcgc                                                      <  1:4534845/91‑1 (MQ=255)
                                   acCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCa                                                     >  2:205785/1‑91 (MQ=255)
                                     cccAACTCGCATTTCCTGTACAATTTCTTTTCGAAATTCCCTGTGTTGGGGCCTCTGCCTAGAGTTGAAAAAAGCTTGTTGCAAGCGCAtt                                                   >  2:3713782/1‑91 (MQ=11)
                                      ccAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCAttt                                                  <  2:1486122/91‑1 (MQ=255)
                                      ccAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCAttt                                                  >  1:4974110/1‑91 (MQ=255)
                                      ccAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCAttt                                                  >  1:2396912/1‑91 (MQ=255)
                                      ccAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATGCTCTGTGTTGGGGCCTCTGACTAGAGGTGGAAAAAGCTTGTTGCAAGCGCAttt                                                  <  1:2969087/91‑1 (MQ=11)
                                       cAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGCCTAGAGTTGAAAAAAGCTTGTTGCAAGCGCAtttt                                                 >  2:5833199/1‑91 (MQ=255)
                                         aCTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTttca                                               <  1:4575396/91‑1 (MQ=255)
                                          cTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTttcat                                              <  1:5991973/91‑1 (MQ=255)
                                            cGCATTGTCTGTAGAATTTCTCTTCGAAATTCTCTGTGTTGGGGCCTCTTACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTttcattc                                            <  2:467225/91‑1 (MQ=14)
                                              cATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAg                                          >  2:2096730/1‑91 (MQ=255)
                                              cATGGCCTGTAGAATTTCGTTTCGAAATTCTCTGTGTTGGGGCCTCTGAATAAAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAg                                          <  2:735109/91‑1 (MQ=14)
                                               aTTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGt                                         >  2:4015648/1‑91 (MQ=255)
                                               aTTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGt                                         <  1:1725594/91‑1 (MQ=255)
                                               aTTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGt                                         <  2:1928565/91‑1 (MQ=255)
                                                ttGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTc                                        <  2:4741543/91‑1 (MQ=255)
                                                ttGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTc                                        <  2:1096722/91‑1 (MQ=255)
                                                ttGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTc                                        >  1:3598230/1‑91 (MQ=255)
                                                  gCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCaa                                      <  1:991839/91‑1 (MQ=255)
                                                   ccTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAac                                     <  2:4881667/91‑1 (MQ=255)
                                                    cTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAact                                    >  1:772554/1‑91 (MQ=255)
                                                      gTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAactac                                  >  2:2577581/1‑91 (MQ=255)
                                                      gTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAactac                                  >  2:3391020/1‑91 (MQ=255)
                                                       tAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAactact                                 <  2:3502611/91‑1 (MQ=255)
                                                          aaTTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGc                              <  2:5263350/91‑1 (MQ=255)
                                                           aTTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCa                             <  2:6421737/91‑1 (MQ=255)
                                                           aTTTCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCa                             <  2:3212343/91‑1 (MQ=255)
                                                             ttCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaat                           <  2:6953578/91‑1 (MQ=255)
                                                             ttCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaat                           >  2:3087217/1‑91 (MQ=255)
                                                              tCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGTTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaata                          >  2:2824876/1‑91 (MQ=255)
                                                              tCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaata                          <  1:1073683/91‑1 (MQ=255)
                                                              tCTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaata                          >  2:2758736/1‑91 (MQ=255)
                                                               cTTTTCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaatat                         <  1:5610538/91‑1 (MQ=255)
                                                                ttttCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaatata                        >  1:122809/1‑91 (MQ=255)
                                                                ttttCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaatata                        >  2:662957/1‑91 (MQ=255)
                                                                   tCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaatataata                     <  2:955733/91‑1 (MQ=255)
                                                                   tCGAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCaatataata                     >  2:1789685/1‑91 (MQ=255)
                                                                     gAAATTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAAtatt                   >  1:2345689/1‑91 (MQ=255)
                                                                      aaaTTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAAtatta                  <  1:4084313/91‑1 (MQ=255)
                                                                       aaTTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAAtattat                 <  2:5927136/91‑1 (MQ=255)
                                                                       aaTTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAAtattat                 <  1:1554778/91‑1 (MQ=255)
                                                                        aTTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATtata                <  2:248521/91‑1 (MQ=255)
                                                                        aTTCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATtata                >  1:675281/1‑91 (MQ=255)
                                                                         ttCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAg               <  1:6265505/91‑1 (MQ=255)
                                                                         ttCTCTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAg               >  2:6348560/1‑91 (MQ=255)
                                                                          tctcTGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGa              >  2:6183448/1‑91 (MQ=255)
                                                                            tctGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGAcc            <  1:5053434/91‑1 (MQ=255)
                                                                             ctGTGTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGAccc           >  2:1274744/1‑91 (MQ=255)
                                                                              tgtgTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCt          >  1:4729814/1‑91 (MQ=255)
                                                                              tgtgTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCt          >  2:3800566/1‑91 (MQ=255)
                                                                              tgtgTTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCt          >  1:3371851/1‑91 (MQ=255)
                                                                               gtgtTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTa         <  1:741537/91‑1 (MQ=255)
                                                                                tgtTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAg        <  2:2875921/91‑1 (MQ=255)
                                                                                 gtTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAgg       <  1:2504231/91‑1 (MQ=255)
                                                                                 gtTGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAgg       <  1:3511859/91‑1 (MQ=255)
                                                                                  ttGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGa      <  2:5860654/91‑1 (MQ=255)
                                                                                  ttGGGGCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGa      >  1:226319/1‑91 (MQ=255)
                                                                                    ggggCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACa    >  1:2253534/1‑91 (MQ=255)
                                                                                    ggggCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACa    <  1:4314490/91‑1 (MQ=255)
                                                                                    ggggCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACa    <  1:4394620/91‑1 (MQ=255)
                                                                                      ggCCTCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACAtt  <  1:958437/91‑1 (MQ=255)
                                                                                          |                                                                                      
AATTTCTTTTCGAAATTCTCTGTGTTGGGGCCCACACCCCAACTCGCATTGCCTGTAGAATTTCTTTTCGAAATTCTCTGTGTTGGGGCCCCTGACTAGAGTTGAAAAAAGCTTGTTGCAAGCGCATTTTCATTCAGTCAACTACTAGCAATATAATATTATAGACCCTAGGACATT  >  USA300TCH1516_ALE/848548‑848724

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: