Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,294,114 C→T 10.1% intergenic (+15/+25) pdxH → / ← PP_1130 pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,294,1140CT10.1% 191.7 / 9.5 69intergenic (+15/+25)pdxH/PP_1130pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base T (3/4);  total (31/38)
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.11e-01

CTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCC  >  NC_002947/1293966‑1294260
                                                                                                                                                    |                                                                                                                                                   
cTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGgccgcc                                                                                                                                                     <  8:127332/149‑1 (MQ=255)
  gAGCATTGTGGGGGAGACCGTGTGGTGCCCGCGCGCCTCGAGTTCTGGCAGGGGCTGGCAAGTCGTCTGCATGCTCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCaa                                                                                                                                                  <  8:233660/149‑1 (MQ=255)
    gCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                >  2:58665/1‑149 (MQ=255)
      aTTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCgg                                                                                                                                              >  7:59729/1‑149 (MQ=255)
           gggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCa                                                                                                                                         <  8:317131/149‑1 (MQ=255)
           gggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCa                                                                                                                                         <  3:217264/149‑1 (MQ=255)
              ggATACCGTTTGTTGCCCGAGCGCATCTAGTTCTGGCAGGGGCGTGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTTGCGGCCCCAgtt                                                                                                                                      <  8:92579/149‑3 (MQ=255)
                 tACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAtt                                                                                                                                   <  3:103055/149‑1 (MQ=255)
                     gTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGGCCCAGCCATTGATc                                                                                                                               <  4:20615/149‑1 (MQ=255)
                     gTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATc                                                                                                                               <  8:59729/149‑1 (MQ=255)
                        tgttgCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                             <  2:43744/148‑1 (MQ=255)
                              ccGAGCGGATCGAGGTCTGTCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAccc                                                                                                                      <  8:221568/149‑1 (MQ=255)
                              ccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAccc                                                                                                                      >  7:353273/1‑149 (MQ=255)
                                     cATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATaa                                                                                                               <  1:58665/149‑1 (MQ=255)
                                      aTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                             <  5:318903/134‑1 (MQ=255)
                                      aTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                             >  6:318903/1‑134 (MQ=255)
                                       tCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                             >  5:208840/1‑149 (MQ=255)
                                       tCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                             >  3:343424/1‑149 (MQ=255)
                                          aGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                          >  6:132534/1‑149 (MQ=255)
                                              cTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgccgc                                                                                                       <  2:334894/148‑1 (MQ=255)
                                               tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAc                                                                                                              <  6:91330/140‑1 (MQ=255)
                                               tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAc                                                                                                              >  5:91330/1‑140 (MQ=255)
                                               tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCt                                                                                                     >  6:395094/1‑149 (MQ=255)
                                                 gCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAACCGGCCCCAGACATTGATCAAGCTCCACCGGATAACCTgccgcagag                                                                                                   >  8:167306/1‑145 (MQ=255)
                                                 gCAGGGGCGGGCAAGTCGGCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTc                                                                                                    <  8:54716/148‑1 (MQ=255)
                                                    ggggCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCgg                                                                                                                                              <  6:88549/103‑1 (MQ=255)
                                                                 cgtcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAcccc                                                                                    <  7:199279/148‑1 (MQ=255)
                                                                  gtcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGc                                                                                  <  8:16160/149‑1 (MQ=255)
                                                                     tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                 >  7:169531/1‑131 (MQ=255)
                                                                     tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                 <  8:169531/131‑1 (MQ=255)
                                                                     tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCt                                                                                               <  8:81815/133‑1 (MQ=255)
                                                                     tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCt                                                                                               >  7:81815/1‑133 (MQ=255)
                                                                       cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCa                                                                                <  5:28715/146‑1 (MQ=255)
                                                                          gATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCt                                                                          >  6:343768/1‑149 (MQ=255)
                                                                                   aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                       <  6:310397/127‑1 (MQ=255)
                                                                                   aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                       >  5:310397/1‑127 (MQ=255)
                                                                                     cTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAcc                                                               >  4:67990/1‑149 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAcccc                                                                                    >  3:309123/1‑127 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAcccc                                                                                    <  4:309123/127‑1 (MQ=255)
                                                                                         cGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                 >  3:317508/1‑111 (MQ=255)
                                                                                         cGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                 <  4:317508/111‑1 (MQ=255)
                                                                                          ggCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtctc                                                          <  6:208840/149‑1 (MQ=255)
                                                                                               gtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTGGCGCTCCACCCACCCCGCCAGGTGCTTACGCCCCCCCTTCTCCTGcc                                                     >  4:59689/1‑149 (MQ=255)
                                                                                                  gACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCcgc                                                   <  8:353273/148‑1 (MQ=255)
                                                                                                   aCGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCTAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCCCCTTCTCCTGCcgcg                                                  >  1:308402/1‑149 (MQ=255)
                                                                                                    cGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGcc                                                <  3:67990/149‑1 (MQ=255)
                                                                                                      ggCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTg                                              <  5:395094/149‑1 (MQ=255)
                                                                                                       gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgc                                             >  4:187515/1‑149 (MQ=255)
                                                                                                       gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgc                                             >  2:146153/1‑149 (MQ=255)
                                                                                                         aGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgcgc                                           <  4:343424/149‑1 (MQ=255)
                                                                                                                gCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCCCCTGCTCCTGCCGCGCCTGCGCGATCCtt                                    >  5:48750/1‑149 (MQ=255)
                                                                                                                 cGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGc                                                                                  >  2:261689/1‑102 (MQ=255)
                                                                                                                 cGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGc                                                                                  <  1:261689/102‑1 (MQ=255)
                                                                                                                  gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTc                                  <  2:30496/149‑1 (MQ=255)
                                                                                                                     agCGGCTGGCGCGCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGACTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCa                                <  2:360921/148‑1 (MQ=255)
                                                                                                                      gCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAg                               >  8:139441/1‑148 (MQ=255)
                                                                                                                       cGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCa                             <  5:343768/149‑1 (MQ=255)
                                                                                                                                cccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCCCCTTCTCCTGCCGCGCCTGCGCCCACCGTTCCAGCATGAAcgcgc                    >  3:78551/1‑149 (MQ=255)
                                                                                                                                cccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgc                    <  2:308402/149‑1 (MQ=255)
                                                                                                                                   tcaatggcTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGtt                 >  3:323682/8‑149 (MQ=255)
                                                                                                                                    gAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAATGCGCGttt                >  4:85029/1‑149 (MQ=255)
                                                                                                                                       cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGttttt              >  1:25520/1‑148 (MQ=255)
                                                                                                                                        aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa            <  5:132534/149‑1 (MQ=255)
                                                                                                                                         gcTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa            <  4:323682/147‑1 (MQ=255)
                                                                                                                                          cTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTAc                                                                       <  3:190730/88‑1 (MQ=39)
                                                                                                                                          cTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTAc                                                                       >  4:190730/1‑88 (MQ=39)
                                                                                                                                          cTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg           >  4:31192/1‑148 (MQ=255)
                                                                                                                                          cTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCt                                                                          >  1:175553/1‑85 (MQ=255)
                                                                                                                                          cTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCt                                                                          <  2:175553/85‑1 (MQ=255)
                                                                                                                                              ggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCg      <  7:139441/149‑1 (MQ=255)
                                                                                                                                                  gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGccc  >  5:286577/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCC  >  NC_002947/1293966‑1294260

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