Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,158,443 T→C 5.2% 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 5.2% 647.7 / 16.2 203intergenic (+27/‑297)uvrC/sdhCUvrABC system protein C/Succinate dehydrogenase cytochrome b558 subunit
Reads supporting (aligned to +/- strand):  ref base T (97/95);  new base C (8/3);  total (105/98)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 207 total aligned reads displayed.

TTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGTTT  >  USA300TCH1516_ALE/1158353‑1158533
                                                                                          |                                                                                          
ttttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCAt                                                                                            >  1:1359320/1‑91 (MQ=255)
 tttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATg                                                                                           <  1:1095162/91‑1 (MQ=255)
 tttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATg                                                                                           <  1:530318/91‑1 (MQ=255)
   tAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcc                                                                                         <  2:1622447/91‑1 (MQ=255)
   tAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcc                                                                                         <  1:965853/91‑1 (MQ=255)
    aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                        <  1:1536080/91‑1 (MQ=255)
     aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                                       <  1:3149248/91‑1 (MQ=255)
      aaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCa                                                                                      <  1:1315970/91‑1 (MQ=255)
       aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaa                                                                                     <  1:2862449/91‑1 (MQ=255)
       aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGATTCTCACGcccagt                                                                                     <  1:3010929/91‑4 (MQ=38)
        atatAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGATTCTCACGcccagta                                                                                    <  2:1030040/91‑5 (MQ=37)
           taGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAata                                                                                 >  2:651754/1‑91 (MQ=255)
             ggTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAatata                                                                               >  2:25749/1‑91 (MQ=255)
                 ttCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                                           <  2:3256823/91‑1 (MQ=255)
                 ttCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                                           <  1:2152889/91‑1 (MQ=255)
                   ccTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaat                                                                         <  1:2875638/91‑1 (MQ=255)
                    cTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaata                                                                        <  2:819775/91‑1 (MQ=255)
                      gAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataat                                                                      >  2:1197111/1‑91 (MQ=255)
                        aaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATga                                                                    <  2:725090/91‑2 (MQ=255)
                        aaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGt                                                                    <  1:1217277/91‑1 (MQ=255)
                        aaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGt                                                                    <  1:3224582/91‑1 (MQ=255)
                         aaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGtt                                                                   <  2:1460015/91‑1 (MQ=255)
                         aaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGtt                                                                   >  1:476452/1‑91 (MQ=255)
                             tCGCAAAGAACGTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACaa                                                               <  1:3514138/91‑1 (MQ=255)
                             tCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACaa                                                               >  2:1146087/1‑91 (MQ=255)
                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                              >  2:2391804/1‑91 (MQ=255)
                              cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTAATGCCCCAAAATATCTGAAATAATGTTACAat                                                              >  1:3325771/1‑91 (MQ=255)
                                cAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAatat                                                            <  1:2707833/91‑1 (MQ=255)
                                cAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAatat                                                            >  1:780262/1‑91 (MQ=255)
                                 aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                           <  1:3057730/91‑1 (MQ=255)
                                  aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                                          >  2:3857239/1‑91 (MQ=255)
                                    gAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGtt                                                        >  1:1827152/1‑91 (MQ=255)
                                         tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                                   >  2:2458768/1‑91 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATCTATGAAATAATGTTACAATATGGTTAACaaa                                                  <  1:3442984/91‑1 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  >  2:3097107/1‑91 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  <  2:240672/91‑1 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  <  1:1433234/91‑1 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  <  2:1480550/91‑1 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  <  2:1647698/91‑1 (MQ=255)
                                          aCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaaa                                                  >  1:1757263/1‑91 (MQ=255)
                                            gtGAACAATTGCATAAATAAATTACTGGGCGTGAGACTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  2:3840249/90‑1 (MQ=255)
                                            ctGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  2:1557293/90‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  1:3666001/91‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                >  1:734932/1‑91 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  2:396220/91‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  1:1673878/91‑1 (MQ=255)
                                            aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                                <  2:2830480/91‑1 (MQ=255)
                                             tGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGtt                                               >  1:1906479/1‑91 (MQ=255)
                                              gAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTa                                              <  2:3537619/91‑1 (MQ=255)
                                              gAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTa                                              >  2:3751021/1‑91 (MQ=255)
                                               aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                             <  1:898196/91‑1 (MQ=255)
                                                   aTTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttt                                         >  2:521816/1‑91 (MQ=255)
                                                   aTTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttt                                         >  2:1920337/1‑91 (MQ=255)
                                                     tGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttttt                                       >  2:2033690/1‑91 (MQ=255)
                                                      gCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTg                                      >  1:831193/1‑91 (MQ=255)
                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                                     >  1:3555549/1‑91 (MQ=255)
                                                       cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                                     >  2:1442001/1‑91 (MQ=255)
                                                        ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                                    >  1:3373377/1‑91 (MQ=255)
                                                         taaataaaTTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAAt                                   >  1:962408/1‑91 (MQ=255)
                                                         taaataaaTTACTGGGCGTGAGAATCATTCTCATGCCACAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAAt                                   >  1:3269808/1‑91 (MQ=255)
                                                               aatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGtt                                               <  1:3479342/73‑1 (MQ=255)
                                                               aatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTTTTTGTACGt                             >  2:3205005/1‑91 (MQ=255)
                                                                atTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATCATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                            >  1:3092977/1‑91 (MQ=255)
                                                                atTACTGGGCATGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTg                            <  2:1727196/91‑1 (MQ=255)
                                                                 ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg                           >  1:463353/1‑91 (MQ=255)
                                                                   aCTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGTATGTTTGTGGaa                         >  2:3310011/1‑91 (MQ=255)
                                                                   aCTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa                         >  2:130428/1‑91 (MQ=255)
                                                                   aCTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa                         <  1:1621353/91‑1 (MQ=255)
                                                                     tGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAgg                       <  1:2197701/91‑2 (MQ=255)
                                                                     tGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc                       >  2:2104615/1‑91 (MQ=255)
                                                                     ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc                       >  1:3382127/2‑91 (MQ=38)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                     <  1:3210225/91‑1 (MQ=255)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                     <  2:268744/91‑1 (MQ=255)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                     <  1:2544830/91‑1 (MQ=255)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                     <  2:67145/91‑1 (MQ=255)
                                                                       ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt                     <  1:192887/91‑1 (MQ=255)
                                                                         cGTGAGAATCAGTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGttt                   <  1:3519591/91‑1 (MQ=255)
                                                                             agaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTa               >  1:2450216/1‑91 (MQ=255)
                                                                             agaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTa               <  2:113527/91‑1 (MQ=255)
                                                                              gaATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaa              >  2:2752078/1‑91 (MQ=255)
                                                                              gaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaa              <  2:3373377/91‑1 (MQ=255)
                                                                               aatgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa             >  1:281484/5‑91 (MQ=255)
                                                                               aatgaTTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa             >  2:1443727/5‑91 (MQ=255)
                                                                               aaTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa             <  2:2156184/91‑1 (MQ=255)
                                                                               aaTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa             <  1:2954415/91‑1 (MQ=255)
                                                                               aaTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaaa             <  2:3770545/91‑1 (MQ=255)
                                                                                aTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAAt            >  2:2407454/1‑91 (MQ=255)
                                                                                   aTTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgag         >  1:2435315/1‑91 (MQ=255)
                                                                                   aTTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATgag         >  1:2661682/1‑91 (MQ=255)
                                                                                     tctcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGaa       >  2:1107117/1‑91 (MQ=255)
                                                                                      ctcATGCCCCAGAAGATATGAAAGCATGTTACAATATGGTTAAGAAAGTTTTTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGaaa      <  2:3603668/91‑1 (MQ=255)
                                                                                      ctcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGaaa      >  1:2015056/1‑91 (MQ=255)
                                                                                       tcATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAg     >  1:937214/1‑91 (MQ=255)
                                                                                        cATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGt    >  2:339436/1‑91 (MQ=255)
                                                                                        cATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGt    <  2:1410534/91‑1 (MQ=255)
                                                                                        cATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGt    <  1:432571/91‑1 (MQ=255)
                                                                                        cATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGt    >  1:2343440/1‑91 (MQ=255)
                                                                                         aTGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGtt   <  1:1854263/91‑1 (MQ=255)
                                                                                          tGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGttt  >  2:1616269/1‑91 (MQ=255)
                                                                                          tGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGttt  >  1:2950167/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGAGAAAGTTT  >  USA300TCH1516_ALE/1158353‑1158533

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