Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,612,973 C→G 6.2% L358L (CTC→CTG hyfR → DNA‑binding transcriptional activator HyfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,612,9730CG 6.2% 323.5 / 9.6 96L358L (CTC→CTGhyfRDNA‑binding transcriptional activator HyfR
Reads supporting (aligned to +/- strand):  ref base C (43/47);  new base G (4/2);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

TGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCA  >  NC_000913/2612826‑2613118
                                                                                                                                                   |                                                                                                                                                 
tGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCca                                                                                                                                                 >  1:1516102/1‑150 (MQ=255)
   cGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTGgagca                                                                                                                                              <  2:331509/150‑1 (MQ=255)
       ggCAATGCCGATGCCTGGCGTAGCACGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTa                                                                                                                                          <  1:1109154/150‑1 (MQ=255)
        gCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAg                                                                                                                                         <  1:38671/150‑1 (MQ=255)
         cAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGa                                                                                                                                        <  1:1317636/150‑1 (MQ=255)
          aaTGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGat                                                                                                                                       >  1:1289310/1‑150 (MQ=255)
          aaTGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGat                                                                                                                                       <  2:609771/150‑1 (MQ=255)
            tGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGatat                                                                                                                                     <  1:748380/150‑1 (MQ=255)
                gATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATtgtg                                                                                                                                 >  1:1423853/1‑150 (MQ=255)
                gATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATtgtg                                                                                                                                 >  1:870451/1‑150 (MQ=255)
                 aTGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTgg                                                                                                                                <  2:676573/150‑1 (MQ=255)
                 aTGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTgg                                                                                                                                <  1:1754236/150‑1 (MQ=255)
                    ccTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGa                                                                                                                             >  1:269036/1‑150 (MQ=255)
                    ccTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATAGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGa                                                                                                                             >  1:475931/1‑150 (MQ=255)
                     cTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGaa                                                                                                                            >  2:1502919/1‑150 (MQ=255)
                          gTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCg                                                                                                                       <  2:810247/150‑1 (MQ=255)
                          gTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCg                                                                                                                       <  2:1204361/150‑1 (MQ=255)
                             gCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACa                                                                                                                    <  1:243155/150‑1 (MQ=255)
                              cATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGGAGCTCCTTTCGAATATGGGCATCGGTGACATTATCTATCAAAGACAGGCAATGGAAGACCTGCTCCAGAAGGTAGATATTGTGGCGAAGAGCTACAg                                                                                                                   <  2:743574/150‑1 (MQ=255)
                              cATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAg                                                                                                                   >  2:1494209/1‑150 (MQ=255)
                                tGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTa                                                                                                                 >  1:243612/1‑150 (MQ=255)
                                 gACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTAc                                                                                                                >  2:1652248/1‑150 (MQ=255)
                                 gACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTAc                                                                                                                >  2:773283/1‑150 (MQ=255)
                                  aCCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACg                                                                                                               >  2:1205377/1‑150 (MQ=255)
                                   ccGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACgg                                                                                                              <  2:1423853/150‑1 (MQ=255)
                                        ttGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTg                                                                                                         <  2:870451/150‑1 (MQ=255)
                                         tGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGa                                                                                                        <  1:773283/150‑1 (MQ=255)
                                          gCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGAt                                                                                                       >  2:1092154/1‑150 (MQ=255)
                                           cAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGAtt                                                                                                      >  1:1084968/1‑150 (MQ=255)
                                            aGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGAttt                                                                                                     <  1:1753994/150‑1 (MQ=255)
                                             ggAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTtg                                                                                                    >  1:881065/1‑150 (MQ=255)
                                                aaGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtgg                                                                                                 <  1:857042/150‑1 (MQ=255)
                                                aaGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtgg                                                                                                 <  1:1091146/150‑1 (MQ=255)
                                                aaGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtgg                                                                                                 <  2:1516102/150‑1 (MQ=255)
                                                aaGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtgg                                                                                                 <  1:1695963/150‑1 (MQ=255)
                                                aaGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtgg                                                                                                 <  1:187090/150‑1 (MQ=255)
                                                 aGTTTGCAGCAAGAAAACCACCAGCTTTGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGAAATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGGGACAGTACGGTGTTGATTTgtggt                                                                                                <  2:1846894/150‑1 (MQ=255)
                                                  gTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTgtggtg                                                                                               <  2:1756988/150‑1 (MQ=255)
                                                      gcagcaAGAAAACCACCAGCTTAGCGAGCAGCACCTTTCGAATCTGGGCATCCGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAAc                                                                                           >  2:1876295/1‑150 (MQ=255)
                                                         gcaAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCgg                                                                                        >  1:506949/1‑150 (MQ=255)
                                                          caAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGa                                                                                       >  2:1935420/1‑150 (MQ=255)
                                                           aaGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCGTCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGaa                                                                                      >  2:674186/1‑150 (MQ=255)
                                                           aaGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGaa                                                                                      >  2:1447429/1‑150 (MQ=255)
                                                             gAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCGTCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACt                                                                                    <  1:674186/150‑1 (MQ=255)
                                                               aaaCCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTgg                                                                                  <  2:1289310/150‑1 (MQ=255)
                                                                 accaccAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCa                                                                                <  1:228827/150‑1 (MQ=255)
                                                                 accaccAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCa                                                                                >  1:953435/1‑150 (MQ=255)
                                                                 accaccAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCa                                                                                <  2:471258/150‑1 (MQ=255)
                                                                   caccaGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCaaa                                                                              <  2:243612/150‑1 (MQ=255)
                                                                      caGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAagag                                                                           >  1:433117/1‑150 (MQ=255)
                                                                       aGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAgg                                                                          <  1:1531715/150‑1 (MQ=255)
                                                                           tAGCGAGCAGCTCCTTTCTAATCTGGGCATCGGGTACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGAt                                                                      <  2:881065/150‑1 (MQ=255)
                                                                              cGAGCAGCTCATTTAGAATCTGGGAATCGGTGACATTATCTATAAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGATGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGc                                                                   >  2:1708491/1‑150 (MQ=255)
                                                                               gAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGcc                                                                  <  1:1175185/150‑1 (MQ=255)
                                                                               gAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGcc                                                                  <  2:1206998/150‑1 (MQ=255)
                                                                                agcagcTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCa                                                                 >  2:576403/1‑150 (MQ=255)
                                                                                    gcTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGc                                                             <  1:1944980/150‑1 (MQ=255)
                                                                                     cTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCg                                                            >  1:1358160/1‑150 (MQ=255)
                                                                                     cTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCg                                                            >  1:1487093/1‑150 (MQ=255)
                                                                                      tCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGa                                                           <  1:423357/150‑1 (MQ=255)
                                                                                       ccTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCTACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGAt                                                          <  1:828245/150‑1 (MQ=255)
                                                                                        cTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATc                                                         >  1:1829466/1‑150 (MQ=255)
                                                                                        cTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATc                                                         >  1:30859/1‑150 (MQ=255)
                                                                                         tttCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATcc                                                        >  2:414127/1‑150 (MQ=255)
                                                                                         tttCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATcc                                                        <  1:1374053/150‑1 (MQ=255)
                                                                                          ttCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCa                                                       <  2:1084968/150‑1 (MQ=255)
                                                                                          ttCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCa                                                       <  1:1935420/150‑1 (MQ=255)
                                                                                           tCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCAt                                                      <  2:981703/150‑1 (MQ=255)
                                                                                             gAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCa                                                    <  1:1779064/150‑1 (MQ=255)
                                                                                              aaTCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCaa                                                   <  1:223476/150‑1 (MQ=255)
                                                                                                tCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACt                                                 <  1:1162736/150‑1 (MQ=255)
                                                                                                  tGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTa                                               >  1:841023/1‑150 (MQ=255)
                                                                                                   gggCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAg                                              >  2:252235/1‑150 (MQ=255)
                                                                                                          ggTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCcgcga                                       >  1:1971165/1‑150 (MQ=255)
                                                                                                           gTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCcgcgac                                      >  2:665987/1‑150 (MQ=255)
                                                                                                                aTTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCgcgacgcgac                                 <  2:1582915/150‑1 (MQ=255)
                                                                                                                  tatctatcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACaa                               >  2:1958859/1‑150 (MQ=255)
                                                                                                                   atctatcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAg                              <  2:433117/150‑1 (MQ=255)
                                                                                                                    tctatcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGc                             <  1:1652248/150‑1 (MQ=255)
                                                                                                                      tatcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCg                           >  2:19809/1‑150 (MQ=255)
                                                                                                                       atcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGc                          >  1:578780/1‑150 (MQ=255)
                                                                                                                       atcAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGc                          <  2:1829466/150‑1 (MQ=255)
                                                                                                                            aGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTc                     >  2:913822/1‑150 (MQ=255)
                                                                                                                               cAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCaaa                  >  2:1673664/1‑150 (MQ=255)
                                                                                                                               cAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCaaa                  >  2:1300968/1‑150 (MQ=255)
                                                                                                                                 ggCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAAt                <  2:1738240/150‑1 (MQ=255)
                                                                                                                                        atctaCCTGCTGGAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTgc         >  2:291859/5‑150 (MQ=255)
                                                                                                                                        atctaCCTGCTGGAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACAGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGAGACGAGACAAGCCGCTGGTCAAAATCAACTgc         >  1:975736/5‑150 (MQ=255)
                                                                                                                                         tctaCCTGCTGGAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTgcg        >  2:456449/4‑150 (MQ=255)
                                                                                                                                          ctaCCTGCTGGAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCcg                                          <  2:992087/113‑1 (MQ=255)
                                                                                                                                          ctaCCTGCTGGAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCcg                                          >  1:992087/3‑115 (MQ=255)
                                                                                                                                          aGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTgcgc       <  1:414127/150‑1 (MQ=255)
                                                                                                                                             ccTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCCTCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGc    >  2:863276/1‑150 (MQ=255)
                                                                                                                                             ccTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGc    >  1:81697/1‑150 (MQ=255)
                                                                                                                                              cTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGcc   >  2:434429/1‑150 (MQ=255)
                                                                                                                                               tGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCa  >  2:1064928/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCCGTTGGCAATGCCGATGCCTGGCGTAGCATGACCGATTTGCAGGAAAGTTTGCAGCAAGAAAACCACCAGCTTAGCGAGCAGCTCCTTTCGAATCTGGGCATCGGTGACATTATCTATCAAAGCCAGGCAATGGAAGACCTGCTCCAGCAGGTAGATATTGTGGCGAAGAGCGACAGTACGGTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCA  >  NC_000913/2612826‑2613118

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