Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,212,276 T→C 6.7% *189W (TGA→TGG) ‡ PP_2816 ← transcriptional regulator NfxB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,212,2760TC 6.7% 218.7 / 4.7 74*189W (TGA→TGG) ‡PP_2816transcriptional regulator NfxB
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (3/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATCCCGCTTGCCGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAACGCCGCGGTGACAT  >  NC_002947/3212131‑3212424
                                                                                                                                                  |                                                                                                                                                     
aTCCCGCTTGCCGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCaa                                                                                                                                                    >  2:261769/1‑149 (MQ=255)
      cTTGCCGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAgg                                                                                                                                                >  5:295555/1‑147 (MQ=255)
      cTTGCCGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAgg                                                                                                                                                <  6:295555/147‑1 (MQ=255)
           cGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCcg                                                                                                                                         >  3:103601/1‑149 (MQ=255)
               cATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTg                                                                                                                                    >  3:341437/1‑149 (MQ=255)
               cATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTg                                                                                                                                    <  7:35697/149‑1 (MQ=255)
                aTCGCATCCACATCGCCGCTGCCCTTGAGCTGGGGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCCTCCTTTGACCCCTATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCCCTATTGAGCGCCcgcgca                                                                                                                                    <  1:281330/149‑3 (MQ=255)
                 tCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTgg                                                                                                                                   >  2:2918/1‑148 (MQ=255)
                   gCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCa                                                                                                                                 <  5:359173/148‑1 (MQ=255)
                      tCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGc                                                                                                                             >  4:123910/1‑149 (MQ=255)
                        cAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGAAGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGccc                                                                                                                           >  7:13983/1‑149 (MQ=255)
                         aGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCt                                                                                                                          <  7:265356/149‑1 (MQ=255)
                            tCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTa                                                                                                                       <  4:278320/149‑1 (MQ=255)
                            tCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTa                                                                                                                       <  1:2918/149‑1 (MQ=255)
                               ccgcTGCCCTTGAGCTGGAGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGa                                                                                                                     >  6:228069/1‑148 (MQ=255)
                                cgcTGCCCTTGAGCTGGTGATCTGGAGCACGCCTCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTc                                                                                                                                                      <  8:140405/115‑1 (MQ=255)
                                cgcTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTc                                                                                                                                                      >  7:140405/1‑115 (MQ=255)
                                cgcTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTc                                                                                                                                                      <  8:145428/115‑1 (MQ=255)
                                cgcTGCCCTTGAGCTGGTGATCTGGAGCACGCCCCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTc                                                                                                                                                      >  7:145428/1‑115 (MQ=255)
                                 gcTGCCCTTGAGCTGGTGAGCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAAc                                                                                                                  <  2:243438/149‑1 (MQ=255)
                                        ttGAGCTGCTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGctct                                                                                                           >  5:300855/1‑149 (MQ=255)
                                           aGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCa                                                                                                        >  8:331772/1‑149 (MQ=255)
                                              tGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCa                                                                                                     >  3:1849/1‑149 (MQ=255)
                                                gTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCAcc                                                                                                   >  4:302993/1‑149 (MQ=255)
                                                gTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGAAATGCCGTGTTGAAACGCGCTCTCCAGGAcc                                                                                                   >  1:217111/1‑149 (MQ=255)
                                                 tGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTACGCGGGCGCCCTATTGAGCGccc                                                                                                                                          >  6:13735/1‑110 (MQ=255)
                                                 tGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTACGCGGGCGCCCTATTGAGCGccc                                                                                                                                          <  5:13735/110‑1 (MQ=255)
                                                    tCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTATGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTCGCAATGACCTGTAGAAACGCGCTCTCCAGCACCGTTg                                                                                               >  7:167264/1‑149 (MQ=255)
                                                        gAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGcaca                                                                                           >  8:8914/1‑149 (MQ=255)
                                                              gccgcTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACActctct                                                                                     >  1:130691/1‑149 (MQ=255)
                                                              gccgcTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACActctct                                                                                     <  7:8914/149‑1 (MQ=255)
                                                               ccgcTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCtt                                                                                    <  4:1849/149‑1 (MQ=255)
                                                                cgcTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTg                                                                                   >  3:49426/1‑149 (MQ=255)
                                                                     gATGACGGCGAGCAATAGCGTTCTTTCACTTCCTTTGACGCAATATGTATTGACGGCTGTTGCCCTGGGGGGGCGCTCAATAGGGCTCCCGCGTTGCAATGGCCTGTAGAAACGCGCTCTCCTGCACGGTTGCACACTCTCTTGgcgc                                                                               <  8:167264/148‑1 (MQ=255)
                                                                     gATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGgcgc                                                                               <  3:302993/148‑1 (MQ=255)
                                                                      aTGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGgcgcgc                                                                             >  4:262770/1‑149 (MQ=255)
                                                                      aTGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGgcgcgc                                                                             >  2:24663/1‑149 (MQ=255)
                                                                           gCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGcccc                                                                        >  6:275842/1‑149 (MQ=255)
                                                                           gCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGcccc                                                                        >  8:284736/1‑149 (MQ=255)
                                                                             cGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGgcgc                                                                                                                                            >  4:284085/1‑80 (MQ=255)
                                                                               aGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGtt                                                                    <  4:292846/149‑1 (MQ=255)
                                                                                  aaTAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTg                                                                                   <  2:280673/131‑1 (MQ=255)
                                                                                  aaTAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTg                                                                                   >  1:280673/1‑131 (MQ=255)
                                                                                    tAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTggcgg                                                               <  4:341437/149‑1 (MQ=255)
                                                                                      gCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTGACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTggcggg                                                              <  6:300855/148‑1 (MQ=255)
                                                                                      gCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGACTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCAcc                                                                                                   <  1:1626/111‑1 (MQ=255)
                                                                                      gCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGACTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCAcc                                                                                                   >  2:1626/1‑111 (MQ=255)
                                                                                         ttCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTc                                                          <  2:130691/149‑1 (MQ=255)
                                                                                                 cttccttTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGAccatgc                                                  >  1:369459/1‑149 (MQ=255)
                                                                                                 cttccttTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGAccatgc                                                  <  4:49426/149‑1 (MQ=255)
                                                                                                  ttccttTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTg                                                                                   >  4:20667/1‑115 (MQ=255)
                                                                                                  ttccttTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGAccatgc                                                  >  2:193916/1‑148 (MQ=255)
                                                                                                    ccttTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGAccatgcca                                                <  5:228069/148‑1 (MQ=255)
                                                                                                             gCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGc                                       <  8:13983/149‑1 (MQ=255)
                                                                                                             gCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGc                                       <  7:331772/149‑1 (MQ=255)
                                                                                                              cAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCg                                      <  5:275842/149‑1 (MQ=255)
                                                                                                               aaTAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGa                                     >  3:123652/1‑149 (MQ=255)
                                                                                                                    gTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACa                                         <  4:34858/140‑1 (MQ=255)
                                                                                                                    gTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACa                                         >  3:34858/1‑140 (MQ=255)
                                                                                                                    gTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACg                                 >  4:113884/1‑148 (MQ=255)
                                                                                                                       ttCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACa                             >  2:301285/1‑149 (MQ=255)
                                                                                                                            ggCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCg                        <  3:262770/149‑1 (MQ=255)
                                                                                                                              cTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAAc                                                                                                                  >  7:152028/1‑57 (MQ=255)
                                                                                                                              cTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAAc                                                                                                                  <  8:152028/57‑1 (MQ=255)
                                                                                                                              cTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGt                      >  5:34760/1‑149 (MQ=255)
                                                                                                                                 ttGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCa                    <  3:123910/148‑1 (MQ=255)
                                                                                                                                 ttGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCa                    <  1:261769/148‑1 (MQ=255)
                                                                                                                                      cTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAAcg              >  4:140665/1‑149 (MQ=255)
                                                                                                                                      cTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACGCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAAc                                          >  4:256558/1‑121 (MQ=255)
                                                                                                                                      cTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACCCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAAc                                          <  3:256558/121‑1 (MQ=255)
                                                                                                                                       cacgcgGGCGCCCTATTGAGCGCCCGCGGCGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAAcg              >  1:22374/3‑149 (MQ=255)
                                                                                                                                       cacgcgGGCGCCCTATTGAGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAAcgc             >  5:11351/3‑149 (MQ=255)
                                                                                                                                         cgcgGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAAcgc             <  7:284736/147‑1 (MQ=255)
                                                                                                                                                gcTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAACGCCGCGGTGa     >  7:91256/1‑148 (MQ=255)
                                                                                                                                                  tCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAACGCCGCGGTGACAt  >  7:243745/1‑149 (MQ=255)
                                                                                                                                                  tCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAACGCCGCGGTGACAt  >  6:79710/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCCCGCTTGCCGACCATCGCATCCAGATCGCCGCTGCCCTTGAGCTGGTGATCTGGAGCACGCCGCTCGATGACGCCGAGCAATAGCGTTCATTCCCTTCCTTTGACGCAATAGGTCTTCACGGCTGTTGCCCTGCGCGGGCGCTCAATAGGGCGCCCGCGTGGCAATGCCCTGTAGAAACGCGCTCTCCAGCACCGTTGCACACTCTCTTGGCGCGCGCCCCCGTTGGCGGGCGTCGACCATGCCATGCAACAGCGATACGAACATCTCGGTCAGCAACGCCGCGGTGACAT  >  NC_002947/3212131‑3212424

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: