Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,034,326 A→G 57.2% N358S (AAC→AGC)  hyaB → hydrogenase 1, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,034,3260AG57.2% 85.8 / 164.5 147N358S (AAC→AGC) hyaBhydrogenase 1, large subunit
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base G (49/35);  total (80/67)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 150 total aligned reads displayed.

CGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCTGAATGA  >  NC_000913/1034227‑1034425
                                                                                                   |                                                                                                   
cgGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCaa                                                                                                     <  1:2308000/100‑1 (MQ=255)
    gATTAACGGCGACTTCAACACTGTGCTGCCAGTGGATTGGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCaa                                                                                                     <  2:2175499/96‑1 (MQ=255)
      ttAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCa                                                                                               <  1:1754266/100‑1 (MQ=255)
      ttAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATcn                                                                                               <  1:2988407/100‑2 (MQ=255)
      ttAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCa                                                                                               <  1:3523510/100‑1 (MQ=255)
      ttAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCa                                                                                               <  2:2711331/100‑1 (MQ=255)
        aaCGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAgg                                                                                             >  1:1285882/1‑100 (MQ=255)
         aCGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGt                                                                                            >  2:1534421/1‑100 (MQ=255)
         aCGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTCGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCCGCGATCAGGt                                                                                            >  2:1618528/1‑100 (MQ=255)
          cGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTc                                                                                           >  1:2933979/1‑100 (MQ=255)
               aCTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCg                                                                                      <  2:3153688/100‑1 (MQ=255)
               aCTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCg                                                                                      <  1:1437354/100‑1 (MQ=255)
               aCTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCg                                                                                      <  2:1738473/100‑1 (MQ=255)
               aCTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCg                                                                                      <  1:2533394/100‑1 (MQ=255)
                cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGt                                                                                     >  2:885831/1‑100 (MQ=255)
                cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGt                                                                                     >  2:3569796/1‑100 (MQ=255)
                cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGt                                                                                     >  2:2982245/1‑100 (MQ=255)
                cTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGt                                                                                     >  1:356805/1‑100 (MQ=255)
                  tCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCa                                                                                   >  1:2099187/1‑100 (MQ=255)
                     acaaTGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATcc                                                                                <  1:407906/100‑1 (MQ=255)
                         tgtgCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTc                                                                            >  2:3337181/1‑100 (MQ=255)
                         tgtgCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTCTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTc                                                                            <  1:3386882/100‑1 (MQ=255)
                             ctgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATg                                                                        <  2:213636/100‑1 (MQ=255)
                             ctgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATg                                                                        <  1:129030/100‑1 (MQ=255)
                              tgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATgg                                                                       <  1:1912097/100‑1 (MQ=255)
                              tgcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATgg                                                                       <  2:3829425/100‑1 (MQ=255)
                              tgcCAGTGGACGTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGAt                                                                         <  2:3149351/98‑1 (MQ=255)
                               gcCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGc                                                                      >  1:1231577/1‑100 (MQ=255)
                                ccAGTGGATTTGGTTGATCCGCAGCATGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCa                                                                     <  1:584821/100‑1 (MQ=255)
                                ccAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCa                                                                     >  2:1498598/1‑100 (MQ=255)
                                 cAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCAt                                                                    <  2:4053711/100‑1 (MQ=255)
                                  aGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATc                                                                   >  2:3121965/1‑100 (MQ=255)
                                  aGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATc                                                                   <  1:3412117/100‑1 (MQ=255)
                                   gTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCa                                                                  <  1:305451/100‑1 (MQ=255)
                                     ggATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCAcc                                                                <  2:149316/100‑1 (MQ=255)
                                     ggATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCAcc                                                                <  1:2498593/100‑1 (MQ=255)
                                     ggATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCAcc                                                                <  2:149315/100‑1 (MQ=255)
                                      gATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCg                                                               >  2:2619172/1‑100 (MQ=255)
                                         ttGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGAcc                                                            <  1:2619188/100‑1 (MQ=255)
                                         ttGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGAcc                                                            <  1:45783/100‑1 (MQ=255)
                                           ggTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCg                                                          >  1:867866/1‑100 (MQ=255)
                                           ggTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCg                                                          >  1:2697591/1‑100 (MQ=255)
                                             ttGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTg                                                        >  2:837886/1‑100 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCgggg                                                       >  1:51805/1‑100 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTgg                                                       >  1:3584067/1‑100 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTgg                                                       >  2:3885422/1‑100 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTgg                                                       >  1:1298623/1‑100 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTgg                                                       <  1:3062603/100‑1 (MQ=255)
                                              tGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTgg                                                       >  1:1691579/1‑100 (MQ=255)
                                               gATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTCTCGATATCCCAGCGATCAGGTCGGCCGTCATCCGTTCGATGGCATCACCGCCCCGTGGt                                                      >  2:1013132/1‑100 (MQ=255)
                                               gATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGt                                                      <  1:279063/100‑1 (MQ=255)
                                                aTCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTa                                                     <  1:2586542/100‑1 (MQ=255)
                                                aTCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTa                                                     >  1:3910810/1‑100 (MQ=255)
                                                 tCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTAc                                                    >  1:1409343/1‑100 (MQ=255)
                                                   cGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACaa                                                  >  2:3361012/1‑100 (MQ=255)
                                                   cGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGAGATCAAGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACaa                                                  >  2:3765882/1‑100 (MQ=255)
                                                    gcagcaGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAAc                                                 >  2:3749187/1‑100 (MQ=255)
                                                     cagcagGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAAcc                                                <  2:2753792/100‑1 (MQ=255)
                                                       gcagGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAAcccc                                              >  1:3616739/1‑100 (MQ=255)
                                                          ggTGCAGGGGTTTGTCGTCCCCGCCTGGTATCGATATCCCAACGATCAGGGCGGGCGTCATCCGTTCGACGGCATCACCGACCCGTGGTACAACCCCGGc                                           <  2:3261810/100‑1 (MQ=255)
                                                            tGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGa                                         >  2:687280/1‑100 (MQ=255)
                                                             gCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGAt                                        >  2:3028451/1‑100 (MQ=255)
                                                              cAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATg                                       <  1:2198275/100‑1 (MQ=255)
                                                              cAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATg                                       >  2:3463359/1‑100 (MQ=255)
                                                               aGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGt                                      >  1:2183821/1‑100 (MQ=255)
                                                               aGGAGTTTGTCGACCACGCATGGTATCGATATAGCAGCGATGACATCGGGCGTCATCCGTTCGATGGCCACACCGACCCGTGGTACAACCCCGGCgtgg                                       >  2:4128673/1‑96 (MQ=255)
                                                                agAGTTTGTCGACCACGCCGGGTATCGAGATCCCAACGATCAGGTCGGGCGTCAGCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGAt                                        <  2:2851026/96‑1 (MQ=255)
                                                                 gAGTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCa                                    <  1:3121984/100‑1 (MQ=255)
                                                                  aGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCaa                                   >  1:2970031/1‑100 (MQ=255)
                                                                   gTTTGTCGACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCaaa                                  >  1:3953570/1‑100 (MQ=255)
                                                                    tttGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAg                                 <  2:1177287/100‑1 (MQ=255)
                                                                       gTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCa                              >  2:1433988/1‑100 (MQ=255)
                                                                       ccgatctaCGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCa                              <  1:2087694/93‑1 (MQ=255)
                                                                          gACCACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCg                           <  1:3740837/100‑1 (MQ=255)
                                                                          gACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCg                           >  2:2573391/1‑100 (MQ=255)
                                                                          gACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCg                           >  2:1385205/1‑100 (MQ=255)
                                                                            ccACGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGAt                         <  1:837892/100‑1 (MQ=255)
                                                                             cACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATa                        <  1:3297146/100‑1 (MQ=255)
                                                                              aCGCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCaagatcgg                       >  2:2087681/1‑93 (MQ=255)
                                                                              aCGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATAc                       >  1:2663969/1‑100 (MQ=255)
                                                                              aCGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATAc                       <  2:711517/100‑1 (MQ=255)
                                                                                gCCTGGTATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCa                     >  1:1863119/1‑100 (MQ=255)
                                                                                gCCTGGTATCGATACCCCAGCGATCAGGGCGGGCGTCAGCCGTTCGATGGCATCACCGACCCGTCGCACAACCCCGGCGATGTCAAAGGCAGCGCTAcaa                     >  2:1183852/1‑98 (MQ=255)
                                                                                 ccTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCaa                    >  1:1747606/1‑100 (MQ=255)
                                                                                  cgggTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAAc                   <  1:624796/98‑1 (MQ=255)
                                                                                     gTATCGATATCCAAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACGCGTGGTACAACCCCGGCGATGTCAAAGGCAGAGATACCAACa                  >  2:2040459/1‑98 (MQ=255)
                                                                                      tATCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTc               >  1:4099751/1‑100 (MQ=255)
                                                                                      tATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTc               <  2:101718/100‑1 (MQ=255)
                                                                                      tATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTc               <  1:3749211/100‑1 (MQ=255)
                                                                                        tCGATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTcag             >  2:3838115/1‑100 (MQ=255)
                                                                                          gATATCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTcagc            >  2:2926705/1‑99 (MQ=255)
                                                                                            tatCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCagcagc         <  2:516760/100‑1 (MQ=255)
                                                                                             atCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCt        >  1:1486418/1‑100 (MQ=255)
                                                                                              tCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtg       <  1:2310658/100‑1 (MQ=255)
                                                                                              tCCCAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtg       <  1:3628453/100‑1 (MQ=255)
                                                                                              tCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtg       <  1:1006540/100‑1 (MQ=255)
                                                                                                tcAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtgaa     <  1:3138154/99‑1 (MQ=255)
                                                                                                ccAGCGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtgaa     >  2:2473375/1‑100 (MQ=255)
                                                                                                ccAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCtgaa     >  1:15832/1‑100 (MQ=255)
                                                                                                 cAACGATCAGGTCGGGCGTCATGCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGCTACCAACATTCAGCAGCtgaat    <  2:2402839/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGGTGATTAACGGCGACTTCAACAATGTGCTGCCAGTGGATTTGGTTGATCCGCAGCAGGTGCAGGAGTTTGTCGACCACGCCTGGTATCGATATCCCAACGATCAGGTCGGGCGTCATCCGTTCGATGGCATCACCGACCCGTGGTACAACCCCGGCGATGTCAAAGGCAGCGATACCAACATTCAGCAGCTGAATGA  >  NC_000913/1034227‑1034425

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: