Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,780,409 A→C 5.2% Y106* (TAT→TAG arcD_2 ← Arginine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,780,4090AC 5.2% 791.8 / 28.0 229Y106* (TAT→TAGarcD_2Arginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base A (117/100);  new base C (3/9);  total (120/109)
Fisher's exact test for biased strand distribution p-value = 7.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 232 total aligned reads displayed.

AGTAACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGCATAACTA  >  USA300TCH1516_ALE/2780320‑2780499
                                                                                         |                                                                                          
agtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCAt                                                                                           <  1:3225544/91‑1 (MQ=255)
agtaACGAGGCGAAAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAATGTTACCTACTGCTGACATCAATAGTGTTGCAt                                                                                           <  2:3887337/91‑1 (MQ=255)
 gtaACGAGGCGACAATAACCCTTGGTGATGCCTTGCCTCCTGTAAAAAGCGCGAAAAAGTTACCTACTGCTCCCAGCAACAGTGTTGCATa                                                                                          >  2:2547480/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTTCTGCTGACATCAATAGTGTTGCATa                                                                                          >  2:609552/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATa                                                                                          >  1:3545059/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATa                                                                                          >  1:370361/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATa                                                                                          >  1:1669890/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATa                                                                                          >  1:581353/1‑91 (MQ=255)
 gtaACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATa                                                                                          >  1:1465489/1‑91 (MQ=255)
  taACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAg                                                                                         >  2:832701/1‑91 (MQ=255)
   aaCGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAgg                                                                                        >  1:1151921/1‑91 (MQ=255)
   aaCGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAgg                                                                                        >  1:1151395/1‑91 (MQ=255)
   aaCGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAgg                                                                                        <  1:2432652/91‑1 (MQ=255)
    aCGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGc                                                                                       <  1:1191113/91‑1 (MQ=255)
    aCGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGc                                                                                       <  1:3919330/91‑1 (MQ=255)
       aGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCCTATGCAAc                                                                                    <  1:1100254/91‑1 (MQ=255)
       aGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCaat                                                                                    <  2:2238518/91‑2 (MQ=255)
       aGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAAc                                                                                    <  1:455618/91‑1 (MQ=255)
       aGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAATGTTGCCTATGCAAc                                                                                    <  2:2385848/91‑1 (MQ=38)
         gCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACAt                                                                                  <  1:2066642/91‑1 (MQ=255)
           gACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTg                                                                                >  1:3377256/1‑91 (MQ=255)
            aCAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGc                                                                               <  1:457680/91‑1 (MQ=255)
             cAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGcc                                                                              <  1:482269/91‑1 (MQ=255)
              aaTAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCt                                                                             >  1:3073094/1‑91 (MQ=255)
              aaTAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCt                                                                             <  1:2473558/91‑1 (MQ=255)
               aTAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTa                                                                            >  1:1993686/1‑91 (MQ=255)
                   cacTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTaaaaa                                                                        <  1:1614189/91‑1 (MQ=255)
                    acTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAAc                                                                       <  1:3480025/91‑1 (MQ=255)
                    acTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAAc                                                                       >  1:1013382/1‑91 (MQ=255)
                    acTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAAc                                                                       <  2:2681567/91‑1 (MQ=255)
                    acTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAAc                                                                       >  1:1766732/1‑91 (MQ=255)
                      ttGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGc                                                                     <  1:3314197/91‑1 (MQ=255)
                       tGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCt                                                                    >  1:2577291/1‑91 (MQ=255)
                       tGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCt                                                                    <  1:2529575/91‑1 (MQ=255)
                       tGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCt                                                                    <  1:1380069/91‑1 (MQ=255)
                        ggTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTg                                                                   >  1:3740889/1‑91 (MQ=255)
                         gTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTga                                                                  >  1:941388/1‑91 (MQ=255)
                            aTGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGaa                                                               >  1:1471745/1‑91 (MQ=255)
                            aTGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGaa                                                               <  1:3198917/91‑1 (MQ=255)
                             tgtgTTGCCTCCTTTAAAAATCGGGAAAACGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAAc                                                              >  2:3448993/1‑91 (MQ=255)
                             tgtgTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCACCATTGCCTAAAAACGCTGAGAAc                                                              >  1:1306918/1‑91 (MQ=255)
                             tgtgTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAAc                                                              >  2:2261538/1‑91 (MQ=255)
                             tgtgTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAAc                                                              >  2:3791126/1‑91 (MQ=255)
                              gtgtTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAAcc                                                             <  1:254796/91‑1 (MQ=255)
                                  tGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAata                                                         >  1:263715/1‑91 (MQ=255)
                                   gCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAatag                                                        >  2:3911272/1‑90 (MQ=255)
                                    cctcctTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATc                                                       >  1:1683861/1‑91 (MQ=255)
                                     ctcctTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATcc                                                      <  1:567571/91‑1 (MQ=255)
                                      tcctTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATccc                                                     <  2:692863/91‑1 (MQ=255)
                                       cctTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATcccc                                                    <  2:184820/91‑1 (MQ=255)
                                         tttAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCaa                                                  >  1:1825357/1‑91 (MQ=255)
                                         tttAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCaa                                                  >  1:1707149/1‑91 (MQ=255)
                                          ttAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAg                                                 >  2:1550686/1‑91 (MQ=255)
                                          ttAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAg                                                 <  1:3050090/91‑1 (MQ=255)
                                             aaaaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCAc                                              <  1:753077/91‑1 (MQ=255)
                                              aaaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACt                                             >  2:2825351/1‑91 (MQ=255)
                                              aaaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACt                                             >  1:2782234/1‑91 (MQ=255)
                                              aaaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACt                                             >  2:1682647/1‑91 (MQ=255)
                                               aaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTg                                            <  2:1060050/91‑1 (MQ=255)
                                               aaTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTg                                            >  2:1727979/1‑91 (MQ=255)
                                                aTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGa                                           <  1:1803819/91‑1 (MQ=255)
                                                aTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGa                                           <  1:3414593/91‑1 (MQ=255)
                                                aTCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGa                                           <  1:1039074/91‑1 (MQ=255)
                                                 tCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGAt                                          >  1:691288/1‑91 (MQ=255)
                                                  cGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATa                                         >  1:290829/1‑91 (MQ=255)
                                                   gggAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATaa                                        <  1:2714906/91‑1 (MQ=255)
                                                       aaaaGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATcc                                    >  2:2000803/1‑91 (MQ=255)
                                                       aaaaGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATcc                                    >  1:626627/1‑91 (MQ=255)
                                                       aaaaGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATcc                                    <  2:651654/91‑1 (MQ=255)
                                                        aaaGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCt                                   <  2:885653/91‑1 (MQ=255)
                                                          aGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTAc                                 <  2:2045709/91‑1 (MQ=255)
                                                              aCCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACaaaa                             >  2:969336/1‑91 (MQ=255)
                                                               ccTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAAt                            >  2:3183246/1‑91 (MQ=255)
                                                               ccTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAAt                            >  2:2287160/1‑91 (MQ=255)
                                                                cTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAACATCCCCAAGCACTGATAAATCCCACACAATc                           >  2:39312/1‑91 (MQ=255)
                                                                 tACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCACGCACTGATAAATCCTACACAATcc                          >  2:1498165/1‑90 (MQ=255)
                                                                 tACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCg                          <  1:2931269/91‑1 (MQ=255)
                                                                  aCTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGc                         >  2:440876/1‑91 (MQ=255)
                                                                     gctgACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCaa                      <  2:2732625/91‑1 (MQ=255)
                                                                     gctgACAACAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCaa                      <  2:1251676/91‑1 (MQ=255)
                                                                      ctgACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCaaa                     >  1:2169046/1‑91 (MQ=255)
                                                                       tgACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAAt                    <  2:1436476/91‑1 (MQ=255)
                                                                       tgACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAAt                    <  2:3559732/91‑1 (MQ=255)
                                                                        gACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATc                   <  1:710178/91‑1 (MQ=255)
                                                                         aCATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATcc                  <  2:3663510/91‑1 (MQ=255)
                                                                         aCATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATcc                  >  1:160856/1‑91 (MQ=255)
                                                                           aTCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTg                <  2:335618/91‑1 (MQ=255)
                                                                           aTCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTg                >  1:3276485/1‑91 (MQ=255)
                                                                           aTCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTg                >  1:1303582/1‑91 (MQ=255)
                                                                           aTCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTg                <  1:3836072/91‑1 (MQ=255)
                                                                           aTCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTg                >  1:3751740/1‑91 (MQ=255)
                                                                             cAATAGTGTTGCATAGGCAACAATACATAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCt              <  2:420341/91‑1 (MQ=255)
                                                                               aTAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTg            >  1:142977/1‑91 (MQ=255)
                                                                               aTAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTg            >  2:2034112/1‑91 (MQ=255)
                                                                                  gtgtTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTCGAGc         >  2:546055/1‑91 (MQ=255)
                                                                                  atgtTGCCTATGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGc         >  1:3783525/2‑91 (MQ=255)
                                                                                  atgtTGCCTATGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGc         >  1:100432/2‑91 (MQ=255)
                                                                                      tGCATAGGCAACATTGCCTACAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGCATaa     >  2:369764/1‑91 (MQ=255)
                                                                                      tGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGCATaa     >  2:2366632/1‑91 (MQ=255)
                                                                                         aTAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGCATAACta  <  1:2099481/91‑1 (MQ=255)
                                                                                         |                                                                                          
AGTAACGAGGCGACAATAACACTTGGTAATGTGTTGCCTCCTTTAAAAATCGGGAAAAAGTTACCTACTGCTGACATCAATAGTGTTGCATAGGCAACATTGCCTAAAAACGCTGAGAACCAATATCCCCAAGCACTGATAAATCCTACAAAATCGCCAAATCCTGCTTGAGCATAACTA  >  USA300TCH1516_ALE/2780320‑2780499

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