Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 250,267 Δ1 bp 6.2% coding (936/1347 nt) nfeD ← membrane bound peptidase NefD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947250,2670G. 6.2% 206.1 / 7.3 65coding (936/1347 nt)nfeDmembrane bound peptidase NefD
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base . (2/2);  total (29/36)
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.41e-01

GAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAA  >  NC_002947/250120‑250404
                                                                                                                                                   |                                                                                                                                         
gAGTGCCGAGAGTAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCgg                                                                                                                                          <  8:265290/149‑1 (MQ=255)
gAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCgg                                                                                                                                          >  4:300483/1‑149 (MQ=255)
gAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCgg                                                                                                                                          >  6:5008/1‑149 (MQ=255)
gAGTGCCGAGAGCAACGCGAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCgg                                                                                                                                          <  2:67102/149‑1 (MQ=255)
   gtggtcggtggtatctgcctgctgctggcgctgtatgcgctgcagctgttgccggtgagttttgccggggtggcgctgatcctgcttggcatcaccttcatgatCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGa                                                                                                                                       >  5:340901‑M1/105‑149 (MQ=255)
     ggtcggtggtatctgcctgctgctggcgctgtatgcgctgcagctgttgccggtgagttatgccggggtggcgctgatcctgcttggcatcaccttcatgatCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATc                                                                                                                                     >  1:359120‑M1/103‑149 (MQ=255)
     ccGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGAt                                                                                                                                      >  5:180767/1‑148 (MQ=255)
     ccGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGAt                                                                                                                                      >  2:352749/1‑148 (MQ=255)
       gagaGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAAGGATGCCGCCGAAGCCGACCACTCCGAAACTGGGCAAGAAGGCTTCGGCGATctt                                                                                                                                   <  3:47243/149‑3 (MQ=255)
             agcGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGTCCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGt                                                                                                                             <  1:73596/149‑1 (MQ=255)
             agcGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGt                                                                                                                             >  1:210373/1‑149 (MQ=255)
               cgcCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGa                                                                                                                           <  7:89121/149‑1 (MQ=255)
                gcCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGGCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAACCTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGAt                                                                                                                          <  3:240042/149‑1 (MQ=255)
                 ccAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATg                                                                                                                         >  4:370918/1‑149 (MQ=255)
                   aaGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGcc                                                                                                                       <  4:103131/149‑1 (MQ=255)
                         ccGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCa                                                                                                                  <  1:352749/148‑1 (MQ=255)
                         ccGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCa                                                                                                                  <  3:300483/148‑1 (MQ=255)
                           gATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCa                                                                                                                  >  6:321153/1‑146 (MQ=255)
                           gATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCa                                                                                                                  <  5:321153/146‑1 (MQ=255)
                           gATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGa                                                                                                               >  6:377237/1‑149 (MQ=255)
                            aTCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGAt                                                                                                              >  5:287759/1‑149 (MQ=255)
                             tCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAACGCAACGATGCCGCCGAAGCCGCCAACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCAGGAAGGTGATGCCAAGCAGGAt                                                                                                              >  8:266385/1‑148 (MQ=255)
                                   ccAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCTCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAgcgc                                                                                                        >  1:238663/1‑148 (MQ=255)
                                               ctaaaCCGGCAGCCTCGGAGTCGATCAGTAACAATCCGCCTAACACAAAGGCAACGAACCCCCCCAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTACGGCGATCATGAAGGTCAAGCCAAGCAGGATCAGCGCCACCCCGCCaaaa                                                                                           <  3:247744/147‑1 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                                           <  2:210373/149‑1 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCa                                                                                              <  3:282870/146‑1 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCa                                                                                              >  4:282870/1‑146 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                                           <  8:19363/148‑1 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                                           >  7:13888/1‑148 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                                           >  7:19363/1‑148 (MQ=255)
                                               cGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                                           <  8:13888/148‑1 (MQ=255)
                                                gAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAAc                                                                                          <  1:286268/149‑1 (MQ=255)
                                                       ggAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCgg                                                                                   <  3:99112/149‑1 (MQ=255)
                                                       ggAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCgg                                                                                   >  3:244130/1‑149 (MQ=255)
                                                          gCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCaa                                                                                <  5:5008/149‑1 (MQ=255)
                                                                gtgtCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAATCTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTg                                                                          >  6:248511/1‑149 (MQ=255)
                                                                gtgtCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTg                                                                          <  7:106578/149‑1 (MQ=255)
                                                                  gtCGATCATGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCATCCCGGCaaaa                                                                                           <  8:151138/130‑1 (MQ=255)
                                                                  gtCGATCATGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCATCCCGGCaaaa                                                                                           >  7:151138/1‑130 (MQ=255)
                                                                   tCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAg                                                                       >  3:353443/1‑149 (MQ=255)
                                                                   tCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATAATGAAGGTGATTCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCa                                                                        <  6:180767/148‑1 (MQ=255)
                                                                     gATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAgcg                                                                     >  3:113618/1‑149 (MQ=255)
                                                                      aTCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAgcg                                                                     >  4:83860/1‑148 (MQ=255)
                                                                      aTCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAgcg                                                                     >  5:67012/1‑148 (MQ=255)
                                                                         aGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATa                                                                 <  8:140077/149‑1 (MQ=255)
                                                                         aGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATa                                                                 <  5:377237/149‑1 (MQ=255)
                                                                             tCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAgc                                                             <  1:356483/149‑1 (MQ=255)
                                                                                  gcgcCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCcag                                                        >  3:141862/1‑149 (MQ=255)
                                                                                   cgcCGACCACAAAGCCAACGCTGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCcagc                                                       <  8:99025/149‑1 (MQ=255)
                                                                                    gcCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCcagca                                                      >  3:359098/1‑149 (MQ=255)
                                                                                      cGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCcagcagc                                                    <  6:287759/149‑1 (MQ=255)
                                                                                        aCCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCcagcagcag                                                  <  1:265865/149‑1 (MQ=255)
                                                                                         ccACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcagg                                                 <  5:248511/149‑1 (MQ=255)
                                                                                           acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcagg                                                 >  3:123837/1‑147 (MQ=255)
                                                                                              aaGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGTCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGAt                                            <  2:238663/149‑1 (MQ=255)
                                                                                              aaGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGAt                                            <  6:224249/149‑1 (MQ=255)
                                                                                                ggCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATac                                          >  3:1344/1‑149 (MQ=255)
                                                                                                  cAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATacca                                        <  1:49402/149‑1 (MQ=255)
                                                                                                    aCGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATaccacc                                      <  5:269353/149‑1 (MQ=255)
                                                                                                     cGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGAGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCg                                     <  5:814/149‑1 (MQ=255)
                                                                                                           CGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGAccacc                               <  7:139623‑M1/149‑2 (MQ=255)
                                                                                                           CGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACcaca                               <  7:293330‑M1/149‑1 (MQ=255)
                                                                                                                             cGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCt             >  1:95138/1‑149 (MQ=255)
                                                                                                                                aaCTGGGCAAGAAGGCTTCGGCGATCATGAAGGGGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTggg          <  8:293330/149‑1 (MQ=255)
                                                                                                                                     ggCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCAt     >  6:83575/1‑149 (MQ=255)
                                                                                                                                      gCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATg    <  4:141862/149‑1 (MQ=255)
                                                                                                                                      gCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATg    <  4:1344/149‑1 (MQ=255)
                                                                                                                                        aagaagGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGaa  >  5:54268/1‑149 (MQ=255)
                                                                                                                                        aagaagGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCAGGCAATACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCCGCAGGCAGATACCACCGACCACACCGCCCACCGCAGCACCTGGGTTCAGGaa  >  6:172813/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAA  >  NC_002947/250120‑250404

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