Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,886,769 C→G 6.4% A83P (GCA→CCA)  PP_5161 ← lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,886,7690CG 6.4% 249.3 / 6.2 78A83P (GCA→CCA) PP_5161lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base G (3/2);  total (44/34)
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.39e-01

CTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCGG  >  NC_002947/5886621‑5886917
                                                                                                                                                      |                                                                                                                                                    
cTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt‑‑gc                                                                                                                                                      <  2:342451/149‑1 (MQ=255)
      aGGCCGTGGGCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGTAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCAT‑‑GGCTGccc                                                                                                                                                <  3:252592/149‑1 (MQ=255)
        gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGa                                                                                                                                              >  6:391250/1‑149 (MQ=255)
          cGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGTCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGAcc                                                                                                                                            <  4:416457/149‑1 (MQ=255)
          cGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGAc                                                                                                                                             >  4:212913/1‑148 (MQ=255)
           gTGGCCGGGCAGGTCACAGCTGATCACCGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGAc                                                                                                                                             >  7:100086/1‑147 (MQ=255)
              gCCGGGCAGGTCACAGCTGATCACTGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAgg                                                                                                                                         <  8:100086/148‑1 (MQ=255)
                      ggTCACAGGTGATCACGGCGAATCCCTGCTTGAGTGCCCACTCGATCACATGTCTGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAAcccc                                                                                                                                <  3:66592/149‑1 (MQ=255)
                          acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTg                                                                                                                            >  5:292546/1‑149 (MQ=255)
                           caGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt‑‑gcctgc                                                                                                                                                  <  4:231094/126‑1 (MQ=255)
                           caGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt‑‑gcctgc                                                                                                                                                  >  3:231094/1‑126 (MQ=255)
                           caGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                           <  7:332460/149‑1 (MQ=255)
                           caGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                           <  7:180901/149‑1 (MQ=255)
                           caGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATCTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                           >  3:273264/1‑149 (MQ=255)
                              cTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                         >  8:281593/1‑148 (MQ=255)
                              cTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTg                                                                                                                        <  6:15653/149‑1 (MQ=255)
                               tGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCACGCTGCATg                                                                                                                        >  4:102996/1‑148 (MQ=255)
                                 aTCACGGTGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                     <  3:274310/149‑1 (MQ=255)
                                  tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAg                                                                                                                    <  6:292546/149‑1 (MQ=255)
                                         gAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACcca                                                                                                             <  4:3095/149‑1 (MQ=255)
                                              ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccagccag                                                                                                        >  6:143600/1‑149 (MQ=255)
                                                tGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt‑‑gc                                                                                                                                                      >  4:379997/1‑101 (MQ=255)
                                                tGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt‑‑gc                                                                                                                                                      <  3:379997/101‑1 (MQ=255)
                                                  cTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                    <  5:283964/149‑1 (MQ=255)
                                                     gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                  <  7:205687/148‑1 (MQ=255)
                                                     gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                  <  4:371900/148‑1 (MQ=255)
                                                      aGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGc                                                                                                <  1:198141/149‑1 (MQ=255)
                                                           ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCgg                                                                                             >  4:82407/1‑147 (MQ=255)
                                                           ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCa                                                                                           <  7:397449/149‑1 (MQ=255)
                                                           ccACTCGATCACATGTAGGTAGAGGACCATGTGGTCCTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCAAA‑‑GCCTGCCAGACCAGGACGAACCCCGCTGCCTCAAAGCAACCCACCCACCAGAGTCCCgg                                                                                             >  6:62354/1‑147 (MQ=255)
                                                                  aTCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                                    >  4:399430/1‑149 (MQ=255)
                                                                     acaTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGcc                                                                                 <  5:143600/149‑1 (MQ=255)
                                                                      caTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAAGCCa                                                                                <  5:62354/149‑1 (MQ=255)
                                                                       ctggcGGTACAGGCCCATGTGGGCGTCATAGCCGTGCAGCAGGCACAGCGTCGCGCGTGGCTGCTCGTGCAGCGAT‑‑GCCTGCCAGACCAGGTCGATCGCCGCTTCCTCAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCat                                                                                    <  2:330748/140‑2 (MQ=255)
                                                                       aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCAAACCAGCCAGCTGTGTGCCGCCAGGTTCAGGCCa                                                                                >  7:49353/1‑148 (MQ=255)
                                                                            ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                           <  4:3106/148‑1 (MQ=255)
                                                                              tAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcg                                                                        >  7:398409/1‑149 (MQ=255)
                                                                                   gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGt                                                                   <  3:399430/149‑1 (MQ=255)
                                                                                           ggTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCAGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgc                                                           >  1:96261/1‑149 (MQ=255)
                                                                                            gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAg                                                                 <  2:413228/142‑1 (MQ=255)
                                                                                            gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAg                                                                 >  1:413228/1‑142 (MQ=255)
                                                                                             tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCa                                                                                     <  4:140407/121‑1 (MQ=255)
                                                                                             tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCa                                                                                     >  3:140407/1‑121 (MQ=255)
                                                                                                tAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTc                                                               <  5:34202/140‑1 (MQ=255)
                                                                                                tAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTc                                                               >  6:34202/1‑140 (MQ=255)
                                                                                                tAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtg                                                      >  1:105584/1‑149 (MQ=255)
                                                                                                 aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgt                                                     >  3:199685/1‑149 (MQ=255)
                                                                                                 aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgt                                                     >  3:116139/1‑149 (MQ=255)
                                                                                                    aGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGcc                                                  >  6:170736/1‑149 (MQ=255)
                                                                                                    aGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGcc                                                  >  8:303875/1‑149 (MQ=255)
                                                                                                     gCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGcc                                                  <  5:317710/148‑1 (MQ=255)
                                                                                                      ccGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGa                                                >  5:378824/1‑149 (MQ=255)
                                                                                                            agcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                          >  6:150521/1‑149 (MQ=255)
                                                                                                            agcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                          <  6:351931/149‑1 (MQ=255)
                                                                                                             gcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                  >  3:113441/1‑92 (MQ=255)
                                                                                                             gcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                  <  4:113441/92‑1 (MQ=255)
                                                                                                              cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                          <  1:147014/149‑1 (MQ=255)
                                                                                                              cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                        <  4:199685/149‑1 (MQ=255)
                                                                                                                 gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctggc                                     >  3:420913/1‑149 (MQ=255)
                                                                                                                 gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctggc                                     >  3:2418/1‑149 (MQ=255)
                                                                                                                   aCAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctggc                                     >  1:14920/1‑147 (MQ=255)
                                                                                                                    cAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAg                                  <  1:344167/149‑1 (MQ=255)
                                                                                                                     aGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgt                                 >  6:30109/1‑149 (MQ=255)
                                                                                                                       cGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtgg                               <  7:253332/149‑1 (MQ=255)
                                                                                                                         gcaggcatGGCTGCCCGAGCAGTCAA‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtggt                              >  5:279429/8‑148 (MQ=255)
                                                                                                                            cgcTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAAc                          <  4:420913/149‑1 (MQ=255)
                                                                                                                               atGGCTGCCCGAACAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGACATAGAAGCGCTGGTAGTCATGCGCTTGTGCGGACAGTGGCTGGCGAGTGGTGAACGgc                       >  3:35084/2‑149 (MQ=255)
                                                                                                                                  gctgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACa                                              >  7:297460/1‑123 (MQ=255)
                                                                                                                                  gctgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACa                                              <  8:297460/123‑1 (MQ=255)
                                                                                                                                  gctgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcgc                     >  3:123386/1‑148 (MQ=255)
                                                                                                                                   ctggtCGGGCAGGCAT‑‑GGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcgc                     >  7:383950/1‑147 (MQ=255)
                                                                                                                                    tggtCGGGCAGGCAT‑‑GGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                  >  2:420558/4‑149 (MQ=255)
                                                                                                                                    tgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                  >  8:176978/1‑149 (MQ=255)
                                                                                                                                    tgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                  >  4:323470/1‑149 (MQ=255)
                                                                                                                                    tgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                  >  6:427628/1‑149 (MQ=255)
                                                                                                                                        cGGGCAGGCAT‑‑GGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg              <  1:420558/149‑1 (MQ=255)
                                                                                                                                         gggCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCggg             >  1:211801/1‑149 (MQ=255)
                                                                                                                                          ggCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt            >  7:305997/1‑149 (MQ=255)
                                                                                                                                           gCAGGCAT‑‑GGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTg           >  6:185062/1‑149 (MQ=255)
                                                                                                                                             aGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgcg         <  7:281593/149‑1 (MQ=255)
                                                                                                                                              gCCAA‑‑GCCTGCCCGACCAGGTCGAACCCCGCAGCCTGAAAGCTACCCAGCCAGCTGAGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGaa                           >  7:239484/1‑130 (MQ=255)
                                                                                                                                               ccAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgcg         >  4:357202/1‑147 (MQ=255)
                                                                                                                                                 aT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAgg     >  6:117073/1‑149 (MQ=255)
                                                                                                                                                  t‑‑gcctgccCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGc    <  6:378824/149‑1 (MQ=255)
                                                                                                                                                      cctgccCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCgg  >  6:65815/1‑149 (MQ=255)
                                                                                                                                                      cctgccCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCgg  >  8:264188/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCGG  >  NC_002947/5886621‑5886917

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