Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,292,894 A→C L24R (CTG→CGG)  dsbE ← periplasmic thioredoxin of cytochrome c‑type biogenesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,292,8940AC98.4% 425.8 / ‑2.2 126L24R (CTG→CGG) dsbEperiplasmic thioredoxin of cytochrome c‑type biogenesis
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (60/64);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 127 total aligned reads displayed.

GGATTGTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCAGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTGCGG  >  NC_000913/2292798‑2292990
                                                                                                |                                                                                                
ggATTGTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCGATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCt                                                                                               >  2:3030249/1‑100 (MQ=255)
     gTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCcac                                                                                          >  1:61785/1‑100 (MQ=255)
     gTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCcac                                                                                          <  2:1214907/100‑1 (MQ=255)
     gTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCcac                                                                                          >  1:3069534/1‑100 (MQ=255)
      tCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCcaca                                                                                         >  1:832883/1‑100 (MQ=255)
       ccAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcag                                                                                        <  1:424806/100‑1 (MQ=255)
       ccAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcag                                                                                        <  2:2440962/100‑1 (MQ=255)
       ccAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcag                                                                                        >  2:2537273/1‑100 (MQ=255)
        cAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcagc                                                                                       <  2:1679440/100‑1 (MQ=255)
         aGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcagca                                                                                      >  1:700887/1‑100 (MQ=255)
         aGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcagca                                                                                      >  2:1918333/1‑100 (MQ=255)
         aGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCAcagca                                                                                      >  1:3063518/1‑100 (MQ=255)
         aGTGATTCGAGACGAAACNTCGGCACAGGCTGGCCAATGCGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTNCGCGCCNGCTGCcagaaga                                                                                      >  1:3272138/1‑95 (MQ=255)
            gATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAgcg                                                                                   >  1:432807/1‑100 (MQ=255)
                  agaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCgg                                                                             >  2:2604559/1‑100 (MQ=255)
                  agaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCgg                                                                             <  1:3082488/100‑1 (MQ=255)
                  agaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCgg                                                                             <  1:1867806/100‑1 (MQ=255)
                  agaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCgg                                                                             >  2:1908873/1‑100 (MQ=255)
                   gaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGc                                                                            <  1:1177033/100‑1 (MQ=255)
                   gaCGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGc                                                                            >  2:1732163/1‑100 (MQ=255)
                      gAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAAt                                                                         >  1:1019616/1‑100 (MQ=255)
                        aaCTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCg                                                                       <  1:3279450/100‑1 (MQ=255)
                         aCTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGc                                                                      <  1:1335592/100‑1 (MQ=255)
                         aCTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGc                                                                      <  1:1337228/100‑1 (MQ=255)
                           cgcGGGACAGGCTGGCCAGTGAGCGCCGCTTGCAGATTGGTCGGAGCATCCCCGCGGTCATTACGCGGCCGTTGCCACGGCAGCGCCGCGGCAATCGCCa                                                                    <  2:2155151/98‑1 (MQ=255)
                            tCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAg                                                                   >  1:1884616/1‑100 (MQ=255)
                             cGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAgg                                                                  <  2:2229089/100‑1 (MQ=255)
                             cGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAgg                                                                  <  1:1295547/100‑1 (MQ=255)
                             cGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAgg                                                                  <  1:2140324/100‑1 (MQ=255)
                             cGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAgg                                                                  <  2:2690806/100‑1 (MQ=255)
                             cGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAgg                                                                  <  1:3658656/100‑1 (MQ=255)
                                cacaGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAg                                                               >  1:1892161/1‑100 (MQ=255)
                                 acaGGCTTGCCAATGAGCGCCGATTCCAGCTTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGa                                                              >  1:2958267/1‑100 (MQ=255)
                                   aGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATa                                                            <  2:2413311/100‑1 (MQ=255)
                                    ggCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATaa                                                           <  1:2544921/100‑1 (MQ=255)
                                    ggCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATaa                                                           >  1:3404610/1‑100 (MQ=255)
                                    ggCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATaa                                                           <  1:2357322/100‑1 (MQ=255)
                                    ggCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATaa                                                           <  1:241338/100‑1 (MQ=255)
                                    ggCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATaa                                                           <  2:176407/100‑1 (MQ=255)
                                     gCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAAt                                                          >  1:3590660/1‑100 (MQ=255)
                                      cTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATc                                                         <  1:846207/100‑1 (MQ=255)
                                         gCCAATGAGCGCCGATTCCAGACTGTTCGTATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCCGGAAGATAATCAAc                                                      >  2:1353601/1‑100 (MQ=255)
                                            aaTGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGa                                                   <  2:1625812/100‑1 (MQ=255)
                                              tGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAAt                                                 <  2:1193803/100‑1 (MQ=255)
                                                  cgcTGATTCCAGATGGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCGGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAAc                                             <  2:1775037/100‑1 (MQ=255)
                                                  cgcCGATTCCGGAGGGGCCGGACCACCCCCTGGGGCATTCCGCGCCCGCGGCCACAGCAGCGCCGGGGCAATCGCCAGGAAGATATTCAACGGAATTacc                                             <  2:1338358/100‑3 (MQ=255)
                                                   gcCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCGACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACa                                            <  2:1112930/100‑1 (MQ=255)
                                                   gcCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACa                                            >  1:183946/1‑100 (MQ=255)
                                                   gcCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACa                                            >  1:1512326/1‑100 (MQ=255)
                                                    ccGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACaa                                           >  2:226837/1‑100 (MQ=255)
                                                     cGATTCCAGATTGGTCGCATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAAt                                          <  1:3550037/100‑1 (MQ=255)
                                                      gATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATa                                         <  1:2405940/100‑1 (MQ=255)
                                                         tCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACtt                                      <  2:3633225/100‑1 (MQ=255)
                                                          ccAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACttt                                     >  1:1073208/1‑100 (MQ=255)
                                                          ccAGATTGGTCGGATCATCCCCTTCGGCATTACGCACCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACttt                                     >  1:75456/1‑100 (MQ=255)
                                                           cAGATTGGTCGGATCATCCCATTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTg                                    >  2:1805331/1‑100 (MQ=255)
                                                             gATTGGTCGGATCATCCCCTTCGGCATTACGCTCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCATGAAGATAATCAACGGAATTAACAATACTTTgcg                                  >  2:1091869/1‑100 (MQ=255)
                                                             gATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTgcg                                  <  1:2389413/100‑1 (MQ=255)
                                                             gATTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGCAATTAACAATACTTTgcg                                  >  1:3827182/1‑100 (MQ=255)
                                                              aTTGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTgcgc                                 <  1:3809928/100‑1 (MQ=255)
                                                                tGGTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCtt                               >  1:2164857/1‑100 (MQ=255)
                                                                 ggTGGGATCATCCCCTTTGGCCTTACGCGCCCGCAGCCACAGAAGCGCCGCGGCAATCCCCAGGAAAATAATCTACGGAATTAACAATACTTTGCGCtac                              >  2:3474482/1‑98 (MQ=255)
                                                                 ggTCGGATCATCCCCTTCGGCCTTATGCGCCCGCTGACACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTc                              <  2:3788647/100‑1 (MQ=255)
                                                                 ggTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTc                              <  1:1800850/100‑1 (MQ=255)
                                                                 ggTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTc                              >  2:1608938/1‑100 (MQ=255)
                                                                 ggTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTc                              <  1:2438626/100‑1 (MQ=255)
                                                                  gTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCa                             >  1:3081563/1‑100 (MQ=255)
                                                                  gTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCa                             <  2:1478106/100‑1 (MQ=255)
                                                                  gTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCa                             >  1:2768365/1‑100 (MQ=255)
                                                                  gTCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCa                             >  1:256574/1‑100 (MQ=255)
                                                                   tccgctcatcCCCTGCGGCATGCCGCGCCCGCTTCTACAGCAGCCTCTCGGCAATCGCCAGGAAGATAATCAACGTAATTAACAATACTTTGCGCTTCAt                            <  2:3369011/95‑1 (MQ=255)
                                                                   tCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCAt                            >  2:3922511/1‑100 (MQ=255)
                                                                   tCGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCAt                            >  2:671378/1‑100 (MQ=255)
                                                                    cGGATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATa                           >  2:804646/1‑100 (MQ=255)
                                                                     ggATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGCATTAACACTACTTTGCGCTTCATAc                          >  2:1750799/1‑100 (MQ=255)
                                                                     ggATCATCCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATAc                          <  1:2250425/100‑1 (MQ=255)
                                                                       atcatcCCCTTCGGCATTACGCGCCAGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACgg                        <  1:1636818/100‑1 (MQ=255)
                                                                        tcatcCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGc                       <  1:1620842/100‑1 (MQ=255)
                                                                        tcatcCCCTTCGCCATTACCCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGTAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGc                       <  2:3323416/100‑1 (MQ=255)
                                                                          atcCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCt                     >  1:3781735/1‑100 (MQ=255)
                                                                           tcCCCTTCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTc                    >  1:1550212/1‑100 (MQ=255)
                                                                               cTTCGGCATTACGCGCCCGCTGCAAAAGCAGCGCCGCGGCAATCGCCAGGAAGATCATCAACGGAATTAACAATCCTTTGCGCTTCATACGGCCTCCGgc                <  2:3023895/100‑1 (MQ=255)
                                                                                ttCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGgcg               <  2:611728/100‑1 (MQ=255)
                                                                                ttCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGgcg               >  2:1211563/1‑100 (MQ=255)
                                                                                 tCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGgcgc              >  1:1276707/1‑100 (MQ=255)
                                                                                 tCGGCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACAGCCTCCGgcgc              >  1:1140922/1‑100 (MQ=255)
                                                                                    gCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGt           >  2:1439956/1‑100 (MQ=255)
                                                                                    gCATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGt           <  1:1493866/100‑1 (MQ=255)
                                                                                     cATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGtt          <  1:2359013/100‑1 (MQ=255)
                                                                                     cATTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGtt          >  1:488746/1‑100 (MQ=255)
                                                                                      aTTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGttt         <  1:3699882/100‑1 (MQ=255)
                                                                                      aTTACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGttt         >  1:2362351/1‑100 (MQ=255)
                                                                                        tACGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGttttt       <  2:1955255/100‑1 (MQ=255)
                                                                                         aCGCGCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGtttttt      >  1:2194689/1‑100 (MQ=255)
                                                                                          cgcgCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTg     <  2:3100537/100‑1 (MQ=255)
                                                                                          cgcgCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTg     >  1:1767560/1‑100 (MQ=255)
                                                                                          cgcgCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTg     <  1:3830975/100‑1 (MQ=255)
                                                                                          cgcgCCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAAGTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTg     <  1:629885/100‑1 (MQ=255)
                                                                                           gcgcCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGCAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTGc    <  2:3634395/100‑1 (MQ=255)
                                                                                             gcCCGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTGCgg  >  2:2147188/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGATTGTCCAGTGATTCGAGACGAAACTTCGGCACAGGCTTGCCAATGAGCGCCGATTCCAGATTGGTCGGATCATCCCCTTCGGCATTACGCGCCAGCTGCCACAGCAGCGCCGCGGCAATCGCCAGGAAGATAATCAACGGAATTAACAATACTTTGCGCTTCATACGGCCTCCGGCGCAGTTTTTTGCGG  >  NC_000913/2292798‑2292990

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