Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,070,757 G→A 8.0% intergenic (+325/‑29) 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,7570GA 8.0% 178.0 / 3.3 62intergenic (+325/‑29)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (3/2);  total (33/29)
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.28e-01

TTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAGCTCAATGC  >  NC_002947/3070616‑3070904
                                                                                                                                             |                                                                                                                                                   
ttAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                              <  1:198232/149‑1 (MQ=255)
 tAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                              >  2:66546/1‑148 (MQ=255)
 tAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCACCGTGGCGATGGCCgg                                                                                                                                              <  1:66546/148‑1 (MQ=255)
       ttCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGAtcctcc                                                                                                                                       >  6:337757/1‑149 (MQ=255)
       ttCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGAtcctcc                                                                                                                                       >  7:267103/1‑149 (MQ=255)
           aCCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTg                                                                                                                                     >  2:343019/1‑147 (MQ=255)
           aCCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTg                                                                                                                                     >  7:357505/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                                                                                                               <  5:337757/149‑1 (MQ=255)
                               gaagaGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCtga                                                                                                               <  8:362940/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                                                                                  >  6:15748/1‑149 (MQ=255)
                                                            ccGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                                  >  1:95171/1‑149 (MQ=255)
                                                            ccGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTgggg                                                                                  >  5:62965/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:273390/148‑1 (MQ=255)
                                                                                    gCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGa                                                           <  5:79294/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)
                                                                                                                                            ggTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAGCTCAATGc  <  7:126253/149‑1 (MQ=255)
                                                                                                                                            ggTGGCCGGATCCTCCTGGAGTATGCGACTGTGAACCTGATTGATCTCCAGGAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGGTGTTGTACGCCGGGGGCATCATCGAGCTCAATGc  <  7:189779/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TTAATGCTTCTACCGGAATCCTCCGGCACAAGAAGAGGACTTGCCCTATGTGGACTAAACCCGCCTATACCGACCTGCGCATTGGCTTCGAAGTGACGATGTATTTCGCCAACCGTTGATGACGGCCCGGCCATCGCCACGGTGGCCGGATCCTCCTGGAGTATGCGCCTGTGAACCTGATTGATCGCCAGCAGGCCCTTGCCTTGGGGCGAGGCTTGCGCCTGGACTGGGAGCCACGGAAGGCTTGCCACGTGCTGTTGTACGCCGGGGGCATCATCGAGCTCAATGC  >  NC_002947/3070616‑3070904

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