Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,409 A→G 5.9% D555G (GAC→GGC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,4090AG 5.9% 289.9 / 6.3 101D555G (GAC→GGC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base A (51/44);  new base G (3/3);  total (54/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 104 total aligned reads displayed.

TGGAGGGTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGT  >  NC_002947/1905262‑1905557
                                                                                                                                                    |                                                                                                                                                     
tGGAGGGTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGac                                                                                                                                                      <  5:442034/149‑1 (MQ=255)
tGGAGGGTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGa                                                                                                                                                       >  5:113280/1‑148 (MQ=255)
   acggTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGTGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACAcc                                                                                                                                                   <  6:238277/147‑1 (MQ=255)
      gTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGTCCTGCAGCGGCACCAAc                                                                                                                                                <  5:418698/149‑1 (MQ=255)
         tCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGTCCTGCAGCGGCACCAACgc                                                                                                                                              <  5:156775/148‑2 (MQ=255)
          ccGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACCacc                                                                                                                                            >  7:354397/1‑148 (MQ=255)
          ccGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGacac                                                                                                                                           >  2:238501/1‑149 (MQ=255)
                tttCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCa                                                                                                                                      >  8:68735/1‑149 (MQ=255)
                tttCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCa                                                                                                                                      >  2:155362/1‑149 (MQ=255)
                    aGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGc                                                                                                                                  <  4:358601/149‑1 (MQ=255)
                      cGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCAt                                                                                                                                     <  7:187235/144‑1 (MQ=255)
                      cGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCAt                                                                                                                                     >  8:187235/1‑144 (MQ=255)
                      cGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGcc                                                                                                                                 >  3:337548/1‑148 (MQ=255)
                         aCTGCTCGAGCGCCATTGCAAGGGCCATCTGGGGTGTGACATCTTCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCAGATCATCGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACa                                                                                                                             >  2:120415/1‑149 (MQ=255)
                         aCTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACa                                                                                                                             >  4:198678/1‑149 (MQ=255)
                         aCTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACa                                                                                                                             >  3:403852/1‑149 (MQ=255)
                         aCTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACa                                                                                                                             >  4:240500/1‑149 (MQ=255)
                           tGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTACAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATg                                                                                                                           >  5:227670/1‑149 (MQ=255)
                             cTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATg                                                                                                                           >  1:411440/1‑147 (MQ=255)
                             cTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCa                                                                                                                         <  2:411440/149‑1 (MQ=255)
                              tCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCCACATTCTCGCCTCGTGCTGGAGCCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATg                                                                                                                           >  5:300245/1‑146 (MQ=255)
                                gAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGa                                                                                                                       >  1:60399/1‑148 (MQ=255)
                                         ggCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCAc                                                                                                             >  4:16843/1‑149 (MQ=255)
                                          gCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCAc                                                                                                             <  6:113280/148‑1 (MQ=255)
                                          gCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCAc                                                                                                             <  3:200829/148‑1 (MQ=255)
                                          gCAAGGGCCATCTGGGGGGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCAcc                                                                                                            <  8:233123/149‑1 (MQ=255)
                                            aaGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTa                                                                                                          >  4:56768/1‑149 (MQ=255)
                                                        gggTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGAGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCggtggt                                                                                                 <  1:120415/147‑1 (MQ=255)
                                                        gggTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACg                                                                                                                                               >  8:161048/1‑100 (MQ=255)
                                                        gggTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACg                                                                                                                                               <  7:161048/100‑1 (MQ=255)
                                                        gggTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGcc                                                                                              >  1:420436/1‑149 (MQ=255)
                                                         ggTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcg                                                                                             >  2:282592/1‑149 (MQ=255)
                                                          gtgtGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcgc                                                                                            <  4:375680/149‑1 (MQ=255)
                                                           tgtgACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcgcg                                                                                           >  1:400379/1‑149 (MQ=255)
                                                             tgACATCGGCACGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGGCCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCa                                                                                         <  8:139095/149‑1 (MQ=255)
                                                              gACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcg                                                                                             <  8:277626/144‑1 (MQ=255)
                                                              gACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcg                                                                                             >  7:277626/1‑144 (MQ=255)
                                                               aCATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCAtt                                                                                       >  5:144862/1‑149 (MQ=255)
                                                                cATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcgc                                                                                            >  5:207973/1‑143 (MQ=255)
                                                                cATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCcgc                                                                                            <  6:207973/143‑1 (MQ=255)
                                                                 aTCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTcc                                                                                     >  6:349337/1‑149 (MQ=255)
                                                                   cTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCg                                                                                   >  8:279399/1‑149 (MQ=255)
                                                                    tGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCg                                                                                   <  1:155362/148‑1 (MQ=255)
                                                                    tGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGAACCGGCGTTGGTGCCGCTGCAGGACCCCAACGACACCGTCATGGCCAACATGCAGAAAAACGGCACCAATTCGGTGGTGCCGCGCAATCCCgg                                                                                  >  2:212320/1‑149 (MQ=255)
                                                                     gCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggt                                                                                 >  3:19836/1‑149 (MQ=255)
                                                                      cAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtg                                                                                <  7:89722/149‑1 (MQ=255)
                                                                       actcctgcggTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtgg                                                                               <  6:300245/143‑1 (MQ=255)
                                                                         gCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGgtggtg                                                                             <  7:279399/149‑1 (MQ=255)
                                                                           cggcggTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGa                                                                            >  7:350506/1‑148 (MQ=255)
                                                                           cggcggTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGa                                                                            >  8:128445/1‑148 (MQ=255)
                                                                                gTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAc                                                                      <  4:403852/149‑1 (MQ=255)
                                                                                  gggCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGTAATAAccc                                                                    <  6:93134/149‑1 (MQ=255)
                                                                                   ggCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAcccc                                                                   >  2:254251/1‑149 (MQ=255)
                                                                                     cAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCg                                                                                   >  2:244268/1‑131 (MQ=255)
                                                                                     cAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCg                                                                                   <  1:244268/131‑1 (MQ=255)
                                                                                       aCATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGAc                                                               <  1:238501/149‑1 (MQ=255)
                                                                                       aCATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCACTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGa                                                                >  1:405456/1‑148 (MQ=255)
                                                                                                ccTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatc                                                      <  6:25345/149‑1 (MQ=255)
                                                                                                ccTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatc                                                      >  4:18376/1‑149 (MQ=255)
                                                                                                     tGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatcgtcat                                                >  2:399750/1‑149 (MQ=255)
                                                                                                        tGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCgg                                             >  6:216618/1‑149 (MQ=255)
                                                                                                        tGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCgg                                             >  1:333983/1‑149 (MQ=255)
                                                                                                            aCCAGCCGATCATGGACCCGTCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAg                                         <  4:337548/149‑1 (MQ=255)
                                                                                                            aCCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAg                                         <  1:112759/149‑1 (MQ=255)
                                                                                                            aCCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAg                                         >  2:112759/1‑149 (MQ=255)
                                                                                                            aCCAGCCGATAATGGACCCGGCGTTGGTGCCTCTTCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAg                                         <  7:68735/149‑1 (MQ=255)
                                                                                                             ccAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAgg                                        >  3:276500/1‑149 (MQ=255)
                                                                                                              cAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCTCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGt                                       <  6:144862/149‑1 (MQ=255)
                                                                                                              cAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGt                                       >  7:13990/1‑149 (MQ=255)
                                                                                                               aGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTg                                      >  1:435250/1‑149 (MQ=255)
                                                                                                                gCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTgg                                     <  3:276561/149‑1 (MQ=255)
                                                                                                                   gATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGcc                                   >  7:71034/1‑148 (MQ=255)
                                                                                                                     tCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAg                                <  5:332875/149‑1 (MQ=255)
                                                                                                                       aTGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGa                               <  8:350506/148‑1 (MQ=255)
                                                                                                                        tGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGaaa                             <  5:216618/149‑1 (MQ=255)
                                                                                                                          gACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAAt                            <  4:19836/148‑1 (MQ=255)
                                                                                                                           aCCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGCCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATAc                          <  6:441797/149‑1 (MQ=255)
                                                                                                                             ccGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGa                        >  1:302154/1‑149 (MQ=255)
                                                                                                                              cGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGAc                       <  3:56768/149‑1 (MQ=255)
                                                                                                                               ggCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACTGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGAcc                      <  6:227670/149‑1 (MQ=255)
                                                                                                                               ggCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGCTCATCGGCCAGGTGGCCAAGAAATACGAc                       <  1:254251/149‑1 (MQ=255)
                                                                                                                                tcGTTGGTGTCCTGCAGCGGCACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAACCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGAc                       >  8:209202/2‑147 (MQ=255)
                                                                                                                                 cGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTg                    <  2:435250/149‑1 (MQ=255)
                                                                                                                                  gTTGGTGTCCTGCAGCGGCACCAACGACACCTTCATGGCCAACATGCAGAATAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCCGGTGGCCAAGAAATa                           <  7:78713/141‑1 (MQ=255)
                                                                                                                                  gTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTg                    >  2:212890/1‑148 (MQ=255)
                                                                                                                                   ttGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATaa                                                                       <  4:10280/97‑1 (MQ=255)
                                                                                                                                   ttGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATaa                                                                       >  3:10280/1‑97 (MQ=255)
                                                                                                                                   ttGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTa                  <  1:282592/149‑1 (MQ=255)
                                                                                                                                      gTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTaca                <  3:16843/148‑1 (MQ=255)
                                                                                                                                       tGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACAcc              <  2:302154/149‑1 (MQ=255)
                                                                                                                                       tGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACAcc              <  3:240500/149‑1 (MQ=255)
                                                                                                                                        gCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACAcc              >  6:219148/1‑148 (MQ=255)
                                                                                                                                        gCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACAcc              >  6:377701/1‑148 (MQ=255)
                                                                                                                                            cTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACAcc              >  3:198747/1‑144 (MQ=255)
                                                                                                                                            cTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGan         <  2:405456/149‑2 (MQ=255)
                                                                                                                                                aGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCAcc     <  2:60399/149‑1 (MQ=255)
                                                                                                                                                   gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGGGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCgg   >  3:411656/1‑149 (MQ=255)
                                                                                                                                                   gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAAt                            <  4:367542/124‑1 (MQ=255)
                                                                                                                                                   gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAAt                            >  3:367542/1‑124 (MQ=255)
                                                                                                                                                    acaacaACGACACCCTCAATGCCAACATGCACAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGt  <  3:243254/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGAGGGTATCCGCACTTTCAGCGAACTGCTCGAGCGCCATGGCAAGGGCCATCTGGGGTGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGT  >  NC_002947/1905262‑1905557

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