Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,183 G→A 9.9% intergenic (+44/‑155) PP_0663 → / → PP_0664 AsnC family transcriptional regulator/homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947772,1830GA 9.9% 219.6 / 9.6 80intergenic (+44/‑155)PP_0663/PP_0664AsnC family transcriptional regulator/homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base A (4/4);  total (40/40)
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.33e-01

CTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAGGACA  >  NC_002947/772040‑772330
                                                                                                                                                |                                                                                                                                                      
cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                                   >  2:17081/1‑149 (MQ=255)
        cgtaGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATTCCCGGGGGGGAGAACTCGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTTGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGtt                                                                                                                                           <  6:39668/146‑1 (MQ=38)
              gggTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGt                                                                                                                                   >  6:242119/1‑149 (MQ=255)
                 tGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGtttt                                                                                                                                <  4:77386/149‑1 (MQ=255)
                 tGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGtttt                                                                                                                                >  7:124283/1‑149 (MQ=255)
                   cACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGATCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTg                                                                                                                               <  1:29712/148‑1 (MQ=255)
                    aCCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGcc                                                                                                                             >  4:95420/1‑149 (MQ=255)
                    aCCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGcc                                                                                                                             >  5:292513/1‑149 (MQ=255)
                       aTGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTg                                                                                                                          >  5:125085/1‑149 (MQ=255)
                             ggaaggaaATCGCTGCCATGGCCGGGGTGGAGAACACGGTGTCTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGGTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAg                                                                                                                   <  5:306164/149‑1 (MQ=255)
                              gaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAgg                                                                                                                  >  4:258085/1‑149 (MQ=255)
                              gaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAgg                                                                                                                  >  2:77744/1‑149 (MQ=255)
                              gaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAgg                                                                                                                  >  5:56885/1‑149 (MQ=255)
                               aaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGt                                                                                                                 >  6:8845/1‑149 (MQ=255)
                               aaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGt                                                                                                                 >  1:140743/1‑149 (MQ=255)
                                aggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGT‑‑Gctagc                                                                                                                                 >  7:293329/1‑129 (MQ=38)
                                aggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGT‑‑Gctagc                                                                                                                                 <  8:293329/133‑5 (MQ=38)
                                  gaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTgcg                                                                                                              >  6:50350/1‑149 (MQ=255)
                                  gaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTgcg                                                                                                              <  5:242119/149‑1 (MQ=255)
                                  gaaATCGCTGCCATGCCCGGGGGGGAGAACACGGTGACTTGGTGCGTCCTGGCCGCGGTGGGCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTgcg                                                                                                              <  6:51623/149‑1 (MQ=255)
                                        gCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                        <  6:125085/149‑1 (MQ=255)
                                            ccATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGTGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCAc                                                                                                    >  2:84081/1‑149 (MQ=255)
                                                  ccGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGc                                                                                                >  4:88027/1‑147 (MQ=255)
                                                      ggTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGAtt                                                                                          <  2:140743/149‑1 (MQ=255)
                                                        tGGAGAAGACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTgg                                                                                        <  6:292513/149‑1 (MQ=255)
                                                        tGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTgg                                                                                        <  3:95420/149‑1 (MQ=255)
                                                        tGGAGAACACGGGGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTgg                                                                                        <  8:157567/149‑1 (MQ=255)
                                                          gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                       >  1:178188/1‑148 (MQ=255)
                                                          gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                       <  3:88027/148‑1 (MQ=255)
                                                          gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                       <  5:50350/148‑1 (MQ=255)
                                                            gaACACGGTGACTTCGTTCGTCCTGGCCGGGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaa                                                                                    <  3:60632/149‑1 (MQ=255)
                                                            gaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaa                                                                                    >  6:248723/1‑149 (MQ=255)
                                                             aaCACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaat                                                                                   <  5:114604/149‑1 (MQ=255)
                                                               cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaat                                                                                   >  4:61770/1‑147 (MQ=255)
                                                               cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaat                                                                                   <  3:61770/147‑1 (MQ=255)
                                                                acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcggc                                                                                                                                              <  7:158469/90‑1 (MQ=255)
                                                                acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcggc                                                                                                                                              >  8:158469/1‑90 (MQ=255)
                                                                acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                                <  2:178188/149‑1 (MQ=255)
                                                                acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                                >  6:18115/1‑149 (MQ=255)
                                                                  ggTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaata                                                                              >  8:271508/1‑149 (MQ=255)
                                                                   gTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATAc                                                                             <  5:8845/149‑1 (MQ=255)
                                                                      aCTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAg                                                                           <  1:77744/148‑1 (MQ=255)
                                                                      aCTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGc                                                                          <  5:18115/149‑1 (MQ=255)
                                                                            ttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                    >  7:65276/1‑149 (MQ=255)
                                                                               gtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCa                                                                 <  1:215446/149‑1 (MQ=255)
                                                                                 ccTGGCCGCGGTGGTCTGAACCGCGTCTTTGCAGTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                       >  1:197242/1‑142 (MQ=255)
                                                                                 ccTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                       <  2:197242/141‑1 (MQ=255)
                                                                                    ggCCGGGGTGGTCTGAACCGCGTGTTTGCATTGGCTATCACAACGCCGGCTGGGTTCGCGCCCAGGCGGTGTTGT‑‑GCGTTTTGCCCTGGCGCACGGTGCGCGCTTTGCACTGGCCAAATTGGCCAATAAATACAGCGGATTGCCATAtt                                                             <  1:84081/149‑1 (MQ=255)
                                                                                    ggCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAtt                                                             >  5:152554/1‑148 (MQ=255)
                                                                                     gCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTc                                                            <  3:258085/148‑1 (MQ=255)
                                                                                       cgcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAAc                                                         <  6:152554/149‑1 (MQ=255)
                                                                                           gtggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTgg                                                      >  4:136691/1‑148 (MQ=255)
                                                                                                cTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGAGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGt                                                  >  1:91961/1‑149 (MQ=255)
                                                                                                 tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagt                                                                                   >  5:242093/1‑111 (MQ=255)
                                                                                                 tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagt                                                                                   <  6:242093/113‑3 (MQ=255)
                                                                                                   aaCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                              <  5:248723/148‑1 (MQ=255)
                                                                                                      cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGct                                          >  2:73729/1‑149 (MQ=255)
                                                                                                      cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCAGATTCAACTGGCAGTGCACGGct                                          >  8:208862/1‑149 (MQ=255)
                                                                                                       gcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGctc                                         <  1:73729/149‑1 (MQ=255)
                                                                                                           cTTTGCATGGGCTAGCACAACGCCGCCGGCGTGCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTAc                                     <  8:124850/149‑1 (MQ=255)
                                                                                                             ttGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGc                                   >  8:167003/1‑149 (MQ=255)
                                                                                                              tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTAc                                     <  1:17081/146‑1 (MQ=255)
                                                                                                                cATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                                 <  7:271508/148‑1 (MQ=255)
                                                                                                                   ttGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTc                              <  6:56885/149‑1 (MQ=255)
                                                                                                                      gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCt                            >  7:141520/1‑148 (MQ=255)
                                                                                                                      gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTc                           >  3:62274/1‑149 (MQ=255)
                                                                                                                        acgCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGA‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCATTGCACTCGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCATCTGGCAGTGCACGGATCTTACGCTAGTCCTCAt                         >  4:241710/3‑149 (MQ=38)
                                                                                                                           cacaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGa                      >  8:33256/1‑149 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                            >  6:177501/1‑124 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                            <  5:177501/124‑1 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                    >  1:92845/1‑148 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGCGCGAACCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCAt                         >  3:30448/1‑143 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGCGCGAACCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCAt                         <  4:30448/143‑1 (MQ=255)
                                                                                                                                cgccgcCTGGGTTCGCGCCCAGGCGGCGTTTT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa                  <  1:274415/148‑1 (MQ=255)
                                                                                                                                cgccgcCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa                  >  2:274415/1‑148 (MQ=255)
                                                                                                                                  ccgccTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCAGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg               >  1:104706/1‑149 (MQ=255)
                                                                                                                                   cgccTGGGCGCGAACCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGt                                                  >  1:276684/1‑113 (MQ=39)
                                                                                                                                   cgccTGGGCGCGAACCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGt                                                  <  2:276684/113‑1 (MQ=39)
                                                                                                                                         ggTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTgga         <  7:167003/148‑1 (MQ=255)
                                                                                                                                           ttCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacag      <  7:33256/149‑1 (MQ=255)
                                                                                                                                               cgcCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacaggaca  <  3:136691/149‑1 (MQ=255)
                                                                                                                                                gcCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg               >  3:199514/1‑135 (MQ=255)
                                                                                                                                                gcCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg               <  4:199514/135‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGT‑‑GCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAGGACA  >  NC_002947/772040‑772330

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: