Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,012 C→T 9.6% P423S (CCC→TCC) ‡ nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,0120CT 9.6% 204.4 / 7.8 73P423S (CCC→TCC) ‡nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (5/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CGTGCTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGC  >  NC_002947/1904864‑1905159
                                                                                                                                                    |                                                                                                                                                   
cGTGCTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGc                                                                                                                                                     <  7:69386/149‑1 (MQ=255)
    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                                                                                                                            >  4:33022/1‑149 (MQ=255)
                         gACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTg                                                                                                                            >  6:107732/1‑149 (MQ=255)
                          acaacaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGc                                                                                                                           <  5:143434/149‑1 (MQ=255)
                            aacaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGCCCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCca                                                                                                                         >  4:232497/1‑149 (MQ=255)
                             acaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCcac                                                                                                                        <  5:107732/149‑1 (MQ=255)
                             acaGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCcac                                                                                                                        <  1:224360/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                                                                                                                   >  6:26609/1‑139 (MQ=255)
                                            cccTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGt                                                                                                                   <  5:26609/139‑1 (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                                                                                                                                   >  8:30833/1‑106 (MQ=255)
                                                             cAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCt                                                                                                                                   <  7:30833/106‑1 (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                                                    <  7:65502/149‑1 (MQ=255)
                                                                                                 cTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAggg                                                    <  8:264272/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                                                                                                                           <  7:144054/66‑1 (MQ=255)
                                                                                                             cAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGc                                                                                                                           >  8:144054/1‑66 (MQ=255)
                                                                                                                ggCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCcgc                                     >  1:363819/1‑149 (MQ=255)
                                                                                                                          cgcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTg                            <  6:188217/148‑1 (MQ=255)
                                                                                                                           gcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCg                          >  5:223363/1‑149 (MQ=255)
                                                                                                                           gcCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCg                          >  2:338317/1‑149 (MQ=255)
                                                                                                                               gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCggtgg                      <  4:75537/149‑1 (MQ=255)
                                                                                                                               gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACAAAGGCCGCGACCGGCTGCggtgg                      <  6:57166/149‑1 (MQ=255)
                                                                                                                                 ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                    >  7:227051/1‑149 (MQ=255)
                                                                                                                                 ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                    <  4:357779/149‑1 (MQ=255)
                                                                                                                                 ccTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTACAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGtt                    >  5:242027/1‑149 (MQ=255)
                                                                                                                                  cTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTg                   >  6:103176/1‑149 (MQ=255)
                                                                                                                                    gggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCgg                               <  6:118804/135‑1 (MQ=255)
                                                                                                                                    gggCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCgg                               >  5:118804/1‑135 (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                <  3:38523/140‑1 (MQ=255)
                                                                                                                                            tcgcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTgcgc                >  4:38523/3‑142 (MQ=255)
                                                                                                                                             agcgcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGc        <  3:60444/148‑1 (MQ=255)
                                                                                                                                              gcgcTGTCGGGCAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCt       >  4:265495/1‑149 (MQ=255)
                                                                                                                                                gcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc      >  6:374652/1‑148 (MQ=255)
                                                                                                                                                gcTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTc      >  6:188338/1‑148 (MQ=255)
                                                                                                                                                    cccGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaagc  >  3:108445/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTGCTCGGCGCGGTGCTGGTCGGCGACAACAGCTATTACGACACCCTGTTGCAGTACGCCCAGAACGGCATCGCCCTGCCGGCCGACCCGGCAGCGCTGATCTTGCCGCAAGGCGGTGGTGCGCCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGCACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGC  >  NC_002947/1904864‑1905159

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: ‑