Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,070,749 T→C 7.9% intergenic (+317/‑37) aldB‑II → / → pqqD‑II aldehyde dehydrogenase/coenzyme PQQ synthesis protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,070,7490TC 7.9% 174.5 / 4.8 63intergenic (+317/‑37)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (3/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAG  >  NC_002947/3070601‑3070896
                                                                                                                                                    |                                                                                                                                                   
gTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCAt                                                                                                                                                     >  3:327453/1‑149 (MQ=255)
  ccGGCCCCGTGAATTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAAGCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACTATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATc                                                                                                                                                    <  4:81560/148‑1 (MQ=255)
               ttAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                      <  1:198232/149‑1 (MQ=255)
                tAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                      <  1:66546/148‑1 (MQ=255)
                tAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                      >  2:66546/1‑148 (MQ=255)
                      ttCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGAtcctcc                                                                                                                               >  6:337757/1‑149 (MQ=255)
                      ttCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGAtcctcc                                                                                                                               >  7:267103/1‑149 (MQ=255)
                          aCCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTg                                                                                                                             >  7:357505/1‑147 (MQ=255)
                          aCCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTg                                                                                                                             >  2:343019/1‑147 (MQ=255)
                            cGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATTACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGt                                                                                                                         >  7:79114/1‑149 (MQ=255)
                            cGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAg                                                                                                                          >  4:93613/1‑148 (MQ=255)
                            cGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGt                                                                                                                         <  8:357505/149‑1 (MQ=255)
                                      cgggtcaAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCtgt                                                                                                               <  6:64282/144‑1 (MQ=255)
                                       ggCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCtgtg                                                                                                              <  4:327453/149‑1 (MQ=255)
                                        gCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCt                                                                                                                 >  3:368500/1‑145 (MQ=255)
                                        gCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCt                                                                                                                 <  4:368500/145‑1 (MQ=255)
                                          acaAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAAc                                                                                                           >  2:146007/1‑149 (MQ=255)
                                             agaagaGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                      >  8:197666/1‑119 (MQ=255)
                                              gaagaGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCtga                                                                                                       <  8:362940/149‑1 (MQ=255)
                                              gaagaGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCtga                                                                                                       <  5:337757/149‑1 (MQ=255)
                                               aagaGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCtgat                                                                                                      <  2:161691/149‑1 (MQ=255)
                                                   ggACTTGCCCTATGTGTACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCtgattgat                                                                                                  <  5:335323/149‑1 (MQ=255)
                                                      cTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGCTGACGGCCCGGCCATCGCCACGGTGGCCGGATCCGCCTGGCGGAGGCGCCTGGGAACCGGCTGGATCGc                                                                                               >  8:81153/1‑149 (MQ=255)
                                                          cccTATGTGGACTAAACCAGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCAGATTGATCGCcag                                                                                            >  6:23292/1‑148 (MQ=255)
                                                             tATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAgg                                                                                        >  2:367843/1‑149 (MQ=255)
                                                                 tgGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCcct                                                                                    >  5:242258/1‑149 (MQ=255)
                                                                  ggACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCca                                                                                             >  5:352640/1‑139 (MQ=255)
                                                                  ggACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCca                                                                                             <  6:352640/139‑1 (MQ=255)
                                                                   gACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCccttg                                                                                  >  1:153871/1‑149 (MQ=255)
                                                                   gACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCcctt                                                                                   <  8:267103/148‑1 (MQ=255)
                                                                          cccGCCTATACCGACCTGCGCATTGGCTGCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTggg                                                                           <  5:265743/149‑1 (MQ=255)
                                                                           ccGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                          >  1:95171/1‑149 (MQ=255)
                                                                           ccGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                          >  5:62965/1‑149 (MQ=255)
                                                                           ccGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                          >  6:15748/1‑149 (MQ=255)
                                                                           ccGCCTATACCGACCTGCGCATTGGCATCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                          >  1:207359/1‑149 (MQ=255)
                                                                             gCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCccttgcctt                                                                              <  5:250494/143‑1 (MQ=255)
                                                                             gCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCccttgcctt                                                                              >  6:250494/1‑143 (MQ=255)
                                                                                   aCCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCtt                                                                  >  8:262130/1‑149 (MQ=255)
                                                                                   aCCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCtt                                                                  >  2:269233/1‑149 (MQ=255)
                                                                                   aCCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCtt                                                                  >  6:79294/1‑149 (MQ=255)
                                                                                         cTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGTCATCGCCACGGGGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCt                                                            <  3:202342/149‑1 (MQ=255)
                                                                                           gcgcATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTg                                                           >  5:298369/1‑148 (MQ=255)
                                                                                            cgcATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTg                                                           <  1:343019/147‑1 (MQ=255)
                                                                                             gcATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGAc                                                        <  1:269233/149‑1 (MQ=255)
                                                                                             gcATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGAc                                                        <  7:262130/149‑1 (MQ=255)
                                                                                             gcATTGGCTGCGCCGGGACGAGGTAGTTCGCCCCCCGTGGCTGCCGGCCCGGCCATCGCCACGGGGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGCTTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGAc                                                        <  7:81153/149‑1 (MQ=255)
                                                                                                   gCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGa                                                   <  1:146007/148‑1 (MQ=255)
                                                                                                   gCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGa                                                   <  5:79294/148‑1 (MQ=255)
                                                                                                   gCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGa                                                   <  5:273390/148‑1 (MQ=255)
                                                                                                        gAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACg                                             >  8:216095/1‑149 (MQ=255)
                                                                                                         aaGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACg                                             >  6:24012/1‑148 (MQ=255)
                                                                                                              aCGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAgg                                        <  5:15748/148‑1 (MQ=255)
                                                                                                                gATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATGCTCCTGGAGTATGCGCCTGTGAAGCTCATTGATCTCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCTCCTTGACTGGGGGCCACGCAAGGCtt                                     >  8:189779/1‑149 (MQ=255)
                                                                                                                gATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTg                                                           >  3:252269/1‑127 (MQ=255)
                                                                                                                gATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTg                                                           <  4:252269/127‑1 (MQ=255)
                                                                                                                gATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCtt                                     <  2:207359/149‑1 (MQ=255)
                                                                                                                         cGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCGTGTAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGAGTGGGAGCCACGGAAGGCTTGCCACGTg                             <  5:23292/148‑1 (MQ=255)
                                                                                                                             aaCCGTTGATGACGGCACGGCCATCGCCACGGTGGCCGGATCCTCCTGGACTAAGCGCCTGTGAACCAGATTGATCGCCAGCAGGCCCTTGCCTTCGGGCGAGGCTTGCGCCTGGACAGGGAGCCACTGAAGGCTTCCCACGAGCagtt                        >  4:363825/1‑149 (MQ=255)
                                                                                                                              aCCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCtgttg                       >  8:126253/1‑149 (MQ=255)
                                                                                                                                        aCGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCgggg             <  3:93613/149‑1 (MQ=255)
                                                                                                                                        aCGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCggg              <  6:62965/148‑1 (MQ=255)
                                                                                                                                           atccggccACCGTGGCGATGGCCGGATCCTCCTGGAGTATGCGCCTGTGCACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGca          >  5:299129/3‑149 (MQ=255)
                                                                                                                                            cccggccATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGca          >  6:109363/1‑148 (MQ=255)
                                                                                                                                                   aTCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAg  <  3:157409/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAG  >  NC_002947/3070601‑3070896

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: ‑