Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,093 C→G 6.9% L450V (CTC→GTC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,0930CG 6.9% 215.8 / 6.6 72L450V (CTC→GTC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base C (29/38);  new base G (2/3);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

GCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCA  >  NC_002947/1904946‑1905238
                                                                                                                                                   |                                                                                                                                                 
gccgACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCt                                                                                                                                                  >  6:120704/1‑149 (MQ=255)
    aCCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCt                                                                                                                                              <  4:17546/149‑1 (MQ=255)
    aCCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCt                                                                                                                                              >  2:192743/1‑149 (MQ=255)
        ggCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGAGGTCGCTGcag                                                                                                                                          <  4:172410/149‑3 (MQ=255)
        ggCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCgg                                                                                                                                          >  5:125071/1‑149 (MQ=255)
        ggCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCCGCAGCGAGGTCGCTGcag                                                                                                                                          <  5:114513/149‑3 (MQ=255)
        ggCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCCGCAGCGAGGTCGCTGcag                                                                                                                                          >  6:114513/1‑147 (MQ=255)
          cAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCCGCAGCGAGGTCGCTGcagcc                                                                                                                                        <  3:159515/149‑5 (MQ=255)
             cgcTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAg                                                                                                                                     <  1:34377/149‑1 (MQ=255)
                   tCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGc                                                                                                                               >  7:97432/1‑149 (MQ=255)
                   gcTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGc                                                                                                                               <  8:176817/148‑1 (MQ=255)
                       gCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTg                                                                                                                                <  1:159494/144‑1 (MQ=255)
                        ccGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGc                                                                                                                                            >  5:50413/1‑131 (MQ=255)
                        ccGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGc                                                                                                                                            <  6:50413/131‑1 (MQ=255)
                        ccGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCaa                                                                                                                          <  6:164523/149‑1 (MQ=255)
                         cGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAg                                                                                                                         >  1:247740/1‑149 (MQ=255)
                         cGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAg                                                                                                                         >  1:163343/1‑149 (MQ=255)
                           cAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGc                                                                                                                       <  2:50441/149‑1 (MQ=255)
                            aaGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGcc                                                                                                                      <  7:137260/149‑1 (MQ=255)
                              ggCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCcgc                                                                                                                    <  6:125071/149‑1 (MQ=255)
                                         gcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCg                                                                                                         >  7:168577/1‑149 (MQ=255)
                                          ccGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCGCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCgg                                                                                                        <  6:214141/149‑1 (MQ=255)
                                           cGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggt                                                                                                       <  1:243647/149‑1 (MQ=255)
                                              cccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggtggt                                                                                                    <  6:179952/149‑1 (MQ=255)
                                               ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCAc                                                                                                                             <  3:208568/123‑1 (MQ=255)
                                               ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCAc                                                                                                                             >  4:208568/1‑123 (MQ=255)
                                                cTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGc                                                                                                                               <  5:177783/120‑1 (MQ=255)
                                                cTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGc                                                                                                                               >  6:177783/1‑120 (MQ=255)
                                                   ggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTgcgc                                                                                               <  2:96789/149‑1 (MQ=255)
                                                   ggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTgcgc                                                                                               >  5:23129/1‑149 (MQ=255)
                                                     cgcCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggtgg                                                                                                     >  4:131916/1‑141 (MQ=255)
                                                     cgcCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggtgg                                                                                                     <  3:131916/141‑1 (MQ=255)
                                                      gcCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGc                                                                                            >  1:283314/1‑149 (MQ=255)
                                                       ccGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGcc                                                                                           <  7:157453/149‑1 (MQ=255)
                                                         gtcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCt                                                                                         >  7:34379/4‑149 (MQ=255)
                                                          aTGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTg                                                                                        <  8:97432/149‑1 (MQ=255)
                                                             cgcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGGCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc                                                                                     <  6:23129/149‑1 (MQ=255)
                                                              gcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc                                                                                     >  1:282524/1‑148 (MQ=255)
                                                               cTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaa                                                                                   <  2:282524/149‑1 (MQ=255)
                                                                   ccGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagc                                                                                 >  3:236677/1‑147 (MQ=255)
                                                                    cGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagcaag                                                                              <  7:133555/149‑1 (MQ=255)
                                                                    cGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagcaag                                                                              <  8:143861/149‑1 (MQ=255)
                                                                         gagcgACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTc                                                                         >  7:200549/3‑149 (MQ=255)
                                                                          cgcgACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCg                                                                        <  8:34379/149‑1 (MQ=255)
                                                                             gACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAg                                                                      >  7:11265/1‑148 (MQ=255)
                                                                               ccATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTc                                                                         >  5:13077/1‑143 (MQ=255)
                                                                               ccATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTc                                                                         <  6:13077/143‑1 (MQ=255)
                                                                               ccATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTc                                                                         >  5:8704/1‑143 (MQ=255)
                                                                               ccATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCAc                                                                   <  2:247740/149‑1 (MQ=255)
                                                                               ccATCTGCTCCAGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTc                                                                         <  6:8704/143‑1 (MQ=255)
                                                                                      cTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGAcc                                                                                                                >  7:247926/1‑97 (MQ=255)
                                                                                      cTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGAcc                                                                                                                <  8:247926/97‑1 (MQ=255)
                                                                                          tGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGc                                                        >  3:285192/1‑149 (MQ=255)
                                                                                           gCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGcc                                                       >  6:180111/1‑149 (MQ=255)
                                                                                            ccACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGTTCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGcc                                                       <  2:283314/148‑1 (MQ=255)
                                                                                                      aGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTgg                                             <  5:180111/148‑1 (MQ=255)
                                                                                                       gCAAGGGGGCGGTGTGCGCAGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACAAAGGCCGCGAACGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCTGCGTGGCg                                           >  2:193728/1‑149 (MQ=255)
                                                                                                           gggggCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTgg                                       <  4:236677/149‑1 (MQ=255)
                                                                                                                  gtgtGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAagag                                >  2:165577/1‑149 (MQ=255)
                                                                                                                   tgtgCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGc                               <  5:232422/149‑1 (MQ=255)
                                                                                                                        gccgcGAGCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCGGCGACCGGCTGCGGTGGTTGCGCTGCGCTGCTCAAGCAAGTCTTCGAGCTCGATCTGACGGCCGGCGGCGTGGCGGTGGACAAGAGCCtgtg                          <  2:180386/149‑1 (MQ=255)
                                                                                                                         ccgcGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGc                         <  1:192743/149‑1 (MQ=255)
                                                                                                                               tCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCAc                   >  5:129168/1‑149 (MQ=255)
                                                                                                                                  aCAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACttt                <  4:285192/149‑1 (MQ=255)
                                                                                                                                   cAGAGGCTGCAGCGACCTCCCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGAGCGGCCCTGCTCAAGCAAGTCTACGAGCACGAACTCACCGCCCGCGGCGTGGCGGAGGACAAGAGCCTGTGCGAGAACATTg               >  4:127878/1‑149 (MQ=255)
                                                                                                                                    aGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGc              >  2:262217/1‑149 (MQ=255)
                                                                                                                                    aGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGc              <  8:11265/149‑1 (MQ=255)
                                                                                                                                     gCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAAGTGACGGCCCGCGGCGTGGCGGTTGACAAGAGCCTGTGCGAGCACTTTGcc             <  3:127878/149‑1 (MQ=255)
                                                                                                                                       ggCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTa           <  2:163343/149‑1 (MQ=255)
                                                                                                                                        accgCAGCGAGGTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTa           >  5:44499/4‑148 (MQ=255)
                                                                                                                                                gACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCa  <  1:262217/149‑1 (MQ=255)
                                                                                                                                                gACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACATTGCCTACACCCgg    >  7:285224/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCA  >  NC_002947/1904946‑1905238

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: