Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,311 A→T 6.9% L385Q (CTG→CAG)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,3110AT 6.9% 176.3 / 2.8 58L385Q (CTG→CAG) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base T (2/2);  total (31/27)
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.13e-01

ATAGGTGAAGGTACATGCCAGCGCAGGTACAGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTCGGCTGTG  >  NC_002947/461167‑461459
                                                                                                                                                |                                                                                                                                                    
ttAGGTGAATGTACATTCGAGAGGAGGTACAGCGCTTCGAATTGTCGAGCGTGAGGAAGCTGTTCACTGATCTTCAGCGCTTCCGGCAGAGGTTGCGCCTTGTCATGTGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCt                                                                                                                                                   <  4:230851/147‑1 (MQ=255)
     tGAAGGTACATGCCAGCGCAGTTACAGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCTGCTGCCagt                                                                                                                                             <  8:119942/149‑3 (MQ=255)
               tGCCAGCGCAGGTACAGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTc                                                                                                                                   >  7:118893/1‑149 (MQ=255)
                        aGGTACAGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTgg                                                                                                                          <  8:271827/149‑1 (MQ=255)
                              aGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCt                                                                                                                            >  4:255517/1‑141 (MQ=255)
                              aGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCt                                                                                                                            <  3:255517/141‑1 (MQ=255)
                                            gcgcgcGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGgcgc                                                                                                      >  1:273728/1‑149 (MQ=255)
                                            gcgcgcGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGgcgc                                                                                                      >  5:212014/1‑149 (MQ=255)
                                               cgcgGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCaa                                                                                                   <  1:26891/149‑1 (MQ=255)
                                                     cggAAGCTGGGCCCTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCt                                                                                             <  8:164587/149‑1 (MQ=255)
                                                      ggAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTg                                                                                            <  6:212014/149‑1 (MQ=255)
                                                              gTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCt                                                                                     <  6:312485/148‑1 (MQ=255)
                                                               tCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCACTCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTg                                                                                   <  2:200644/149‑1 (MQ=255)
                                                                 aCTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTggg                                                                                 >  7:76764/1‑149 (MQ=255)
                                                                 aCTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTggg                                                                                 >  5:302164/1‑149 (MQ=255)
                                                                 aCTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  <  8:118893/148‑1 (MQ=255)
                                                                 aCTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  <  6:275226/148‑1 (MQ=255)
                                                                   tGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  <  8:185684/146‑1 (MQ=255)
                                                                   tGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  >  7:185684/1‑146 (MQ=255)
                                                                   tGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTa                                                                               >  1:301219/1‑149 (MQ=255)
                                                                       cTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATc                                                                             >  2:253163/1‑147 (MQ=255)
                                                                        ttCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAg                                                                          >  7:279171/1‑149 (MQ=255)
                                                                         tCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGATTGGCAGCTGCTGCCAGTCgac                                                                                                                                         >  1:208031/1‑83 (MQ=38)
                                                                         tCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGATTGGCAGCTGCTGCCAGTCgac                                                                                                                                         <  2:208031/85‑3 (MQ=38)
                                                                          cAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAg                                                                          <  6:302164/147‑1 (MQ=255)
                                                                            gcgcTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAg                                                                          >  5:124449/1‑145 (MQ=255)
                                                                            gcgcTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAg                                                                          <  6:124449/145‑1 (MQ=255)
                                                                             cgcTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTgg                                                                     >  4:149300/1‑149 (MQ=255)
                                                                              gcTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGc                                                                    >  1:210382/1‑149 (MQ=255)
                                                                                  ccGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCCGCTTGGCGagc                                                                >  4:190455/1‑147 (MQ=255)
                                                                                     gCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  >  7:220972/1‑128 (MQ=255)
                                                                                     gCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTgg                                                                                  <  8:220972/128‑1 (MQ=255)
                                                                                     gCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTg                                                                                   <  5:242703/127‑1 (MQ=255)
                                                                                     gCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTg                                                                                   >  6:242703/1‑127 (MQ=255)
                                                                                         cGGTTGCGCCTTGGCATGGGTGAGGTCAGCCTTCAGCAATGGGGCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGc                                                         <  8:126765/149‑1 (MQ=255)
                                                                                           gTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCa                                                       <  8:279171/149‑1 (MQ=255)
                                                                                                gcCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCt                                                  >  7:90802/1‑149 (MQ=255)
                                                                                                gcCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCt                                                  >  4:144218/1‑149 (MQ=255)
                                                                                                gcCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCt                                                  <  2:301219/149‑1 (MQ=255)
                                                                                                gcCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCt                                                  >  6:42750/1‑149 (MQ=255)
                                                                                                        aTGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcggc                                          >  1:310882/1‑149 (MQ=255)
                                                                                                           gggCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcggcgg                                        <  2:210382/148‑1 (MQ=255)
                                                                                                           gggCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcggcgg                                        <  2:250662/148‑1 (MQ=255)
                                                                                                                ggTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAAc                                  <  8:90802/149‑1 (MQ=255)
                                                                                                                 gTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcggcgg                                        <  1:235588/142‑1 (MQ=255)
                                                                                                                 gTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcggcgg                                        >  2:235588/1‑142 (MQ=255)
                                                                                                                 gTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCaagt                                 >  5:278384/1‑147 (MQ=255)
                                                                                                                  tCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCt                                                                                             >  5:151657/1‑88 (MQ=255)
                                                                                                                  tCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCt                                                                                             <  6:151657/88‑1 (MQ=255)
                                                                                                                    aGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcgg                                           <  4:167687/136‑1 (MQ=255)
                                                                                                                    aGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCcgg                                           >  3:167687/1‑136 (MQ=255)
                                                                                                                     gCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAg                             >  3:40420/1‑149 (MQ=255)
                                                                                                                      ccTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCGATGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCCACGCGAgc                            >  6:13913/1‑149 (MQ=255)
                                                                                                                                   ggcGATTGGCAGCTGCTGCCAGTCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACtgt               >  6:204749/3‑149 (MQ=255)
                                                                                                                                   gTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACtgt               >  2:141401/1‑149 (MQ=255)
                                                                                                                                   gTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACtgt               >  8:197745/1‑149 (MQ=255)
                                                                                                                                    tCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTgg                     >  7:42103/1‑142 (MQ=255)
                                                                                                                                    tCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTgg                     <  8:42103/142‑1 (MQ=255)
                                                                                                                                           gcagcaGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTc         <  1:141401/147‑1 (MQ=255)
                                                                                                                                           gcagcaGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTCgg       >  6:304129/1‑149 (MQ=255)
                                                                                                                                               cagcTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTCGGCtgt   <  7:197745/149‑1 (MQ=255)
                                                                                                                                                agcTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTCGGCtgtg  <  2:310882/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATAGGTGAAGGTACATGCCAGCGCAGGTACAGCGCTGCGAATTGGCGCGCGGGCGGAAGCTGGTCACTGATCTTCAGCGCTGCCGGCAGCGGTTGCGCCTTGGCATGGGGCAGGTCAGCCTTCAGCAATGGGTCGACTGGCAGCAGCTGCCAATCGCCATGGTCGGCACCTGGCAGGCTTTCCAGCAGGGCGCCAAGCGCCTGCAACCCTTGGGTATCCAGCTTGGCGAACGGCTGGCCAGTGCTGTCCCGGCGGGCCAACTCGAGCGCGCTGGCACTGTGCGCTCGGCTGTG  >  NC_002947/461167‑461459

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: