Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,779,810 A→T 6.2% intergenic (‑116/+26) USA300TCH1516_01664 ← / ← rplT hypothetical protein/50S ribosomal protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,779,8100AT 6.2% 685.9 / 24.5 209intergenic (‑116/+26)USA300TCH1516_01664/rplThypothetical protein/50S ribosomal protein L20
Reads supporting (aligned to +/- strand):  ref base A (94/101);  new base T (3/10);  total (97/112)
Fisher's exact test for biased strand distribution p-value = 9.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 226 total aligned reads displayed.

TTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCG  >  USA300TCH1516_ALE/1779720‑1779899
                                                                                          |                                                                                         
ttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACaa                                                                                           <  2:1969058/91‑1 (MQ=255)
 tgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAAt                                                                                          >  1:1143575/1‑91 (MQ=255)
    tctcTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                       <  1:497549/91‑1 (MQ=255)
    tctcTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                       <  2:2440370/91‑1 (MQ=255)
     ttctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                      <  2:2751117/90‑1 (MQ=255)
     ctctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                      <  2:201088/91‑1 (MQ=255)
      tctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAt                                                                                     <  2:2856354/91‑1 (MQ=255)
           aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTTTTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                                <  1:3279093/91‑1 (MQ=38)
           aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                                <  1:771073/91‑1 (MQ=38)
            aTGCATTTTAGCAAATATGTCATATTTTCGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCttt                                                                               >  2:397853/1‑91 (MQ=255)
             tGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttt                                                                              >  2:358556/1‑91 (MQ=255)
               cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtttttt                                                                            <  2:2356791/91‑1 (MQ=38)
               cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                            <  1:948581/91‑1 (MQ=255)
               cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                            >  2:3392836/1‑91 (MQ=255)
                aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCTTTTTta                                                                           <  1:821915/91‑1 (MQ=37)
                aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCTTTTTta                                                                           <  2:2937330/91‑1 (MQ=37)
                 ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                          <  1:3433377/91‑1 (MQ=255)
                 ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                          <  1:2365936/91‑1 (MQ=255)
                   ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                        >  2:1205011/1‑91 (MQ=255)
                   ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGATATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                        <  2:3689153/91‑1 (MQ=255)
                   ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTCGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                        >  2:2948738/1‑91 (MQ=255)
                    tAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAg                                                                       >  2:1540237/1‑91 (MQ=255)
                      gCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtg                                                                     >  1:2912436/1‑91 (MQ=255)
                      gCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtg                                                                     <  2:2395073/91‑1 (MQ=255)
                       cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                    >  2:2220826/1‑91 (MQ=255)
                         aagatGTCGTATGTTAGCGATTAAAAGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtta                                                                  <  2:1599967/88‑1 (MQ=255)
                         aaTATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtta                                                                  <  2:2264583/91‑1 (MQ=255)
                          atatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttat                                                                 <  2:1186505/91‑1 (MQ=255)
                          atatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTATTTATAAGTGttat                                                                 >  2:534857/1‑91 (MQ=255)
                            atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                               <  2:1971166/91‑1 (MQ=255)
                            atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                               <  1:2873960/91‑1 (MQ=255)
                            atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                               >  2:1087291/1‑91 (MQ=255)
                              gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATCAGTGTTattatt                                                             >  2:1670468/1‑91 (MQ=255)
                              gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTattatt                                                             <  2:2086426/91‑1 (MQ=255)
                              gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTattatt                                                             <  1:2808551/91‑1 (MQ=255)
                                cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTAtttt                                                           >  2:3315582/1‑91 (MQ=255)
                                cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTAtttt                                                           >  2:781747/1‑91 (MQ=255)
                                 atatTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTa                                                          >  1:1001267/1‑91 (MQ=255)
                                 atatTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTa                                                          <  1:2565935/91‑1 (MQ=255)
                                    ttttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAg                                                       <  2:3014212/91‑1 (MQ=255)
                                    ttttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAg                                                       <  1:213977/91‑1 (MQ=255)
                                      ttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCa                                                     <  2:2673098/91‑1 (MQ=255)
                                        aGCGATTAAATGTGGCATTTTGCTATCTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                   >  2:2205848/1‑91 (MQ=255)
                                           gATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATcttt                                                <  1:411672/91‑1 (MQ=255)
                                             ttAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttag                                              <  1:1801710/91‑1 (MQ=255)
                                                aaTGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctt                                           >  1:3555309/1‑91 (MQ=255)
                                                   gtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctttag                                        >  2:3052447/1‑91 (MQ=255)
                                                   atgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTGAAAGCATctttagctttag                                        <  2:3452471/90‑1 (MQ=255)
                                                     ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                      <  1:671180/91‑1 (MQ=255)
                                                     ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                      <  2:3319670/91‑1 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                     <  2:886658/91‑1 (MQ=255)
                                                      gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                     <  1:1175461/91‑1 (MQ=255)
                                                       cATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                    <  1:2931725/91‑1 (MQ=255)
                                                         ttttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTa                                  >  1:975119/1‑91 (MQ=255)
                                                          tttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTaa                                 <  1:1515457/91‑1 (MQ=255)
                                                           ttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAAt                                >  1:1764414/1‑91 (MQ=255)
                                                            tGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGAGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAGTAAtt                               >  1:3615144/1‑91 (MQ=255)
                                                             gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                              >  1:779968/1‑91 (MQ=255)
                                                                 tttCAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTGACTAATTGAGCa                          <  1:3000106/91‑1 (MQ=255)
                                                                 tttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCa                          >  1:1740022/1‑91 (MQ=255)
                                                                   taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                        <  2:2056034/91‑1 (MQ=255)
                                                                    aaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAAt                       <  2:130361/91‑1 (MQ=255)
                                                                    aaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAAt                       >  1:574668/1‑91 (MQ=255)
                                                                     aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATACGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                      >  2:3194470/1‑91 (MQ=255)
                                                                     aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAATCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                      >  2:1368718/1‑91 (MQ=255)
                                                                      aaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGc                     <  2:1421196/91‑1 (MQ=255)
                                                                       aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTCTTTTACAGCACCTTTAGCCTTAGTTACTCATTGAGCAAATCCt                    >  2:3460400/1‑91 (MQ=255)
                                                                       aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCt                    <  1:1723318/91‑1 (MQ=255)
                                                                         aataaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCttt                  >  2:1444034/4‑91 (MQ=255)
                                                                          taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt                 >  1:1361026/1‑91 (MQ=255)
                                                                          taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt                 >  1:3684082/1‑91 (MQ=255)
                                                                          taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt                 <  2:1205471/91‑1 (MQ=255)
                                                                            aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc               >  1:3745018/1‑91 (MQ=255)
                                                                             aaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg              >  1:2609463/1‑91 (MQ=255)
                                                                              aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt             >  2:3735376/1‑91 (MQ=255)
                                                                              aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt             >  2:1404288/1‑91 (MQ=255)
                                                                              aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt             <  1:707207/91‑1 (MQ=255)
                                                                               aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc            <  2:831790/91‑1 (MQ=255)
                                                                                agagATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa           >  2:2379464/1‑91 (MQ=38)
                                                                                agagATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa           >  1:2513862/1‑91 (MQ=255)
                                                                                  agaTGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGa         <  2:1726765/91‑1 (MQ=255)
                                                                                    tttaCAGTGCCTTGTCTTTTTTATATGGGTTATTATTTTAAAGCTTCTTTAGCTTTAGTTATTAATTGAGCAAATGCTTTTTCGTCAGaaa       <  2:2693800/88‑1 (MQ=255)
                                                                                     tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGaatt      >  2:3574175/1‑89 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt     <  2:628404/91‑1 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt     <  2:1143575/91‑1 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt     <  2:3002133/91‑1 (MQ=255)
                                                                                      gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt     <  2:3751239/91‑1 (MQ=255)
                                                                                       aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcg    <  1:1633669/91‑3 (MQ=255)
                                                                                       aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcg    <  1:1974938/91‑3 (MQ=255)
                                                                                       aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg    >  1:2589278/1‑91 (MQ=255)
                                                                                       aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg    <  1:1350076/91‑1 (MQ=255)
                                                                                        cAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgc   <  2:172079/91‑4 (MQ=255)
                                                                                        cAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGc   >  2:1894183/1‑91 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgca  <  1:2407226/91‑5 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcgca  <  2:706927/91‑5 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCg  >  2:929595/1‑91 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCg  <  2:448056/91‑1 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCg  <  2:1321032/91‑1 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCg  <  1:2948738/91‑1 (MQ=255)
                                                                                         aaTGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATCGca  <  1:1115605/91‑2 (MQ=255)
                                                                                          |                                                                                         
TTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTGCG  >  USA300TCH1516_ALE/1779720‑1779899

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: