Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,784 C→G 5.9% E794Q (GAG→CAG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,7840CG 5.9% 218.8 / 2.8 68E794Q (GAG→CAG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/38);  new base G (2/2);  total (28/40)
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.70e-01

AACCCCACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGG  >  NC_002947/1725636‑1725930
                                                                                                                                                    |                                                                                                                                                  
aaCCCCACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTc                                                                                                                                                    >  2:17967/1‑149 (MQ=255)
    ccACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCg                                                                                                                                                >  8:275319/1‑149 (MQ=255)
          aCTTTAGCGGATTTTGCCGACCGATGAAGCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTg                                                                                                                                          <  5:212247/149‑1 (MQ=255)
          aCTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTg                                                                                                                                          <  6:309914/149‑1 (MQ=255)
          aCTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTAGCCGAGCTTg                                                                                                                                          >  5:205234/1‑149 (MQ=255)
            tttAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTc                                                                                                                                        <  4:978/149‑1 (MQ=255)
              tAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCGAGCTGGCCGAgggcgtg                                                                                                                                        <  7:114286/147‑7 (MQ=255)
              tAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTc                                                                                                                                        <  7:275319/147‑1 (MQ=255)
              tAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTc                                                                                                                                        <  6:139915/147‑1 (MQ=255)
              tAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCg                                                                                                                                       <  5:179400/148‑1 (MQ=255)
               aGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGt                                                                                                                                      <  1:253747/148‑1 (MQ=255)
                gCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGa                                                                                                                                    <  6:297168/149‑1 (MQ=255)
                 cGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGAt                                                                                                                                   <  7:168442/149‑1 (MQ=255)
                    aTTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAAc                                                                                                                                <  1:17967/149‑1 (MQ=255)
                          ccgaccgaTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGACCAGCTCGCCGAGCTTGTCGTGATAACTCAGct                                                                                                                          <  8:241106/149‑1 (MQ=255)
                                 aTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgcc                                                                                                                      <  2:133702/146‑1 (MQ=255)
                                 aTGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgcc                                                                                                                      >  1:133702/1‑146 (MQ=255)
                                    aaCCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTgcgc                                                                                                                >  6:308601/1‑149 (MQ=255)
                                       cGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAAc                                                                                                                                <  4:82129/130‑1 (MQ=255)
                                       cGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAAc                                                                                                                                >  3:82129/1‑130 (MQ=255)
                                          tcgagGGTCGAAGTTTCCAGCCCGTCGTCCCGGTGCAGGGCCACCTGGCGGATCGACAGGCTCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGc                                                                                                          <  6:334958/147‑1 (MQ=255)
                                             agGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTg                                                                                                       >  2:291230/1‑149 (MQ=255)
                                                 tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGc                                                                                                                           >  3:194524/1‑125 (MQ=255)
                                                 tCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGc                                                                                                                           <  4:194524/125‑1 (MQ=255)
                                                   gAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCg                                                                                                         <  8:42024/141‑1 (MQ=255)
                                                   gAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCg                                                                                                         >  7:42024/1‑141 (MQ=255)
                                                               ccGTCGTCACGGGGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAg                                                                                     <  3:31495/149‑1 (MQ=255)
                                                                  tcgtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCttcagc                                                                                    >  2:238146/1‑142 (MQ=255)
                                                                      cACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTgcgc                                                                                                                <  5:116392/115‑1 (MQ=255)
                                                                      cACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTgcgc                                                                                                                >  6:116392/1‑115 (MQ=255)
                                                                       aCGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTc                                                                              <  1:106547/148‑1 (MQ=255)
                                                                       aCGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTACGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggcg                                                                              <  1:238146/148‑3 (MQ=255)
                                                                           tGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaat                                                                         >  3:27500/1‑149 (MQ=255)
                                                                               gggCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAg                                                                     >  1:16072/1‑149 (MQ=255)
                                                                                     ccTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGt                                                                               <  2:315923/133‑1 (MQ=255)
                                                                                     ccTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGt                                                                               >  1:315923/1‑133 (MQ=255)
                                                                                      cTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTCCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCgg                                                               >  3:19676/1‑148 (MQ=255)
                                                                                       tGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGaa                                                             >  1:102034/1‑149 (MQ=255)
                                                                                         gCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                      <  8:204672/122‑1 (MQ=255)
                                                                                           ggATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCTCGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         <  5:308601/149‑1 (MQ=255)
                                                                                           ggATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         >  3:80009/1‑149 (MQ=255)
                                                                                             aTCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGt                                                       >  7:262446/1‑149 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTTGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCg                                                                <  3:214772/136‑1 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgc                                                                                                                       >  7:161037/1‑81 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGcttgc                                                                                                                       <  8:161037/81‑1 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                 <  7:269946/135‑1 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGc                                                                 >  8:269946/1‑135 (MQ=255)
                                                                                                 aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTAGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCg                                                                >  4:214772/1‑136 (MQ=255)
                                                                                                     gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCga                                               <  4:27500/149‑1 (MQ=255)
                                                                                                      cgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                              <  3:161657/149‑1 (MQ=255)
                                                                                                      cgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaa                                              <  1:291230/149‑1 (MQ=255)
                                                                                                         aGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaagaa                                           <  2:16072/149‑1 (MQ=255)
                                                                                                              ggCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTc                                        <  4:80009/147‑1 (MQ=255)
                                                                                                                  gACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         >  2:310621/1‑126 (MQ=255)
                                                                                                                  gACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         <  1:310621/126‑1 (MQ=255)
                                                                                                                  gACAGTACTCACTTGGCCACGCCCTCGGACAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         <  7:161183/126‑1 (MQ=255)
                                                                                                                  gACAGTACTCACTTGGCCACGCCCTCGGACAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGc                                                         >  8:161183/1‑126 (MQ=255)
                                                                                                                       tACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGtt                               >  3:39116/1‑147 (MQ=255)
                                                                                                                        aCTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCTGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                            >  3:333856/1‑149 (MQ=255)
                                                                                                                         cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCTGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                            <  4:333856/148‑1 (MQ=255)
                                                                                                                         cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGc                                  <  7:280357/142‑1 (MQ=255)
                                                                                                                         cTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGc                                  >  8:280357/1‑142 (MQ=255)
                                                                                                                            gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGGGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCa                        >  5:303644/1‑149 (MQ=255)
                                                                                                                                   gacaagcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGa                 >  6:223761/7‑149 (MQ=255)
                                                                                                                                   cACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGc                                  <  2:11734/132‑1 (MQ=255)
                                                                                                                                       cccTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGa             <  6:303644/149‑1 (MQ=255)
                                                                                                                                        gcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                            <  2:62486/132‑1 (MQ=255)
                                                                                                                                        gcTCGGCGAGCTGGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAg                            >  1:62486/2‑133 (MQ=255)
                                                                                                                                                  cTCGCCGAGGTTTGCGTGATAACTCACCTTGCCGTGCGGGTCGGCTTGCTCCAAGCCTTGCTTTAGCAAAGGCTGTAAGACTTTGCGGCACCTCGTTTTGTGGCAGTACTCCGGGTCCTTCACGCCATGCCGGATCGACAGGCGCTggg  <  6:205234/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACCCCACGGACTTTAGCGGATTTTGCCGACCGATGAACCGTTAGAGGGTCGAAGTTTCCAGCCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGG  >  NC_002947/1725636‑1725930

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