Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,667,239 T→A 6.6% intergenic (+17/+54) nhaP → / ← PP_4975 NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,667,2390TA 6.6% 180.4 / 2.6 61intergenic (+17/+54)nhaP/PP_4975NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CACTGTCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAC  >  NC_002947/5667092‑5667385
                                                                                                                                                   |                                                                                                                                                  
cACTGTCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCttt                                                                                                                                                   >  6:2006/1‑149 (MQ=255)
cACTGTCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGAGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCttt                                                                                                                                                   >  1:297200/1‑149 (MQ=255)
   tGTCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCg                                                                                                                                                >  6:357878/1‑149 (MQ=255)
     tCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTTGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGACGCTTTGCGGc                                                                                                                                              >  4:281072/1‑149 (MQ=255)
     tCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCGGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGc                                                                                                                                              >  3:150498/1‑149 (MQ=255)
     tCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCgg                                                                                                                                               <  5:357878/148‑1 (MQ=255)
     tCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCg                                                                                                                                                <  2:311710/147‑1 (MQ=255)
       gCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCAAGCCTTGCGAAGGCCGGGGGCGCTTTGGGGcc                                                                                                                                             >  8:326556/1‑148 (MQ=255)
        cTACCGCGGGGCGAGGAACGGGACCTGCGGCTGTCGACCCCCTACATCGGGGTGTTGTCGGCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCAAAGCGGcccc                                                                                                                                           <  8:36950/149‑2 (MQ=255)
           ccGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCAAAGCGGccccg                                                                                                                                          <  1:73687/147‑3 (MQ=255)
              cTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcgcg                                                                                                                                     >  1:78413/1‑149 (MQ=255)
                 ggCGAGGAACGGGCCCTGCTGCTGTCGATCACCTACATCGGGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGGGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGt                                                                                                                                  <  5:245426/149‑1 (MQ=255)
                 ggCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTTGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCggg                                                                                                                                   <  3:281072/148‑1 (MQ=255)
                 ggCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcgc                                                                                                                                      <  3:220794/145‑1 (MQ=255)
                 ggCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcgc                                                                                                                                      >  4:220794/1‑145 (MQ=255)
                 ggCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCggg                                                                                                                                   <  3:318619/148‑1 (MQ=255)
                 ggCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCgg                                                                                                                                    <  2:60296/147‑1 (MQ=255)
                    gAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCATGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCCGGAGaa                                                                                                                               >  3:347306/1‑149 (MQ=255)
                       gAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTGTGCGGCCCATTAGCGGGTgtcccc                                                                                                                            >  6:250447/1‑144 (MQ=255)
                       gAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGAAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCATTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAccc                                                                                                                            >  5:113496/1‑149 (MQ=255)
                         aCGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcgcg                                                                                                                                     >  2:105177/1‑138 (MQ=255)
                         aCGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcgcg                                                                                                                                     <  1:105177/138‑1 (MQ=255)
                         aCGGGACCTGCTGCTGTCGATCACCTACATCGTCGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCcgc                                                                                                                          <  5:250447/149‑1 (MQ=255)
                            ggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCGGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaa                                                                                                                       <  7:326556/149‑1 (MQ=255)
                                 tgctgctgTGGTTGATCTTCTTGGTGGTGTTGTGGTCGATTGTGGTGCAGGGGTTGAGCTTTGGGGGGGTGGTTCGCTAGGTCAGCGACCAGCCTTGAGATTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGc                                                                                                                  <  4:347306/149‑1 (MQ=25)
                                   ctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaatgggc                                                                                                                  <  8:106339/147‑6 (MQ=255)
                                   ctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGcc                                                                                                                 <  5:105200/148‑1 (MQ=255)
                                      ctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGca                                                                                                             <  2:297200/149‑1 (MQ=255)
                                         tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGcac                                                                                                            <  2:78413/147‑1 (MQ=255)
                                         tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                          <  2:185564/149‑1 (MQ=255)
                                           gATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCGGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAgg                                                                                                         <  4:150498/148‑1 (MQ=255)
                                            aTCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGc                                                                                                                                              <  8:53695/110‑1 (MQ=255)
                                            aTCACCTACATCGTGGTGTTGTCGACGATTATGGTGCAGGGGTTGAGAATTGGGCGGGTGGTGCACAAGGTAAGCGCCCAGCATTGAAAATGCCGGGGCCGATTTGCGGc                                                                                                                                              >  7:53695/1‑110 (MQ=255)
                                                       cGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcg                                                                                                                                       <  4:259163/106‑1 (MQ=255)
                                                       cGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTcg                                                                                                                                       >  3:259163/1‑106 (MQ=255)
                                                         tggtgTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTAtt                                                                                          <  5:2006/149‑1 (MQ=255)
                                                              ttgtCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCTAAGAGTCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                     >  1:129973/1‑149 (MQ=255)
                                                                   gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCg                                                                                 >  1:333946/1‑148 (MQ=255)
                                                                   gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCg                                                                                 >  6:218179/1‑148 (MQ=255)
                                                                     cgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGAGCCCAGCCTTGAGAATGCCGTGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTCTGCGCATGCCGaa                                                                               >  8:288590/1‑148 (MQ=255)
                                                                           tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                  >  6:361780/1‑139 (MQ=255)
                                                                           tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                  <  5:361780/139‑1 (MQ=255)
                                                                             gggCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCggg                                                                         <  6:113496/144‑1 (MQ=255)
                                                                                cAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGACCCAGCGAAGAGGCCGGCACTGGCACCCAGGGTATATTCGCCTGCCGAATCCGGGAACTgg                                                                   >  8:256703/1‑149 (MQ=255)
                                                                                 accgcttGACCATTGGGCGGCTGCTGCCCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGCGTCAACCCGCCAAGACGCCGGCACACGCAACCAGCGTCTATACGACTCCCGAATcc                                                                            >  2:85353/6‑139 (MQ=255)
                                                                                 aGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATcc                                                                            <  1:85353/139‑1 (MQ=255)
                                                                                   gggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCgg                                                                                                                                               <  7:160922/70‑1 (MQ=255)
                                                                                   gggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCgg                                                                                                                                               >  8:160922/1‑70 (MQ=255)
                                                                                      ttGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                     <  7:229468/125‑1 (MQ=255)
                                                                                      ttGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                     >  8:229468/1‑125 (MQ=255)
                                                                                        gAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                 >  8:315106/1‑111 (MQ=255)
                                                                                        gAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                 <  7:315106/111‑1 (MQ=255)
                                                                                          gCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACg                                                          <  2:333946/148‑1 (MQ=255)
                                                                                           cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACAAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGt                                                         >  5:24106/1‑148 (MQ=255)
                                                                                             ttGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                          <  1:344041/129‑1 (MQ=255)
                                                                                             ttGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                          >  2:344041/1‑129 (MQ=255)
                                                                                                        gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTg                                           >  5:205860/1‑149 (MQ=255)
                                                                                                                        cccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCg                           <  5:218179/149‑1 (MQ=255)
                                                                                                                                  agaATGCCGGGGCCGCTTTTCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCAAAGGAAACAAGGGTCTATTAGAATGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATc                 >  6:20976/1‑149 (MQ=255)
                                                                                                                                    cgcgaatgggCCGCAAAGCGGCCCATTCGCGGGTGACCCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTCATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATcc                >  3:181002/8‑148 (MQ=255)
                                                                                                                                        ccGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc           >  5:306976/1‑149 (MQ=255)
                                                                                                                                        aatgggCCGCAAAGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc           >  4:17833/4‑149 (MQ=255)
                                                                                                                                                 cTTTGCGGCACATTCGCGGGAGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATTGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc  <  7:288590/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACTGTCGCTACCGCTGGGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAC  >  NC_002947/5667092‑5667385

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