Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,016 A→G 9.3% D424G (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,0160AG 9.3% 202.3 / 9.9 75D424G (GAC→GGC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (5/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGT  >  NC_002947/1904868‑1905163
                                                                                                                                                    |                                                                                                                                                   
cTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGa                                                                                                                                                     >  7:110869/1‑149 (MQ=255)
 tCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGa                                                                                                                                                     >  4:32527/1‑148 (MQ=255)
   ggCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAg                                                                                                                                                  >  7:229854/1‑149 (MQ=255)
    gcgcGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGTCGGGCAgc                                                                                                                                                 >  3:241435/1‑149 (MQ=255)
         gTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGTCGGGCAgcgcatc                                                                                                                                            <  4:241435/149‑4 (MQ=255)
               gTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATc                                                                                                                                        <  4:51256/147‑1 (MQ=255)
                   gCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTcc                                                                                                                                  <  2:47842/149‑1 (MQ=255)
                    cGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGTCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCt                                                                                                                                 <  8:120878/149‑1 (MQ=255)
                     gACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTg                                                                                                                                >  6:107732/1‑149 (MQ=255)
                     gACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTg                                                                                                                                >  4:33022/1‑149 (MQ=255)
                      acaacaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGc                                                                                                                               <  5:143434/149‑1 (MQ=255)
                        aacaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGCCCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCca                                                                                                                             >  4:232497/1‑149 (MQ=255)
                         acaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCcac                                                                                                                            <  1:224360/149‑1 (MQ=255)
                         acaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCcac                                                                                                                            <  5:107732/149‑1 (MQ=255)
                            gCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAAc                                                                                                                         <  5:93320/149‑1 (MQ=255)
                             cTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACg                                                                                                                        >  3:15558/1‑149 (MQ=255)
                                ttACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCa                                                                                                                     <  3:32527/149‑1 (MQ=255)
                                  aCGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAg                                                                                                                    >  7:241302/1‑148 (MQ=255)
                                      cacCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTTGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAgg                                                                                                               <  5:264250/149‑1 (MQ=255)
                                      cacCCTGTTGCAGTACGCACAGTACGGCATCGCCCTGCCGGCCGACCGTGCAGTTCTTATTTTGCCTCAAGTCGGTTGTGCTACGGAGCTTGGCGCAGATGCGCTGACCTACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCaatg                                                                                                               <  8:108687/149‑3 (MQ=255)
                                        cccTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGt                                                                                                                       <  5:26609/139‑1 (MQ=255)
                                        cccTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGt                                                                                                                       >  6:26609/1‑139 (MQ=255)
                                           tGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCg                                                                                                          >  5:347277/1‑149 (MQ=255)
                                               gCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGgtgt                                                                                                      >  3:75537/1‑149 (MQ=255)
                                               gCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGgtgt                                                                                                      >  7:33508/1‑149 (MQ=255)
                                                      gCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGcgccgc                                                                                               <  4:66180/149‑1 (MQ=255)
                                                        ccAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGa                                                                                             >  8:65502/1‑149 (MQ=255)
                                                        ccAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGa                                                                                             >  4:139382/1‑149 (MQ=255)
                                                         cAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCt                                                                                                                                       <  7:30833/106‑1 (MQ=255)
                                                         cAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCt                                                                                                                                       >  8:30833/1‑106 (MQ=255)
                                                         cAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGAt                                                                                            >  7:264272/1‑149 (MQ=255)
                                                          aGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATc                                                                                           <  7:2203/149‑1 (MQ=255)
                                                          aGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATc                                                                                           >  3:246441/1‑149 (MQ=255)
                                                          aGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATc                                                                                           <  4:246441/149‑1 (MQ=255)
                                                            aaCGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGa                                                                                         <  8:33508/149‑1 (MQ=255)
                                                            aaCGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGGCAGCAAGGGGGCGGTGTGCGCCGCGATCGa                                                                                         <  2:333048/149‑1 (MQ=255)
                                                                 cATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTgc                                                                                                    <  5:119531/133‑1 (MQ=255)
                                                                 cATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTgc                                                                                                    >  6:119531/1‑133 (MQ=255)
                                                                      cccTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGc                                                                               >  3:357779/1‑149 (MQ=255)
                                                                           ccggccgACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGa                                                                          <  8:241302/149‑1 (MQ=255)
                                                                                     cGGCAGCGCGGATCTGGCCGCAAGGCGGGTGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGCCAGCAACGGGGCGGTGTGCGCCGCGATCGACAGCGTCTGCAGCGACCTCGCTGCg                                                                <  4:285471/149‑1 (MQ=255)
                                                                                           cgcTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAg                                                          <  4:235613/149‑1 (MQ=255)
                                                                                            gcTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAgg                                                         >  5:57166/1‑149 (MQ=255)
                                                                                             cTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAggg                                                        <  8:264272/149‑1 (MQ=255)
                                                                                             cTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAggg                                                        <  7:65502/149‑1 (MQ=255)
                                                                                              tGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGc                                                       >  6:194421/1‑149 (MQ=255)
                                                                                                  cTTGCCGCAAGGGGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCGCGACAGCGCGACCATCTGGTGCTGCAACAACGTCAGCAAGGGGGCGGTGTGGGCCGCGATCGACAGCGGGTGCAGCGAGCTCGCTGCGGTCAAGGTCTGc                                                    <  8:110869/148‑1 (MQ=255)
                                                                                                   ttGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTAAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCAc                                                  >  6:308370/1‑149 (MQ=255)
                                                                                                     gCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCAcc                                                 <  3:232497/148‑1 (MQ=255)
                                                                                                         cAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGc                                                                                                                               >  8:144054/1‑66 (MQ=255)
                                                                                                         cAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGc                                                                                                                               <  7:144054/66‑1 (MQ=255)
                                                                                                            ggCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCcgc                                         >  1:363819/1‑149 (MQ=255)
                                                                                                                      cgcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTg                                <  6:188217/148‑1 (MQ=255)
                                                                                                                       gcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCg                              >  2:338317/1‑149 (MQ=255)
                                                                                                                       gcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCg                              >  5:223363/1‑149 (MQ=255)
                                                                                                                           gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggtgg                          <  4:75537/149‑1 (MQ=255)
                                                                                                                           gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACAAAGGCCGCGACCGGCTGCggtgg                          <  6:57166/149‑1 (MQ=255)
                                                                                                                             ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                        <  4:357779/149‑1 (MQ=255)
                                                                                                                             ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                        >  7:227051/1‑149 (MQ=255)
                                                                                                                             ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTACAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                        >  5:242027/1‑149 (MQ=255)
                                                                                                                              cTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTg                       >  6:103176/1‑149 (MQ=255)
                                                                                                                                gggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCgg                                   >  5:118804/1‑135 (MQ=255)
                                                                                                                                gggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCgg                                   <  6:118804/135‑1 (MQ=255)
                                                                                                                                 ggCGCCGATGCGCTGCCCGACTGTGCGACCATCTGCTCCTGCCAGATCGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCTACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGTCTGCGGTGGTTgcgc                    <  8:229854/149‑1 (MQ=255)
                                                                                                                                    gcCGATGCGCGGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGc                 <  8:227051/149‑1 (MQ=255)
                                                                                                                                    atggtcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGc                 >  7:12095/7‑149 (MQ=255)
                                                                                                                                       gtcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCt              >  7:198608/4‑149 (MQ=255)
                                                                                                                                        tcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTgcgc                    >  4:38523/3‑142 (MQ=255)
                                                                                                                                        tcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTgcgc                    <  3:38523/140‑1 (MQ=255)
                                                                                                                                         agcgcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGc            <  3:60444/148‑1 (MQ=255)
                                                                                                                                          gcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCt           >  4:265495/1‑149 (MQ=255)
                                                                                                                                            gcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc          >  6:188338/1‑148 (MQ=255)
                                                                                                                                            gcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc          >  6:374652/1‑148 (MQ=255)
                                                                                                                                                cccGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagc      >  3:108445/1‑148 (MQ=255)
                                                                                                                                                 ccGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagcaa    <  6:347277/149‑1 (MQ=255)
                                                                                                                                                  cGAAAGCGCGACCATCTGCTCCTGCGACTACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCTGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagcaag   <  6:223363/149‑1 (MQ=255)
                                                                                                                                                   gACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGt  >  2:15087/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGT  >  NC_002947/1904868‑1905163

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: