Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,886,764 C→G 5.8% W84C (TGG→TGC)  PP_5161 ← lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,886,7640CG 5.8% 263.4 / 3.8 83W84C (TGG→TGC) PP_5161lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (41/36);  new base G (3/2);  total (45/38)
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.08e-01

CGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCG  >  NC_002947/5886619‑5886910
                                                                                                                                                 |                                                                                                                                                    
cgcTGGACAGGCCGTGGCGGGGCAGGTCTCTGCTGATCACGGCGAAACCCTGCTTGAGTGCCGACTAGATCACATGTCGTTAGAGGCCCATGTGGGCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                   <  5:121571/149‑1 (MQ=255)
cgcTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                   >  3:425755/1‑149 (MQ=255)
  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                                                                                           <  4:371900/148‑1 (MQ=255)
                                                       gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                           <  7:205687/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                                                          >  1:413228/1‑142 (MQ=255)
                                                                                              gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAg                                                          <  2:413228/142‑1 (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                                       <  8:297460/123‑1 (MQ=255)
                                                                                                                                    gctgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACa                                       >  7:297460/1‑123 (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           >  6:427628/1‑149 (MQ=255)
                                                                                                                                      tgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg           >  8:176978/1‑149 (MQ=255)
                                                                                                                                      tgctCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg           >  4:323470/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)
                                                                                                                                                 |                                                                                                                                                    
CGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAT‑‑GCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCG  >  NC_002947/5886619‑5886910

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