Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,421,749 C→T R14R (CGC→CGT yfcH → putative NAD‑dependent nucleotide‑sugar epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,7490CT98.6% 263.4 / ‑4.1 73R14R (CGC→CGTyfcHputative NAD‑dependent nucleotide‑sugar epimerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/40);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGT  >  NC_000913/2421601‑2421895
                                                                                                                                                    |                                                                                                                                                  
gTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCa                                                                                                                                                  <  1:848868/151‑1 (MQ=255)
       aTGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGAtt                                                                                                                                           >  1:1015962/1‑151 (MQ=255)
       aTGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGAtt                                                                                                                                           >  1:675323/1‑151 (MQ=255)
         gCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTcc                                                                                                                                         <  2:1238517/151‑1 (MQ=255)
         gCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTcc                                                                                                                                         >  2:737836/1‑151 (MQ=255)
            gAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCcg                                                                                                                                        <  2:750867/149‑1 (MQ=255)
              aGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGtt                                                                                                                                    >  1:731595/1‑151 (MQ=255)
               gTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGttt                                                                                                                                   >  1:1116521/1‑151 (MQ=255)
                 ggAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTg                                                                                                                                  <  1:406890/150‑1 (MQ=255)
                 ggAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTg                                                                                                                                  <  1:613205/150‑1 (MQ=255)
                   agaTCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTg                                                                                                                               <  1:317065/151‑1 (MQ=255)
                    gaTCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTgg                                                                                                                              >  2:444655/1‑151 (MQ=255)
                            aaCTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGc                                                                                                                      <  2:144859/151‑1 (MQ=255)
                              cTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGcc                                                                                                                     >  2:839080/1‑150 (MQ=255)
                                gCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATc                                                                                                                  >  2:914775/1‑151 (MQ=255)
                                gCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATc                                                                                                                  <  1:624233/151‑1 (MQ=255)
                                     cgctgcgctACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAAt                                                                                                              <  2:731926/150‑1 (MQ=255)
                                        tgcgctACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACg                                                                                                          >  1:1186012/1‑151 (MQ=255)
                                          cgctACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACgg                                                                                                         >  1:879286/1‑150 (MQ=255)
                                          cgctACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGt                                                                                                        <  1:303389/151‑1 (MQ=255)
                                           gctACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACgg                                                                                                         >  2:512171/1‑149 (MQ=255)
                                               cGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGa                                                                                                   <  1:560503/151‑1 (MQ=255)
                                                gCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGa                                                                                                   <  2:664389/150‑1 (MQ=255)
                                                gCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGAc                                                                                                  <  2:211477/151‑1 (MQ=255)
                                                 ccGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGa                                                                                                   >  2:533361/1‑149 (MQ=255)
                                                   gATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGAcgcg                                                                                               <  1:444026/151‑1 (MQ=255)
                                                    aTTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGt                                                                                              >  1:1022151/1‑151 (MQ=255)
                                                    aTTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGt                                                                                              >  2:538426/1‑151 (MQ=255)
                                                         gTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAAccc                                                                                         <  2:1158410/151‑1 (MQ=255)
                                                          tATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCg                                                                                        >  2:1193618/1‑151 (MQ=255)
                                                           aTCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGc                                                                                       <  2:353685/151‑1 (MQ=255)
                                                            tCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGc                                                                                       >  1:1154/1‑150 (MQ=255)
                                                                 gACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGaaa                                                                                  <  2:32550/150‑1 (MQ=255)
                                                                  aCCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGc                                                                                >  2:684918/1‑151 (MQ=255)
                                                                  aCCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGc                                                                                >  2:761368/1‑151 (MQ=255)
                                                                     ttAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGcc                                                                               >  2:245321/1‑149 (MQ=255)
                                                                       aaGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGt                                                                            <  2:925519/150‑1 (MQ=255)
                                                                        aGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTc                                                                          >  2:251402/1‑151 (MQ=255)
                                                                               aGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTtctc                                                                   >  2:1246364/1‑151 (MQ=255)
                                                                               aGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTtctc                                                                   <  1:295132/151‑1 (MQ=255)
                                                                                gCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCg                                                                  <  1:1032567/151‑1 (MQ=255)
                                                                                gCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCg                                                                  >  1:770088/1‑151 (MQ=255)
                                                                                      aTCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTc                                                            >  2:877381/1‑151 (MQ=255)
                                                                                        cGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCg                                                           <  1:217811/150‑1 (MQ=255)
                                                                                              atCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTgtc                                                     >  1:1150369/1‑148 (MQ=255)
                                                                                              atCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGac                                                     >  2:389314/1‑150 (MQ=255)
                                                                                                 gggAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTa                                                 <  1:638117/151‑1 (MQ=255)
                                                                                                    aGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATgg                                              <  2:1223306/151‑1 (MQ=255)
                                                                                                     ggCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGc                                             >  2:260478/1‑151 (MQ=255)
                                                                                                     ggCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGc                                             >  2:443365/1‑151 (MQ=255)
                                                                                                         gCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAgg                                         >  2:290049/1‑151 (MQ=255)
                                                                                                         gCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAg                                          <  2:1262275/150‑1 (MQ=255)
                                                                                                           aTGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGc                                       <  2:849563/151‑1 (MQ=255)
                                                                                                           aTGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGc                                       >  2:327576/1‑151 (MQ=255)
                                                                                                             gAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCtt                                     <  2:204691/151‑1 (MQ=255)
                                                                                                                  tAGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCg                                 >  1:412696/1‑150 (MQ=255)
                                                                                                                   aGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGa                                <  1:1163194/150‑1 (MQ=255)
                                                                                                                   aGTGATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGAt                               <  1:830595/151‑1 (MQ=255)
                                                                                                                      gATCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCa                             <  1:975280/150‑1 (MQ=255)
                                                                                                                       aTCACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCa                             >  2:820873/1‑149 (MQ=255)
                                                                                                                         cACCGGAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAg                          <  2:162014/150‑1 (MQ=255)
                                                                                                                              gAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCt                    <  1:468040/151‑1 (MQ=255)
                                                                                                                              gAGGGACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCt                    <  2:1119654/151‑1 (MQ=255)
                                                                                                                                  gACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAAc                <  2:564478/151‑1 (MQ=255)
                                                                                                                                  gACAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAAc                <  1:914723/151‑1 (MQ=255)
                                                                                                                                   aCAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACg               <  2:732443/151‑1 (MQ=255)
                                                                                                                                    cAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACgg              <  2:319921/151‑1 (MQ=255)
                                                                                                                                    cAGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACgg              <  2:1054020/151‑1 (MQ=255)
                                                                                                                                     aGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACgg              >  2:1026543/1‑150 (MQ=255)
                                                                                                                                     aGGATTAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGc             <  2:1146938/151‑1 (MQ=255)
                                                                                                                                          tAATTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTg         >  1:219366/1‑150 (MQ=255)
                                                                                                                                           aaTTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGAt       <  1:361926/151‑1 (MQ=255)
                                                                                                                                            aTTGGTCGTCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATg      <  2:737627/151‑1 (MQ=255)
                                                                                                                                                gtcgtcATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGt  >  2:208888/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTGACGTATGCTGAAGTGGAGATCGATAAACTGCACGCGCTGCGCTACGCCGATTCGGTATCCATGACCTTAAGCTGGCAGCGTTAATCGCCATATCGGGAGGCTGCATGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGT  >  NC_000913/2421601‑2421895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: