Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,906,456 G→T 6.5% V53V (GTG→GTT)  nirD → nitrite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,906,4560GT 6.5% 264.6 / 3.6 84V53V (GTG→GTT) nirDnitrite reductase
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base T (2/4);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGAACACGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGG  >  NC_002947/1906318‑1906597
                                                                                                                                          |                                                                                                                                             
tgtgAACACGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAAcc                                                                                                                                        >  5:313998/1‑146 (MQ=255)
tgtgAACACGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAAcc                                                                                                                                        <  6:313998/146‑1 (MQ=255)
tgtgAACACGCAAGAAAACTGGCAGACGGTGTGCCAGCCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGTTGGACAACcgcg                                                                                                                                     >  4:372252/1‑149 (MQ=255)
   gaccacGCAAGAAACCGGGCAGACGGTGTGCCAGGCACAAGACCTCGTGGCCGAATCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTTGTCCACcgcgtac                                                                                                                                  <  2:304960/146‑4 (MQ=255)
    aaCACGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTTGTCCACcgcgtaca                                                                                                                                 <  7:54647/149‑5 (MQ=255)
     acacGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCa                                                                                                                                >  4:280026/1‑149 (MQ=255)
       acGCAATAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACg                                                                                                                              >  2:273652/1‑149 (MQ=255)
       acGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACg                                                                                                                              >  7:181946/1‑149 (MQ=255)
         gCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCt                                                                                                                            >  6:101358/1‑149 (MQ=255)
           aaGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTcc                                                                                                                          >  3:343689/1‑149 (MQ=255)
            aGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAAcc                                                                                                                                        <  1:90418/134‑1 (MQ=255)
            aGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAAcc                                                                                                                                        >  2:90418/1‑134 (MQ=255)
              aaaaCTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCt                                                                                                                            >  4:387796/1‑144 (MQ=255)
              aaaaCTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCt                                                                                                                            <  3:387796/144‑1 (MQ=255)
                 ccTGGCACACCGTCTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGc                                                                                                                     >  1:222976/2‑148 (MQ=255)
                 aCTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGcc                                                                                                                    <  6:146048/149‑1 (MQ=255)
                 aCTGGCAGACGGTGGGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTTGACAACCGCGACCCACGCTCCGGTGcc                                                                                                                    <  5:237866/149‑1 (MQ=255)
                  cTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCa                                                                                                                   <  7:287239/149‑1 (MQ=255)
                   tGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCaa                                                                                                                  >  1:136054/1‑149 (MQ=255)
                   tGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGACCAGCCGCTGTACGAGGTTGACAACCGCGACCCACGCTCCGGTGCCaa                                                                                                                  >  8:130986/1‑149 (MQ=255)
                   tGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTATACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGCCACCCGCGACCCACGCTCCGGTGGCaa                                                                                                                  >  4:47018/1‑149 (MQ=255)
                         aCGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATatcat                                                                                                            <  8:296417/149‑1 (MQ=255)
                          cGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATatcatc                                                                                                           <  1:273652/149‑1 (MQ=255)
                              gtgCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTc                                                                                                       <  4:343689/149‑1 (MQ=255)
                                 ccAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTg                                                                                                    <  3:280026/149‑1 (MQ=255)
                                 ccAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTg                                                                                                    >  3:312470/1‑149 (MQ=255)
                                   aGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTggg                                                                                                  >  5:40180/1‑149 (MQ=255)
                                   aGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTCACGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTggg                                                                                                  >  1:330894/1‑149 (MQ=255)
                                     gCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCt                                                                                                >  2:331740/1‑149 (MQ=255)
                                      cacaAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGTGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGGTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTTGGCtg                                                                                               <  4:323356/149‑1 (MQ=255)
                                       acaAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCtgg                                                                                              >  3:95236/1‑149 (MQ=255)
                                       acaAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCtgg                                                                                              <  4:312470/149‑1 (MQ=255)
                                        caAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCtggt                                                                                             >  5:389850/1‑149 (MQ=255)
                                        caAGACCTGGTGGCCGACTCCGGCGTCGTGGGTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTTGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCtggt                                                                                             <  5:319054/149‑1 (MQ=255)
                                            aCCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGc                                                                                         <  2:407208/149‑1 (MQ=255)
                                             gcTGGTGCCCGTCCTCGGCGTCGTGGTTTGGCTCGAGGTGGCCCAGTTGGCGGTTTTCTACCTGCCGGGGCTGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCa                                                                                        <  4:304996/148‑1 (MQ=255)
                                               tggtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTg                                                                                                    <  7:301532/135‑1 (MQ=255)
                                               tggtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTg                                                                                                    >  8:301532/1‑135 (MQ=255)
                                                 gtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    <  3:128950/149‑1 (MQ=255)
                                                 gtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    >  2:308005/1‑149 (MQ=255)
                                                 gtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    <  6:40180/149‑1 (MQ=255)
                                                 gtggCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    >  2:302879/1‑149 (MQ=255)
                                                     ccGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCaa                                                                                <  8:181946/149‑1 (MQ=255)
                                                      cGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAg                                                                               >  1:225696/1‑149 (MQ=255)
                                                          tCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACTACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGt                                                                             >  7:180318/1‑147 (MQ=255)
                                                            cGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    >  7:32103/1‑138 (MQ=255)
                                                            cGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    <  8:32103/138‑1 (MQ=255)
                                                            cGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGt                                                                         >  5:21760/1‑149 (MQ=255)
                                                             ggCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGtt                                                                        <  5:101358/149‑1 (MQ=255)
                                                                 tcgtGGGTTGGCTCGATGGGGCCCAGGGGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtg                                                                    <  2:264457/149‑1 (MQ=255)
                                                                  cgtGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtgg                                                                   <  1:308005/149‑1 (MQ=255)
                                                                  cgtGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtgg                                                                   <  1:302879/149‑1 (MQ=255)
                                                                      gTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtggtgg                                                                >  3:266888/1‑148 (MQ=255)
                                                                            cTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGcgcc                                                         >  6:375030/1‑149 (MQ=255)
                                                                            cTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGcgcc                                                         <  1:331740/149‑1 (MQ=255)
                                                                               gATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGct                                                      <  2:222976/149‑1 (MQ=255)
                                                                                aTGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    >  5:362177/1‑118 (MQ=255)
                                                                                aTGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAgtgt                                                                                    <  6:362177/118‑1 (MQ=255)
                                                                                   gggCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTAc                                                  <  2:205944/149‑1 (MQ=255)
                                                                                     gCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACaa                                                <  6:389850/149‑1 (MQ=255)
                                                                                      cccAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAag                                               <  4:266888/149‑1 (MQ=255)
                                                                                         aGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAagcag                                            <  5:375030/149‑1 (MQ=255)
                                                                                              gCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACtt                                       >  4:374988/1‑149 (MQ=255)
                                                                                                gttgttCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAagcagccatt                                       >  7:163713/1‑143 (MQ=255)
                                                                                                gttgttCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCa                                     <  2:136054/149‑1 (MQ=255)
                                                                                                    ttCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTCGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCct                                 >  6:165003/1‑149 (MQ=255)
                                                                                                      cTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCctgc                               >  3:189330/1‑149 (MQ=255)
                                                                                                       tACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGCCCCACGCTCCGGTGACAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGAGAGTTGGTGGGCGCAGCGCCGCGCATCAAGCACAACTTCAGCctgcg                              >  8:365183/1‑148 (MQ=255)
                                                                                                         ccTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGa                            <  3:374988/149‑1 (MQ=255)
                                                                                                          cTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCATGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAg                           >  6:113679/1‑149 (MQ=255)
                                                                                                          cTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGct                                                      >  6:382726/1‑122 (MQ=255)
                                                                                                          cTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGct                                                      <  5:382726/122‑1 (MQ=255)
                                                                                                           tGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtggtgg                                                                <  4:138872/111‑1 (MQ=255)
                                                                                                           tGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTtggtggtgg                                                                >  3:138872/1‑111 (MQ=255)
                                                                                                              cGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGGCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTggg                       <  4:324278/149‑1 (MQ=255)
                                                                                                                        ccAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTTAGTGGGGAGTGTCTgg             <  6:399066/149‑1 (MQ=255)
                                                                                                                         cAGCCGCTGTACGCGGTGGACAACCGCGACCGACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTgga            <  4:95236/149‑1 (MQ=255)
                                                                                                                         cAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCGGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTgga            <  3:372252/149‑1 (MQ=255)
                                                                                                                         cAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTgga            <  2:330894/149‑1 (MQ=255)
                                                                                                                         cAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTg              >  2:396112/1‑147 (MQ=255)
                                                                                                                             cGCTGTACGCGGTGGACTACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGAt        <  8:180318/149‑1 (MQ=255)
                                                                                                                             cGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGAt        <  4:189330/149‑1 (MQ=255)
                                                                                                                                  gtcgcgGTTGTCCACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGAGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCgg   >  4:154891/3‑149 (MQ=255)
                                                                                                                                    cgcgGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCggg  >  4:94239/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                             
TGTGAACACGCAAGAAAACTGGCAGACGGTGTGCCAGGCACAAGACCTGGTGGCCGACTCCGGCGTCGTGGTTTGGCTCGATGGGGCCCAGGTGGCGTTGTTCTACCTGCCGGGGCAGGCCCAGCCGCTGTACGCGGTGGACAACCGCGACCCACGCTCCGGTGCCAATATCATCGGTCGTGGGCTGGTGGGCAGTGTGCAAGGTGAGTTGGTGGTGGCAGCGCCGCTCTACAAGCAGCACTTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGG  >  NC_002947/1906318‑1906597

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