Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,158,435 C→G 9.6% intergenic (+19/‑305) 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,4350CG 9.6% 750.1 / 53.4 228intergenic (+19/‑305)uvrC/sdhCUvrABC system protein C/Succinate dehydrogenase cytochrome b558 subunit
Reads supporting (aligned to +/- strand):  ref base C (106/100);  new base G (7/15);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 234 total aligned reads displayed.

CTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGA  >  USA300TCH1516_ALE/1158345‑1158525
                                                                                          |                                                                                          
cTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATc                                                                                            >  1:1126907/1‑91 (MQ=255)
cTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATc                                                                                            >  2:3830005/1‑91 (MQ=255)
cTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATc                                                                                            <  2:1948464/91‑1 (MQ=255)
  tGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTGGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATGAt                                                                                          <  2:3112300/91‑1 (MQ=255)
   gAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGAtt                                                                                         <  2:647649/91‑1 (MQ=255)
   gAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGAtt                                                                                         <  1:557474/91‑1 (MQ=255)
    aaGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATGATtc                                                                                        <  1:3288185/91‑1 (MQ=255)
    aaGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCATGAGAATGATtc                                                                                        <  2:2986373/91‑1 (MQ=255)
     aGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATtct                                                                                       <  1:2416606/91‑1 (MQ=255)
       aTTTTAAAAATATAGGTATTCCTGACACCGTCGCAAAGAACCCACATGAACAATTGCATCAATAAATTACTGGGCGTGAGActcctcctcc                                                                                     >  2:3743756/1‑90 (MQ=255)
         tttAAAAATATTGGTATTCCTGAAAACGTCGCAAAGAACCTTCATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATg                                                                                   <  1:408798/91‑1 (MQ=255)
          ttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGc                                                                                  >  1:3770176/1‑91 (MQ=255)
          ttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGc                                                                                  <  1:1946672/91‑1 (MQ=255)
          ttAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGc                                                                                  >  1:2825097/1‑91 (MQ=255)
          ggaaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATGGCATAAATAAATTACTGGGCATGAGAATGATTCTCACGc                                                                                  <  2:20122/89‑1 (MQ=255)
            aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGATCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                <  2:765763/91‑1 (MQ=255)
            aaaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGccc                                                                                >  1:1744087/1‑91 (MQ=255)
             aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                               <  1:993088/91‑1 (MQ=255)
             aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                               <  2:1500696/91‑1 (MQ=255)
             aaaaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGcccc                                                                               >  1:3831665/1‑91 (MQ=255)
               aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaa                                                                             <  1:833195/91‑1 (MQ=255)
               aaTATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaa                                                                             >  1:1383629/1‑91 (MQ=255)
                atatAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaa                                                                            >  2:3917654/1‑91 (MQ=255)
                 tataGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCaaaa                                                                           <  2:2589933/91‑1 (MQ=255)
                        aTTCCTGAAAACGTCGCAGAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                                    <  1:860066/91‑1 (MQ=255)
                        aTTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGa                                                                    <  2:2771631/91‑1 (MQ=255)
                         ttCCGGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGaa                                                                   <  1:3665485/91‑1 (MQ=255)
                           ccTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaat                                                                 >  2:311453/1‑91 (MQ=255)
                            cTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaata                                                                >  1:3905468/1‑91 (MQ=255)
                             tGAAATCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataa                                                               <  2:3644389/91‑1 (MQ=255)
                             tGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataa                                                               <  2:2639174/91‑1 (MQ=255)
                             tGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataa                                                               >  1:2915120/1‑91 (MQ=255)
                              gAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAaataat                                                              >  2:3333387/1‑91 (MQ=255)
                               aaaaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATCACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATg                                                             >  2:3052634/1‑91 (MQ=255)
                                 aaCGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGtt                                                           <  1:3232298/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)
                                      cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                      >  1:3109301/1‑91 (MQ=255)
                                      cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                      >  2:2557016/1‑91 (MQ=255)
                                      cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                      >  2:3713536/1‑91 (MQ=255)
                                      cGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAat                                                      >  2:860216/1‑91 (MQ=255)
                                       gCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAata                                                     >  2:3232116/1‑91 (MQ=255)
                                         aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                   <  2:3536040/91‑1 (MQ=255)
                                         aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                   <  2:3169410/91‑1 (MQ=255)
                                         aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                   <  2:1442583/91‑1 (MQ=255)
                                         aaaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATg                                                   <  1:1841734/91‑1 (MQ=255)
                                          aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                                  >  2:2298884/1‑91 (MQ=255)
                                          aaGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATgg                                                  >  2:1243659/1‑91 (MQ=255)
                                           aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                                 >  2:1665923/1‑91 (MQ=255)
                                           aGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGt                                                 <  2:2843159/91‑1 (MQ=255)
                                              aCCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTaa                                              >  1:39234/1‑91 (MQ=255)
                                               ccTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAAc                                             >  2:503545/1‑91 (MQ=255)
                                                cTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACa                                            >  2:628464/1‑91 (MQ=255)
                                                 tACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACaa                                           >  2:2392225/1‑91 (MQ=255)
                                                    aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                        <  2:2201330/91‑1 (MQ=255)
                                                    aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                        <  2:345523/91‑1 (MQ=255)
                                                    aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                        <  1:1003962/91‑1 (MQ=255)
                                                    aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                        <  2:2784554/91‑1 (MQ=255)
                                                    aTGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATGCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGt                                        <  2:1303583/91‑1 (MQ=255)
                                                     tGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGtt                                       >  2:1001970/1‑91 (MQ=255)
                                                      gAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTa                                      >  1:2647212/1‑91 (MQ=255)
                                                       aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                     <  2:169607/91‑1 (MQ=255)
                                                       aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                     >  1:2176290/1‑91 (MQ=255)
                                                       aaCAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAt                                     >  1:1614028/1‑91 (MQ=255)
                                                         cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttt                                   >  2:2925515/1‑91 (MQ=255)
                                                         cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttt                                   <  1:635784/91‑1 (MQ=255)
                                                         cAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttt                                   >  2:3599701/1‑91 (MQ=255)
                                                            ttGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttttt                                >  2:2798882/1‑91 (MQ=255)
                                                            ttGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAtttttt                                <  2:3368461/91‑1 (MQ=255)
                                                             tGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTAttttttt                               >  2:306594/1‑91 (MQ=255)
                                                               cATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGa                             >  2:1221051/1‑91 (MQ=255)
                                                                ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGaa                            >  2:1475604/1‑91 (MQ=255)
                                                                ataaataaATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTGTGaa                            <  1:3122542/91‑1 (MQ=255)
                                                                  aaataaatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATg                          <  2:1122000/91‑1 (MQ=255)
                                                                   aataaatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGt                         >  1:2859773/1‑91 (MQ=255)
                                                                     taaatTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTAc                       >  1:1158935/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)
                                                                         ttACTGGGCGTGAGAATCATTCTCATGCCCCAGAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTCTTTGAATGTACGTgg                   <  1:3367695/91‑1 (MQ=255)
                                                                         ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg                   <  1:2781663/91‑1 (MQ=255)
                                                                         ttACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTgg                   >  1:3602192/1‑91 (MQ=255)
                                                                           ttggggCATGAGAATGATTCTCACGCCCCAATATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTCGAATGTACGTAGaa                 >  1:3564139/4‑91 (MQ=25)
                                                                           ttggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa                 >  1:3286116/4‑91 (MQ=38)
                                                                           aCTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGaa                 <  2:2112858/91‑1 (MQ=255)
                                                                            cTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAg                >  2:1875376/1‑91 (MQ=255)
                                                                             tGGGCGTGAGAATCATTCTCATGTCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc               >  2:1117260/1‑91 (MQ=255)
                                                                             ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc               >  1:765977/2‑91 (MQ=38)
                                                                             ggggCATGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGc               >  2:1750655/2‑91 (MQ=38)
                                                                              gggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCg              <  1:1611440/91‑1 (MQ=255)
                                                                              gggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCg              <  1:1431779/91‑1 (MQ=255)
                                                                               ggCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGt             <  2:1011617/91‑1 (MQ=255)
                                                                                 cGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGttt           <  1:2998447/91‑1 (MQ=255)
                                                                                   tGAGAATGATTCTCACGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTc         >  1:2735648/1‑91 (MQ=255)
                                                                                      gaATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTaa      <  2:966542/91‑1 (MQ=255)
                                                                                        aTCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTCAACAAAGTTATTTTTTTGAATGTACGTGGCAGCGTTTTCTAAAt    >  2:3399303/1‑91 (MQ=255)
                                                                                         tCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATg   >  1:593146/1‑91 (MQ=255)
                                                                                         tCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATg   <  1:14474/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTTGAAGATTTTAAAAATATAGGTATTCCTGAAAACGTCGCAAAGAACCTACATGAACAATTGCATAAATAAATTACTGGGCGTGAGAATCATTCTCATGCCCCAAAATATATGAAATAATGTTACAATATGGTTAACAAAGTTATTTTTTTGAATGTACGTGGAAGCGTTTTCTAAATGA  >  USA300TCH1516_ALE/1158345‑1158525

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: