Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,538,533 T→A 5.2% intergenic (+23/+98) PP_4871 → / ← PP_4872 lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,538,5330TA 5.2% 222.2 / 2.7 77intergenic (+23/+98)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base A (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

TTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAC  >  NC_002947/5538387‑5538677
                                                                                                                                                  |                                                                                                                                                
tttCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTcg                                                                                                                                                >  4:71375/1‑149 (MQ=255)
tttCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTcg                                                                                                                                                >  7:220377/1‑149 (MQ=255)
   cGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCgg                                                                                                                                             >  3:250916/1‑149 (MQ=255)
     tGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTc                                                                                                                                           <  3:141588/149‑1 (MQ=255)
       cgcCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTcc                                                                                                                                          <  7:58595/148‑1 (MQ=255)
       cgcCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCGAAGCGGTccc                                                                                                                                         <  5:28202/149‑2 (MQ=255)
         cccGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCttt                                                                                                                                       >  1:228906/1‑149 (MQ=255)
          ccGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTc                                                                                                                                      <  1:350502/149‑1 (MQ=255)
              aaCACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTcg                                                                                                                                     >  2:335524/1‑146 (MQ=255)
              aaCACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTcg                                                                                                                                     <  1:335524/146‑1 (MQ=255)
              aaCACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCATGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGa                                                                                                                                  <  1:82345/149‑1 (MQ=255)
                 acCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGaca                                                                                                                                >  4:34487/1‑148 (MQ=255)
                 acCGGGCCATGAAGCTGCTGCGGCAGCAGCCGGAAGAGCTGCTGGAAGGGATGGACAGCATTGTGGCGCCCGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTAAGACTATTGGGACCGCTTCGCGGTCCATTAGCGaca                                                                                                                                >  8:217186/1‑148 (MQ=255)
                    gggCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAg                                                                                                                            <  6:309866/149‑1 (MQ=255)
                     ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAgg                                                                                                                           <  3:309154/149‑1 (MQ=255)
                     ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAgg                                                                                                                           <  2:228906/149‑1 (MQ=255)
                     ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAg                                                                                                                            >  1:82430/1‑148 (MQ=255)
                     ggCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAg                                                                                                                            <  2:82430/148‑1 (MQ=255)
                        gATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCg                                                                                                                        >  1:195963/1‑149 (MQ=255)
                         gtccggctgctGCCCCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGc                                                                                                                       >  6:181493/6‑149 (MQ=255)
                            ttGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTc                                                                                                                     <  3:199613/148‑1 (MQ=255)
                              gctgctGGGGCATCAGCCGGACGATCTTCTGGTATGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAATTTGGTTGTCAATCATCCTTACTTAATCTGTTTTTTTTCGGACCGCTTCGCGGTCCTTTcg                                                                                                                                     <  7:235117/130‑1 (MQ=255)
                              gctgctGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTcg                                                                                                                                     >  8:235117/1‑130 (MQ=255)
                              gctgctGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCca                                                                                                                  >  1:25633/1‑149 (MQ=255)
                                tgctgGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCcac                                                                                                                 >  2:17432/1‑148 (MQ=255)
                                     gggCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAg                                                                                                             <  1:17432/147‑1 (MQ=255)
                                            agcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                    <  4:268459/149‑1 (MQ=255)
                                              ccGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTTGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTc                                                                                                                                           >  2:118808/1‑108 (MQ=255)
                                              ccGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                    >  8:316485/1‑147 (MQ=255)
                                              ccGGACGAGCTGCTGGAAGGGATGGACAGCATTGTGGCGCCGGTTGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTc                                                                                                                                           <  1:118808/108‑1 (MQ=255)
                                               cGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAg                                                                                                 >  6:158947/1‑149 (MQ=255)
                                                ggACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGc                                                                                                >  8:318301/1‑149 (MQ=255)
                                                   cGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGc                                                                                             <  1:215562/149‑1 (MQ=255)
                                                   cGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGc                                                                                             <  2:152931/149‑1 (MQ=255)
                                                      gctgctGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                          <  3:71375/149‑1 (MQ=255)
                                                        tgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCGCCCACAAAGGCAATGGGGGAGCGGCCgtgt                                                                                        >  3:21220/1‑149 (MQ=255)
                                                        ggctgGAAGGGATGGACAGCTTTGTGGCGCCGGGCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtgt                                                                                        <  3:10886/148‑1 (MQ=255)
                                                          ctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtgt                                                                                        <  2:220426/147‑1 (MQ=255)
                                                             gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                   >  1:165627/1‑149 (MQ=255)
                                                             gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTc                                                                                     <  6:137311/147‑1 (MQ=255)
                                                                ggtaTGGACAGCTGTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACTAGCAGCCTGACTTAAGCTTAGACTTGTCGGACGGCTTCGCGGTCCTTTCGCTTCACAAGGCCGCTCCCACAAAGGCATTGGGGCAGCGGCCTTGTGTCGCGAt                                                                                <  4:213231/146‑1 (MQ=255)
                                                                gggATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGAt                                                                                >  3:244992/1‑149 (MQ=255)
                                                                     ggACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCt                                                                           <  5:138116/149‑1 (MQ=255)
                                                                      gACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCt                                                                           <  1:62279/148‑1 (MQ=255)
                                                                      gACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                                          <  7:318301/149‑1 (MQ=255)
                                                                      gACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTGCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTg                                                                          <  7:333401/149‑1 (MQ=255)
                                                                       aCAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAgg                                                                                                            <  3:337149/114‑1 (MQ=255)
                                                                       aCAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAgg                                                                                                            >  4:337149/1‑114 (MQ=255)
                                                                          gCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCa                                                                        >  6:281264/1‑147 (MQ=255)
                                                                                    gcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCTTTGCAGcccatc                                                            <  4:244992/149‑4 (MQ=255)
                                                                                    gcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGATGAGACTATTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGccccc                                                             >  1:72555/1‑147 (MQ=255)
                                                                                      cGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATa                                                          <  5:281264/149‑1 (MQ=255)
                                                                                        gTCGCGCCTACGGGGGCCCACACGCAGCCTCACTCACGCTGAGACTTTTGGCCCCGCTTCGCGGCCCTTTCGCCACACACGCCCGCCCCCACAACGCCAATGGGGGCGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGt                                                        <  2:195963/149‑1 (MQ=25)
                                                                                        gTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACGGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGt                                                        <  7:316485/149‑1 (MQ=255)
                                                                                            cgcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTACTTTCGCGACACAAGGCCGCACCCACAAAGGAAATGGGGGAGCGGCCTTGAGTCGCGATGGGCTGCAAAGCACCCCCAATAGTCGa                                                     >  7:57462/1‑148 (MQ=255)
                                                                                              ccTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCCAGAGTCGATCa                                                  >  6:234871/1‑149 (MQ=255)
                                                                                                   gTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                          >  8:86866/1‑104 (MQ=255)
                                                                                                   gTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                          <  7:86866/104‑1 (MQ=255)
                                                                                                           aCAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAAt                                     <  2:25633/149‑1 (MQ=255)
                                                                                                                  gcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCCATAGTCGCTCAGTTAACGCGGAAGGACAGGc                              >  4:88478/1‑149 (MQ=255)
                                                                                                                       cattaaGGTGAGACTTTTGGGACCGGTTCGCGGTCCTTTCGCGAAACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGa                           <  3:34487/143‑1 (MQ=255)
                                                                                                                       actgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGAcg                         >  8:86370/1‑149 (MQ=255)
                                                                                                                         tgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGAcgcg                       >  5:147775/1‑149 (MQ=255)
                                                                                                                           aaGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGc                     >  1:280812/1‑149 (MQ=255)
                                                                                                                               tGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttct                 >  6:371301/1‑149 (MQ=255)
                                                                                                                               tGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCCCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGctact                 >  2:202298/1‑149 (MQ=255)
                                                                                                                                 agaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttg               <  7:194846/149‑1 (MQ=255)
                                                                                                                                  gaCTTTTGGGACCGCTTCGCGGTCCTTTCGGGTCACAATGCCGGTCCCCCCAAGGCAATGGGGGTGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttg               <  4:250916/148‑1 (MQ=38)
                                                                                                                                    gcgaaaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgttg            >  7:213225/7‑149 (MQ=39)
                                                                                                                                     ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGc           >  6:163616/1‑149 (MQ=255)
                                                                                                                                       ttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGGCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGcc          <  2:280812/148‑1 (MQ=255)
                                                                                                                                        aaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgtt             <  7:271514/142‑1 (MQ=255)
                                                                                                                                        aaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgtt             >  8:271514/3‑144 (MQ=255)
                                                                                                                                           gACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGc     >  5:191490/1‑149 (MQ=255)
                                                                                                                                            aCCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGGAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGcc    >  1:263058/1‑149 (MQ=255)
                                                                                                                                            aCCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGc     >  1:198089/1‑148 (MQ=255)
                                                                                                                                              cGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAc  <  6:147775/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAC  >  NC_002947/5538387‑5538677

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: