Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,528,353 T→A 6.6% intergenic (+28/‑48) purB → / → PP_4017 adenylosuccinate lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,528,3530TA 6.6% 202.9 / 2.2 70intergenic (+28/‑48)purB/PP_4017adenylosuccinate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base A (3/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGC  >  NC_002947/4528209‑4528499
                                                                                                                                                 |                                                                                                                                                  
gcgcTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAACGCCCGGCCAGGcc                                                                                                                                                <  7:4493/149‑1 (MQ=255)
      ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGc                                                                                                                                            <  2:65315/147‑1 (MQ=255)
      ctgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGt                                                                                                                                          <  1:113388/149‑1 (MQ=255)
           cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                     >  4:51452/1‑149 (MQ=255)
           cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTa                                                                                                                                     >  2:141684/1‑149 (MQ=255)
                  gACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCg                                                                                                                               <  5:202343/148‑1 (MQ=255)
                    cGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc                                                                                                                            >  5:86683/1‑149 (MQ=255)
                        cTGGACATGCCAGCCGAAGCCAAGGCCGATCTCAAGCAGCTCACCCCCTCAACCTACATCGGTAACGCGGCAGCCCAGACCAAACGCATCTAAGCTCAGCGTCGCGTACAACGCCCGGCCTCGCCGCGCGTTATTAAACACCTACGaa                                                                                                                         >  6:157390/1‑148 (MQ=255)
                                     ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                            >  6:197331/1‑148 (MQ=255)
                                     ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                           <  6:294420/149‑1 (MQ=255)
                                      cGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCTCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCa                                                                                                          <  4:115014/149‑1 (MQ=255)
                                       gAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCaa                                                                                                         <  2:214480/149‑1 (MQ=255)
                                       gAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCaa                                                                                                         <  5:197331/149‑1 (MQ=255)
                                        aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCCACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                            >  7:230846/1‑145 (MQ=255)
                                        aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCCACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt                                                                                                            <  8:230846/145‑1 (MQ=255)
                                           ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                     >  6:313528/1‑149 (MQ=255)
                                             aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt                                                                                                   >  5:290310/1‑149 (MQ=255)
                                             aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTACCGTACGAAAAATACGTTTTTTCAATAGGtt                                                                                                   >  8:272687/1‑149 (MQ=255)
                                              aGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCAGCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTg                                                                                                  <  4:294048/149‑1 (MQ=255)
                                                  cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCggg                                                                                                                                             <  5:33679/102‑1 (MQ=255)
                                                  cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCggg                                                                                                                                             >  6:33679/1‑102 (MQ=255)
                                                   gAGCTCAAGCAGCTTACCCCCGCTACCTACATCGTTAACGCGGCAGCCCAGACCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt                                                                                             <  5:157390/149‑1 (MQ=255)
                                                   gAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACg                                                                                                                  >  5:9387/1‑128 (MQ=255)
                                                   gAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACg                                                                                                                  <  6:9387/128‑1 (MQ=255)
                                                   gAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt                                                                                             <  1:141684/149‑1 (MQ=255)
                                                      cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGaa                                                                                          <  3:51452/149‑1 (MQ=255)
                                                       tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAAt                                                                                         >  6:297156/1‑149 (MQ=255)
                                                       tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAAt                                                                                         >  7:55913/1‑149 (MQ=255)
                                                         aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATc                                                                                        <  8:71246/148‑1 (MQ=255)
                                                         aaGCAGCTGACCCCCGCTAACTACATCGGTAACGCGGCCTATCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATc                                                                                        <  1:105019/149‑1 (MQ=255)
                                                              gcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATa                                                                                  <  6:86683/149‑1 (MQ=255)
                                                                  aCCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTc                                                                               >  6:184749/1‑148 (MQ=255)
                                                                    ccccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCAc                                                                            >  1:208785/1‑149 (MQ=255)
                                                                     cccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACt                                                                           <  5:297156/149‑1 (MQ=255)
                                                                      ccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTGCACTg                                                                          >  3:167242/1‑149 (MQ=255)
                                                                            cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTGCACTGCAgc                                                                      <  4:167242/147‑1 (MQ=255)
                                                                            cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgctg                                                                    >  2:13853/1‑149 (MQ=255)
                                                                                  tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc                                                              >  7:37849/1‑149 (MQ=255)
                                                                                  tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc                                                              <  6:290310/149‑1 (MQ=255)
                                                                                    ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAGAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcgg                                                            >  8:70539/1‑149 (MQ=255)
                                                                                    ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcg                                                             <  8:55913/148‑1 (MQ=255)
                                                                                    ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcg                                                             <  1:13853/148‑1 (MQ=255)
                                                                                     gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAGAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCggcggc                                                           <  7:70539/149‑1 (MQ=255)
                                                                                     gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTTCTcggcgg                                                            >  6:291169/1‑148 (MQ=255)
                                                                                             gCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCAc                                                                            <  1:80690/124‑1 (MQ=255)
                                                                                             gCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCAc                                                                            >  2:80690/1‑124 (MQ=255)
                                                                                             gCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGc                                                   <  7:155965/149‑1 (MQ=255)
                                                                                                  cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgc                                               >  7:165819/1‑149 (MQ=255)
                                                                                                  cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCAGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgc                                               >  8:54685/1‑149 (MQ=255)
                                                                                                    cAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACTCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAAt                                                                                                        <  1:107086/90‑1 (MQ=255)
                                                                                                            aggCATCTATGCTGAGCGTCGCGTACAAGGCCCGGCCTGGCGGGGCGTTTTTATTCCCGTTCGAAAATTACGTTTTTTCAATAGGTTGAACTTGAATCCTGATACTCCACTGCAGCTGCTCGGCTGCATCTCGGCGCGCGAATTCATgc                                     <  8:165819/147‑1 (MQ=255)
                                                                                                                aTCTAAGCTGAGCGTCGCGTACAACGCACGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGa                                 >  5:312862/1‑149 (MQ=255)
                                                                                                                       cTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATTGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTgg                          >  2:80603/1‑149 (MQ=255)
                                                                                                                       cTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                     >  1:35789/1‑138 (MQ=255)
                                                                                                                       cTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgc                                     <  2:35789/138‑1 (MQ=255)
                                                                                                                       cTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTgg                          >  2:84764/1‑149 (MQ=255)
                                                                                                                        tGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGa                                             <  1:104839/129‑1 (MQ=255)
                                                                                                                        tGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGa                                             >  2:104839/1‑129 (MQ=255)
                                                                                                                            cGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaa                     <  1:80603/149‑1 (MQ=255)
                                                                                                                            cGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaa                     <  1:84764/149‑1 (MQ=255)
                                                                                                                             gTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaag                    <  6:312862/149‑1 (MQ=255)
                                                                                                                                 cgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGc                <  5:184749/149‑1 (MQ=255)
                                                                                                                                     aaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgct            >  7:14814/2‑149 (MQ=255)
                                                                                                                                     aaaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGCCATCTCGGCCCGCGAATTAATGCGCGACTACTGGCAGAAGAAGCCgct            >  1:155613/2‑149 (MQ=255)
                                                                                                                                      aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTGCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGCCTACTGGCAGAAGAAGCCgctg           >  6:202731/1‑149 (MQ=255)
                                                                                                                                      aaCGCCCGGCCAGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAAGATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCaga                      >  1:124269/1‑138 (MQ=255)
                                                                                                                                        cGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGGTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTAGTGGCAGAAGAAGCCgctgct         >  2:129961/1‑149 (MQ=255)
                                                                                                                                           ccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCTTGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGt      <  5:313528/149‑1 (MQ=255)
                                                                                                                                            cGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTa     <  2:253325/149‑1 (MQ=255)
                                                                                                                                                cTGGCAGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAAAATGAATCCTTATACTACACTGCAGCTGCTCGGCGGCAGCTCGGCCCGCGAATTCATGCGCGACTAAAGGAAGAAGAACCAGCTGCTGCTaagg  >  4:210027/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGC  >  NC_002947/4528209‑4528499

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: