Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,294,114 C→T 8.2% 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,1140CT 8.2% 176.6 / 4.3 61intergenic (+15/+25)pdxH/PP_1130pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

ATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCT  >  NC_002947/1293972‑1294262
                                                                                                                                              |                                                                                                                                                    
aTTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCgg                                                                                                                                                >  2:173612/1‑149 (MQ=255)
   gggggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGcc                                                                                                                                              <  4:97968/148‑1 (MQ=255)
       gggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgt                                                                                                                                         <  1:239918/149‑2 (MQ=255)
         gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATTGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCa                                                                                                                                       >  1:28000/1‑149 (MQ=255)
         gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAg                                                                                                                                          <  6:222942/146‑1 (MQ=255)
         gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAg                                                                                                                                          >  5:222942/1‑146 (MQ=255)
                   gttgCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGc                                                                                                                             >  1:32758/1‑149 (MQ=255)
                       cccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATTGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcc                                                                                                                         <  2:28000/149‑1 (MQ=255)
                       cccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCa                                                                                                                                <  1:281265/142‑1 (MQ=255)
                       cccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCa                                                                                                                                >  2:281265/1‑142 (MQ=255)
                       cccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcc                                                                                                                         >  6:195876/1‑149 (MQ=255)
                         cGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATAGTTTGAACTATCTGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcccc                                                                                                                       <  2:167623/149‑1 (MQ=255)
                                        cTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGccc                                                                                                                                             <  4:35076/112‑1 (MQ=255)
                                        cTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGccc                                                                                                                                             >  3:35076/1‑112 (MQ=255)
                                           gCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcg                                                                                                     <  3:173001/149‑1 (MQ=255)
                                           gCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTc                                                                                                      <  5:195876/148‑1 (MQ=255)
                                            cAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCAGATAACCTGCCGCCTcgc                                                                                                    <  1:120443/149‑1 (MQ=255)
                                              ggggCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCgcgc                                                                                                  <  2:32758/149‑1 (MQ=255)
                                                ggcgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCgg                                                                                                                     <  1:165342/128‑1 (MQ=255)
                                                ggcgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCgg                                                                                                                     >  2:165342/1‑128 (MQ=255)
                                                                   tGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCt                                                                                                 <  4:286453/129‑1 (MQ=255)
                                                                   tGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCt                                                                                                 >  3:286453/1‑129 (MQ=255)
                                                                   tGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGc                                                                             <  1:129757/149‑1 (MQ=255)
                                                                    gATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                   >  4:293856/1‑126 (MQ=255)
                                                                    gATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgcg                                                                                                   <  3:293856/126‑1 (MQ=255)
                                                                     aTCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                            <  8:100420/116‑1 (MQ=255)
                                                                     aTCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                            >  7:100420/1‑116 (MQ=255)
                                                                         tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGc                                                                       >  1:219904/1‑149 (MQ=255)
                                                                            gAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGccc                                                                     >  1:60268/1‑148 (MQ=255)
                                                                                tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCt                                                                >  4:278600/1‑149 (MQ=255)
                                                                                tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCACCAACCACCCCGCCAGTTGCTTACGCCCCACCt                                                                >  3:190003/1‑149 (MQ=255)
                                                                                    ggCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGa                                                                                                                                   <  5:100285/78‑1 (MQ=255)
                                                                                    ggCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGa                                                                                                                                   >  6:100285/1‑78 (MQ=255)
                                                                                    ggCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtctc                                                            >  6:138064/1‑149 (MQ=255)
                                                                                     gCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           >  3:57669/1‑149 (MQ=255)
                                                                                     gCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           >  6:161849/1‑149 (MQ=255)
                                                                                       tggtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGcccc                                                                    >  3:35912/1‑138 (MQ=255)
                                                                                       tggtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGcccc                                                                    <  4:35912/138‑1 (MQ=255)
                                                                                         gtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCACCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGcc                                                       >  7:49579/1‑149 (MQ=255)
                                                                                                       gCTGCGAGAGCGGCTGGTGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAc                                         >  3:184546/1‑149 (MQ=255)
                                                                                                       gCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAc                                         >  1:232040/1‑149 (MQ=255)
                                                                                                       gCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAc                                         >  6:190568/1‑149 (MQ=255)
                                                                                                         tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                               <  7:279472/109‑1 (MQ=255)
                                                                                                         tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                               >  8:279472/1‑109 (MQ=255)
                                                                                                         tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCt                                       >  4:198241/1‑149 (MQ=255)
                                                                                                           cGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCttt                                     >  2:206333/1‑149 (MQ=255)
                                                                                                            gagaGCGGCGGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACACCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTc                                    <  7:74395/149‑1 (MQ=255)
                                                                                                                gCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGc                                <  2:232040/149‑1 (MQ=255)
                                                                                                                  ggCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCAt                              <  3:198241/149‑1 (MQ=255)
                                                                                                                    cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGa                            <  1:206333/149‑1 (MQ=255)
                                                                                                                    cTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGa                            <  5:138064/149‑1 (MQ=255)
                                                                                                                       gcgcCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACACCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcg                         >  6:249819/1‑149 (MQ=255)
                                                                                                                           ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgcg                     >  1:100416/1‑149 (MQ=255)
                                                                                                                                 cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTc               <  4:184546/149‑1 (MQ=255)
                                                                                                                                  aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa              <  1:273633/149‑1 (MQ=255)
                                                                                                                                  aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa              <  2:219904/149‑1 (MQ=255)
                                                                                                                                   aCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgc                        >  7:250814/1‑138 (MQ=255)
                                                                                                                                   aCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgc                        <  8:250814/138‑1 (MQ=255)
                                                                                                                                       gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGAt          <  2:60268/148‑1 (MQ=255)
                                                                                                                                       gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATc         <  2:100416/149‑1 (MQ=255)
                                                                                                                                        ggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCg        >  4:207197/1‑149 (MQ=255)
                                                                                                                                              cAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt  >  2:279065/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCT  >  NC_002947/1293972‑1294262

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