Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,158,427 G→A 9.6% intergenic (+11/‑313) 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,4270GA 9.6% 708.0 / 45.2 236intergenic (+11/‑313)uvrC/sdhCUvrABC system protein C/Succinate dehydrogenase cytochrome b558 subunit
Reads supporting (aligned to +/- strand):  ref base G (107/106);  new base A (4/19);  total (111/125)
Fisher's exact test for biased strand distribution p-value = 3.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 237 total aligned reads displayed.

AAGCTACACTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTT  >  USA300TCH1516_ALE/1158337‑1158516
                                                                                          |                                                                                         
aagCTACACTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCg                                                                                           <  2:125718/91‑1 (MQ=255)
  gCTACACTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTg                                                                                         <  1:3350640/91‑1 (MQ=255)
  gCTACACTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTg                                                                                         <  1:444526/91‑1 (MQ=255)
    tACACTTGAAGATTTTAAAAATATAGGTATTCCTGACAACGGCGCACAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTgag                                                                                       >  2:169364/1‑91 (MQ=255)
      cacTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGaa                                                                                     >  2:145750/1‑91 (MQ=255)
      cacTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGaa                                                                                     >  1:3647865/1‑91 (MQ=255)
       acTTGAAGATTTTAAAAATATAGGTATTCCTGAATACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAAt                                                                                    <  2:405693/91‑1 (MQ=255)
       acTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAAt                                                                                    <  2:1139017/91‑1 (MQ=255)
       acTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAAt                                                                                    <  1:840367/91‑1 (MQ=255)
       acTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTCCATAAATAAAATACTGGGCGCGAGaaa                                                                                    >  2:39006/1‑90 (MQ=255)
          tGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTGGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATGAt                                                                                 <  2:3112300/91‑1 (MQ=255)
          tGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAAtgat                                                                                 <  1:1984497/91‑4 (MQ=255)
          tGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAAtgat                                                                                 <  1:80631/91‑4 (MQ=255)
           gAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGAtt                                                                                <  2:647639/91‑1 (MQ=255)
           gAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGAtt                                                                                <  1:557474/91‑1 (MQ=255)
            aaGATTTTAAAAATATAGGTATTCCTGAAAACTTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGTGCATGAGAATGATtc                                                                               <  1:3368673/91‑1 (MQ=255)
            aaGAGTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATGATtc                                                                               <  1:3287567/91‑1 (MQ=255)
             aGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATtct                                                                              <  1:2416794/91‑1 (MQ=255)
               aTTTTAAAAATATAGGTATTCCTGACACCGTCGCAAAGAACCCACATGAACAATTGCATCAATAAATTACTGGGCGTGAGActcctcctcc                                                                            >  2:3743756/1‑90 (MQ=255)
                ttttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCAt                                                                           >  1:2836598/1‑91 (MQ=255)
                  ttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGc                                                                         >  1:3770176/1‑91 (MQ=255)
                  ttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGc                                                                         <  1:1946672/91‑1 (MQ=255)
                   tacaaaTATAGGTACGCCTGAAAACGTCGCAAAGAGCCTACATGAACAATTGCATAAATAAATTACTGTGCGTGAGAATCATTCTCAGGcc                                                                        <  2:3340411/88‑1 (MQ=255)
                   tAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATCCTCATGcc                                                                        >  2:3653999/1‑91 (MQ=255)
                    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                       <  2:2202592/91‑1 (MQ=255)
                    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                       >  1:1744087/1‑91 (MQ=255)
                    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                       >  1:1204385/1‑91 (MQ=255)
                     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                      >  1:3831665/1‑91 (MQ=255)
                     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                      <  2:3189584/91‑1 (MQ=255)
                      aaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATGGCATAAATAAATTACTGGGCATGAGAATGATTCTCACGcccag                                                                     <  2:3054925/91‑3 (MQ=38)
                       aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaa                                                                    >  1:1383629/1‑91 (MQ=255)
                        atatAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaa                                                                   >  2:1623280/1‑91 (MQ=255)
                         tataGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaaa                                                                  <  2:2589933/91‑1 (MQ=255)
                             ggTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAatata                                                              <  1:1133532/91‑1 (MQ=255)
                               tATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATg                                                            <  1:2913434/91‑1 (MQ=255)
                                aTTCCTGAAAACGTCGCAGAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                           <  1:860066/91‑1 (MQ=255)
                                aTTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                           >  2:362406/1‑91 (MQ=255)
                                aTTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                           <  2:2771631/91‑1 (MQ=255)
                                 ttCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                          <  2:1939022/91‑1 (MQ=255)
                                 ttCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                          <  2:3148357/91‑1 (MQ=255)
                                   ccTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaat                                                        >  2:311453/1‑91 (MQ=255)
                                    cTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaata                                                       >  1:3905468/1‑91 (MQ=255)
                                     tGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataa                                                      >  1:3636282/1‑91 (MQ=255)
                                      gAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataat                                                     <  2:359517/91‑1 (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)
                                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCTAAAATATATGAAATAATGTTACAat                                             >  1:3691445/1‑91 (MQ=255)
                                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                             >  1:3109301/1‑91 (MQ=255)
                                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                             >  2:3713536/1‑91 (MQ=255)
                                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                             >  2:860155/1‑91 (MQ=255)
                                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                             >  2:2557016/1‑91 (MQ=255)
                                               gCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAata                                            >  1:2930636/1‑91 (MQ=255)
                                                cAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAatat                                           <  2:459386/91‑1 (MQ=255)
                                                  aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                         <  2:2594869/91‑1 (MQ=255)
                                                   aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                        >  2:1665923/1‑91 (MQ=255)
                                                   aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                        <  2:2843159/91‑1 (MQ=255)
                                                     aaCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTa                                      <  1:297861/91‑1 (MQ=255)
                                                      aCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTaa                                     >  2:3370459/1‑91 (MQ=255)
                                                      aCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTaa                                     >  1:39234/1‑91 (MQ=255)
                                                         tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                  >  2:2392225/1‑91 (MQ=255)
                                                         tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                  >  1:2617838/1‑91 (MQ=255)
                                                         tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                  >  1:32290/1‑91 (MQ=255)
                                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                 <  1:2686765/91‑1 (MQ=255)
                                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                 >  2:3267486/1‑91 (MQ=255)
                                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                 <  2:1022595/91‑1 (MQ=255)
                                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                               <  1:1003962/91‑1 (MQ=255)
                                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                               <  2:2784554/91‑1 (MQ=255)
                                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                               <  2:345523/91‑1 (MQ=255)
                                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                               <  1:3259560/91‑1 (MQ=255)
                                                             tGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGtt                              >  2:1001970/1‑91 (MQ=255)
                                                              gAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAGGTTa                             <  1:2223694/91‑1 (MQ=255)
                                                               aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                            >  1:1614028/1‑91 (MQ=255)
                                                                 cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGtttttt                          >  2:1301982/1‑91 (MQ=255)
                                                                 cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttt                          >  2:2925515/1‑91 (MQ=255)
                                                                  aaTTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttt                         >  2:675655/1‑91 (MQ=255)
                                                                  aaTTGCATAAATAAATTACTGGGCGTGAGAATCATGCTCATGCCCCAAAATATATGAAGGAATGTTACAATATGGTTAACAAAGTTAtttt                         <  2:2939855/91‑1 (MQ=255)
                                                                    ttGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttttt                       >  2:585287/1‑91 (MQ=255)
                                                                    ttGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttttt                       >  1:628834/1‑91 (MQ=255)
                                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                    >  2:2386852/1‑91 (MQ=255)
                                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                    >  2:1221051/1‑91 (MQ=255)
                                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGGTAACAAAGTTATTTTTTTGa                    >  1:2463053/1‑91 (MQ=255)
                                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                   >  1:1764622/1‑91 (MQ=255)
                                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                   >  2:1475604/1‑91 (MQ=255)
                                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                   <  1:973220/91‑1 (MQ=255)
                                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTATGaa                   <  1:3122323/91‑1 (MQ=255)
                                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg           >  2:961563/1‑91 (MQ=255)
                                                                                 ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg          <  2:1184782/91‑1 (MQ=255)
                                                                                 ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg          >  1:3601857/1‑91 (MQ=255)
                                                                                 ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg          <  1:1488600/91‑1 (MQ=255)
                                                                                  tACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGa         <  2:1317360/91‑1 (MQ=255)
                                                                                   aCTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa        <  2:2112858/91‑1 (MQ=255)
                                                                                    cTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAg       >  2:1874185/1‑91 (MQ=255)
                                                                                     tGGGCGTGAGAATCATTCTCATGTCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc      >  2:1117260/1‑91 (MQ=255)
                                                                                     ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc      >  2:1750655/2‑91 (MQ=38)
                                                                                      gggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCg     <  1:746816/91‑1 (MQ=255)
                                                                                      gggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCg     <  1:1431779/91‑1 (MQ=255)
                                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt    <  1:2407120/91‑1 (MQ=255)
                                                                                          |                                                                                         
AAGCTACACTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTT  >  USA300TCH1516_ALE/1158337‑1158516

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