Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,158,443 T→C 7.7% intergenic (+27/‑297) uvrC → / → sdhC UvrABC system protein C/Succinate dehydrogenase cytochrome b558 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,158,4430TC 7.7% 760.3 / 33.9 253intergenic (+27/‑297)uvrC/sdhCUvrABC system protein C/Succinate dehydrogenase cytochrome b558 subunit
Reads supporting (aligned to +/- strand):  ref base T (121/111);  new base C (18/2);  total (139/114)
Fisher's exact test for biased strand distribution p-value = 8.05e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 255 total aligned reads displayed.

TTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGTT  >  USA300TCH1516_ALE/1158353‑1158532
                                                                                          |                                                                                         
   tAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcc                                                                                        <  2:1188191/91‑1 (MQ=255)
   tAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATCCTCATGcc                                                                                        >  2:3653999/1‑91 (MQ=255)
    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                       <  2:2202592/91‑1 (MQ=255)
    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                       >  1:1204385/1‑91 (MQ=255)
    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                       >  1:1483368/1‑91 (MQ=255)
    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                       >  1:1744087/1‑91 (MQ=255)
     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                                      <  1:993135/91‑1 (MQ=255)
     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                                      <  2:1500696/91‑1 (MQ=255)
     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                                      >  1:3831666/1‑91 (MQ=255)
       aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaa                                                                                    >  1:1383629/1‑91 (MQ=255)
        atatAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaa                                                                                   >  2:3917654/1‑91 (MQ=255)
         tataGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaaa                                                                                  <  2:3114254/91‑1 (MQ=255)
         tataGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaaa                                                                                  <  2:2589933/91‑1 (MQ=255)
         tataGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaaa                                                                                  <  2:3731249/91‑1 (MQ=255)
               tATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATg                                                                            >  1:3661977/1‑91 (MQ=255)
                aTTCCTGAAAACGTCGCAGAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                                           <  1:859958/91‑1 (MQ=255)
                aTTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                                           >  2:362406/1‑91 (MQ=255)
                 ttCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                                          <  2:3148357/91‑1 (MQ=255)
                     tGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataa                                                                      >  1:2915120/1‑91 (MQ=255)
                       aaaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATCACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATg                                                                    >  2:3052634/1‑91 (MQ=255)
                        aaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGt                                                                   <  1:2821736/91‑1 (MQ=255)
                         aaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGtt                                                                  <  2:2026614/91‑1 (MQ=255)
                         aaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGtt                                                                  <  1:3232037/91‑1 (MQ=255)
                          aCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTa                                                                 >  2:3011310/1‑91 (MQ=255)
                           cGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTAc                                                                >  1:2325989/1‑91 (MQ=255)
                             tCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACaa                                                              >  1:2935832/1‑91 (MQ=255)
                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCTAAAATATATGAAATAATGTTACAat                                                             >  1:3691445/1‑91 (MQ=255)
                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                             >  2:860155/1‑91 (MQ=255)
                               gCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAata                                                            >  2:3232116/1‑91 (MQ=255)
                                 aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                          <  2:1442583/91‑1 (MQ=255)
                                 aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                          <  2:3536040/91‑1 (MQ=255)
                                  aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                                         <  2:2594869/91‑1 (MQ=255)
                                  aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                                         >  2:1243659/1‑91 (MQ=255)
                                   aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                                        >  2:1665923/1‑91 (MQ=255)
                                   aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                                        <  2:2843159/91‑1 (MQ=255)
                                   aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                                        <  1:532723/91‑1 (MQ=255)
                                    gAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCCTGCCCCCAAATATATGAAATAATGTTACAACATGGtt                                                       >  1:694119/1‑91 (MQ=255)
                                      aCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTaa                                                     >  2:3370459/1‑91 (MQ=255)
                                       ccTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAAc                                                    >  2:503545/1‑91 (MQ=255)
                                         tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                                  >  2:2392225/1‑91 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                               <  2:1094710/91‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                               <  1:1003962/91‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                               <  2:2201040/91‑1 (MQ=255)
                                              gAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTa                                             >  1:2647212/1‑91 (MQ=255)
                                               aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                            >  1:2176290/1‑91 (MQ=255)
                                               aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                            >  1:1614004/1‑91 (MQ=255)
                                                 cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttt                                          >  2:2925515/1‑91 (MQ=255)
                                                  aaTTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttt                                         >  2:675655/1‑91 (MQ=255)
                                                   aTTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttt                                        >  1:913306/1‑91 (MQ=255)
                                                    ttGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttttt                                       >  2:585287/1‑91 (MQ=255)
                                                     tGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttttt                                      >  2:306594/1‑91 (MQ=255)
                                                     tGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAACATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttttt                                      >  1:3944178/1‑91 (MQ=255)
                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                                    >  2:2386852/1‑91 (MQ=255)
                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                                   >  1:1764622/1‑91 (MQ=255)
                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                                   >  1:1673665/1‑91 (MQ=255)
                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTATGaa                                   <  1:3122323/91‑1 (MQ=255)
                                                          aaataaatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATg                                 <  2:1122000/91‑1 (MQ=255)
                                                           aataaatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGt                                >  1:2859773/1‑91 (MQ=255)
                                                             taaatTACTGAGCGTGTGATTCATTCTCATGATCCGTATTATATGAAATAAGGTTACGATATGGTTAACAAAGTTATTTTTGGGAATGTAc                              <  2:454356/91‑1 (MQ=255)
                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                           >  1:2429481/1‑91 (MQ=255)
                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                           >  1:2538758/1‑91 (MQ=255)
                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                           >  1:3229858/1‑91 (MQ=255)
                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                           >  2:961563/1‑91 (MQ=255)
                                                                 ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg                          >  1:3602192/1‑91 (MQ=255)
                                                                  tACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGa                         <  2:1317360/91‑1 (MQ=255)
                                                                   ttggggCATGAGAATGATTCTCACGCCCCAATATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTCGAATGTACGTAGaa                        >  1:3564139/4‑91 (MQ=25)
                                                                   ttggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa                        >  1:3286116/4‑91 (MQ=38)
                                                                    cTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAg                       >  2:1874185/1‑91 (MQ=255)
                                                                     ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc                      >  2:1750655/2‑91 (MQ=38)
                                                                     ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc                      >  1:765977/2‑91 (MQ=38)
                                                                      gggCGTGAGAATCATTCTCATGCCCCAAAAGATATGAAAGAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCg                     <  1:3914391/91‑1 (MQ=255)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                    <  2:1011617/91‑1 (MQ=255)
                                                                            gagaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCt               <  1:2766837/91‑1 (MQ=255)
                                                                            gagaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCt               <  2:3335692/91‑1 (MQ=255)
                                                                             agaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTa              <  1:3057013/91‑1 (MQ=255)
                                                                             agaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGCTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTa              >  2:3496345/1‑91 (MQ=255)
                                                                             agaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGCTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTa              >  2:2788481/1‑91 (MQ=255)
                                                                               aatgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa            >  2:1692229/5‑91 (MQ=255)
                                                                               aatgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa            >  2:3297882/5‑91 (MQ=255)
                                                                               aatgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa            >  1:826082/5‑91 (MQ=255)
                                                                               aaTCATTCTCATGCCACAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa            <  2:3855073/91‑1 (MQ=255)
                                                                                atgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGCATGTACGTGGACGCGTTTTCTACAt           >  2:3408749/4‑91 (MQ=255)
                                                                                atgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAAt           >  1:1908583/4‑91 (MQ=255)
                                                                                aTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAAt           >  1:1101243/1‑91 (MQ=255)
                                                                                aTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAAt           >  1:3354691/1‑91 (MQ=255)
                                                                                aTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAAt           >  1:1442460/1‑91 (MQ=255)
                                                                                 tCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATg          >  1:593146/1‑91 (MQ=255)
                                                                                  cATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATga         >  2:1966246/1‑91 (MQ=255)
                                                                                   aTTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgag        <  2:1812105/91‑1 (MQ=255)
                                                                                   aTTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgag        >  2:980404/1‑91 (MQ=255)
                                                                                    ttCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgaga       >  2:2807999/1‑91 (MQ=255)
                                                                                    ttCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgaga       >  1:3831088/1‑91 (MQ=255)
                                                                                     tctcATGCCACAAAATATGTGAAGTAATGTTACAAAATGGTTAACAAAGTTAGTTTTTTGGATGTACGAGTCAGCATTTTCTAAATTAGaa      <  2:293861/91‑1 (MQ=255)
                                                                                      ctcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGaaa     <  2:3264277/91‑1 (MQ=255)
                                                                                       tcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAg    <  2:2734655/91‑1 (MQ=255)
                                                                                       tcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAg    <  2:148444/91‑1 (MQ=255)
                                                                                       tcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAg    <  1:3396554/91‑1 (MQ=255)
                                                                                       tcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAg    >  1:2061023/1‑91 (MQ=255)
                                                                                        cATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGt   >  2:276186/1‑91 (MQ=255)
                                                                                         aTGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGtt  <  1:981985/91‑1 (MQ=255)
                                                                                          |                                                                                         
TTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGTT  >  USA300TCH1516_ALE/1158353‑1158532

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: