Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,810,289 Δ1 bp 6.3% intergenic (+50/‑119) eno → / → ftsB enolase/cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,810,2890T. 6.3% 195.5 / 7.2 63intergenic (+50/‑119)eno/ftsBenolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base . (2/2);  total (30/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 = 7.34e-01

CGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTT  >  NC_002947/1810146‑1810437
                                                                                                                                                 |                                                                                                                                                        
cGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACt                                                                                                                                                     >  2:177386/1‑149 (MQ=255)
 gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACtt                                                                                                                                                    <  1:177386/149‑1 (MQ=255)
 gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACtt                                                                                                                                                    >  6:121638/1‑149 (MQ=255)
 gACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACtt                                                                                                                                                    <  1:14215/149‑1 (MQ=255)
 gACCGCGTTTCCAAGTACAAACAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGAAGCCGTAGTTGTAGGAACGTTCGTACTGAACtt                                                                                                                                                    >  2:24180/1‑149 (MQ=255)
   ccGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTc                                                                                                                                                  <  2:287888/149‑1 (MQ=255)
     gcgTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCt                                                                                                                                                <  1:191949/149‑1 (MQ=255)
           ccAAGTACAACCCGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGGGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTa                                                                                                                                           <  6:204847/148‑1 (MQ=255)
              aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTTCAGTACGAACGTTCCTAcaactccg                                                                                                                                       <  6:252081/149‑8 (MQ=38)
              aGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAAAGTTCAGTACGAACGTTCCTAcaactccg                                                                                                                                       <  6:246555/149‑8 (MQ=38)
                     ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                               >  2:60207/1‑149 (MQ=255)
                     ccAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGttt                                                                                                                               <  2:332863/149‑1 (MQ=255)
                          tgctgcGCATCGAAGAGGAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGGTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGt                                                                                                                           <  1:60207/149‑1 (MQ=255)
                                      aaGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTgg                                                                                                            >  1:23778/1‑149 (MQ=255)
                                         agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                         >  5:379566/1‑149 (MQ=255)
                                             aCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                     <  5:93290/149‑1 (MQ=255)
                                               tGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACgg                                                                                                                                    <  8:241744/119‑1 (MQ=255)
                                               tGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACgg                                                                                                                                    >  7:241744/1‑119 (MQ=255)
                                                 ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGGTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                                 >  1:334457/1‑149 (MQ=255)
                                                   cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                               <  5:121638/149‑1 (MQ=255)
                                                    gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                                 <  7:180687/146‑1 (MQ=255)
                                                    gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCACGGCTAAGCAAGTGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                                 >  8:180687/1‑146 (MQ=255)
                                                      cAAAGCGGTTTACCGCTGTCGTGCCTAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTCg                                                                                            <  5:273526/149‑1 (MQ=255)
                                                       aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt                                                                                           <  6:379566/149‑1 (MQ=255)
                                                              tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                    >  5:150940/1‑149 (MQ=255)
                                                              tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                    >  2:37959/1‑149 (MQ=255)
                                                              tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTc                                                                                    >  2:43506/1‑149 (MQ=255)
                                                                tACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCg                                                                                   <  5:343471/148‑1 (MQ=255)
                                                                  ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCt                                                                                                              <  2:41648/119‑1 (MQ=255)
                                                                  ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCt                                                                                                              >  1:41648/1‑119 (MQ=255)
                                                                   cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGATATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGTCTTCGTg                                                                                          <  7:185919/138‑1 (MQ=255)
                                                                   cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCCGTCGACGAAGCTCTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                                <  1:37959/149‑1 (MQ=255)
                                                                   cgcgGTCGTGCCGAGGTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGGTTCTGTCGACGAAGCCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAgg                                                                                <  1:43506/149‑1 (MQ=255)
                                                                          gTGCCGAGTTTCGCGGCTAAGCAAGAGATCGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGATAGTACTCAACTTTCCTACGTACAAACGCGTTTCTGTCGACGAAGCCTCGCCACGGCCAGGCATCGTGCTAATCGAGGCCCCGaac                                                                        >  6:364314/1‑148 (MQ=255)
                                                                           tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                       <  8:182581/149‑1 (MQ=255)
                                                                           tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                       <  2:334457/149‑1 (MQ=255)
                                                                             ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                                      >  7:116276/1‑148 (MQ=255)
                                                                                        ggCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                          <  8:116276/149‑1 (MQ=255)
                                                                                          cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                          <  4:98909/147‑1 (MQ=255)
                                                                                          cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                          >  3:98909/1‑147 (MQ=255)
                                                                                              gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCt                                                   >  5:305889/1‑149 (MQ=255)
                                                                                                     aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc                                            >  6:353495/1‑149 (MQ=255)
                                                                                                     aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTGCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc                                            <  6:211271/149‑1 (MQ=255)
                                                                                                        gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTTCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCt                                          <  2:312350/148‑1 (MQ=255)
                                                                                                        gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                     <  4:270199/90‑1 (MQ=255)
                                                                                                        gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                     >  3:270199/1‑90 (MQ=255)
                                                                                                            aaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACGAACGGGTTTCTGTCGACGACGCCTGGGCGCGGCCAGGGTGCGTGGTATTCGAGGCCCCGAAGTAACAAGTACGCGGGCGCCGTTTTAACCTGGAt                                      >  5:195073/1‑148 (MQ=255)
                                                                                                             aaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATaa                                    <  6:150940/149‑1 (MQ=255)
                                                                                                              aGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAAc                                   <  2:23778/149‑1 (MQ=255)
                                                                                                                  agcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                               >  1:202343/1‑149 (MQ=255)
                                                                                                                       ccGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCt                                          <  7:318790/133‑1 (MQ=255)
                                                                                                                       ccGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTAGTGCTATTAGAGGCCCCGAAATAACAAGTAAGCGGCAGCCTTTTTAAACt                                          >  8:318790/1‑133 (MQ=255)
                                                                                                                         ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                         <  7:137692/132‑1 (MQ=255)
                                                                                                                         ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCACGGCCAGGCATAGTGGGATTAGAGGCCCCGAAATAAAAAATAACCGGCAACCATTTTTAAcag                                          >  8:137692/1‑130 (MQ=255)
                                                                                                                            tacgTAGGAAAGTTCAGTACGAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTcc                     >  7:81272/4‑149 (MQ=255)
                                                                                                                                tAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAATTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                 >  6:282951/1‑149 (MQ=255)
                                                                                                                                tAGGAAAGTTCAGTACGAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                 >  5:110353/1‑149 (MQ=255)
                                                                                                                                   gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg              >  7:209220/1‑149 (MQ=255)
                                                                                                                                   gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg              <  6:305889/149‑1 (MQ=255)
                                                                                                                                   gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGACGCGCAGTCCCTATTgg              >  2:244930/1‑149 (MQ=255)
                                                                                                                                           cGTACTGAACTTTCATACTTACTAACTGGTTTCGGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTc        <  6:195073/148‑1 (MQ=255)
                                                                                                                                              tACTGAACTTTCCTACGTACTAACGGGTGTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGAGCAGTCCCTATTGGTTGTTCCt      >  1:22048/1‑149 (MQ=255)
                                                                                                                                                cTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGTGCAGTCCCTATTGGTTGTTCCTCGt   >  7:37827/1‑149 (MQ=255)
                                                                                                                                                 tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt  >  3:140459/1‑149 (MQ=255)
                                                                                                                                                 tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGACGCGCAGTCCCTATTGGTTGTTCCTCGtt  <  1:244930/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
CGACCGCGTTTCCAAGTACAACCAGCTGCTGCGCATCGAAGAGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTT  >  NC_002947/1810146‑1810437

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: