Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,277,324 C→G 5.1% G16R (GGT→CGT) ‡ PP_1116 ← resolvase family site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,277,3240CG 5.1% 246.7 / 3.2 78G16R (GGT→CGT) ‡PP_1116resolvase family site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCA  >  NC_002947/1277177‑1277471
                                                                                                                                                   |                                                                                                                                                   
tCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAcc                                                                                                                                                     <  8:339785/148‑1 (MQ=255)
tCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCt                                                                                                                                                    >  1:32577/1‑149 (MQ=255)
tCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCt                                                                                                                                                    >  7:247046/1‑149 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                   >  5:205329/1‑148 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                   >  3:27477/1‑148 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                   >  5:316801/1‑148 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                   >  2:365108/1‑148 (MQ=255)
       ttttGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGActct                                                                                                                                             >  5:99630/1‑149 (MQ=255)
         ttGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCAAGGTGActctg                                                                                                                                            <  6:212477/148‑2 (MQ=255)
          tGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTg                                                                                                                                           <  1:95126/148‑1 (MQ=255)
           gCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                          <  7:252252/148‑1 (MQ=255)
            cACCTTCTTGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTa                                                                                                                                        <  4:279967/149‑1 (MQ=255)
            cACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGa                                                                                                                                                 >  3:367514/1‑140 (MQ=255)
            cACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGa                                                                                                                                                 <  4:367514/140‑1 (MQ=255)
            cACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGt                                                                                                                                         <  2:32577/148‑1 (MQ=255)
             aCCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                          >  2:228105/1‑146 (MQ=255)
             aCCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                          >  8:143654/1‑146 (MQ=255)
             aCCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTgg                                                                                                                                          >  8:241255/1‑146 (MQ=255)
               cTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGa                                                                                                                                      >  1:329572/1‑148 (MQ=255)
                  cTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATa                                                                                                                                    <  6:242380/147‑1 (MQ=255)
                   tGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTc                                                                                                                                 >  6:163117/1‑149 (MQ=255)
                    ggACGCTCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTc                                                                                                                                 >  7:82979/1‑148 (MQ=255)
                     gACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTc                                                                                                                                 <  1:128713/147‑1 (MQ=255)
                     gACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCAAAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCg                                                                                                                                >  6:329160/1‑148 (MQ=255)
                      aCGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGaa                                                                                                                              <  6:153502/149‑1 (MQ=255)
                      aCGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGaa                                                                                                                              <  4:7086/149‑1 (MQ=255)
                        gATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATg                                                                                                                            >  6:345416/1‑149 (MQ=255)
                              cAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATg                                                                                                                            <  7:241255/143‑1 (MQ=255)
                              cAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                         <  1:228105/146‑1 (MQ=255)
                                 gagcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                         >  6:230896/1‑143 (MQ=255)
                                 gagcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                         <  5:230896/143‑1 (MQ=255)
                                 gagcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGcaac                                                                                                                    >  8:211718/1‑148 (MQ=255)
                                   gcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCaacaac                                                                                                                 <  5:329160/149‑1 (MQ=255)
                                       gcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGttt                                                                                                             <  7:143654/149‑1 (MQ=255)
                                       gcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGttt                                                                                                             <  1:365108/149‑1 (MQ=255)
                                         ttGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGt                                                                                                           >  1:107636/1‑149 (MQ=255)
                                              cGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCatga                                                                                                      <  6:316801/149‑1 (MQ=255)
                                              cGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCatga                                                                                                      <  8:247046/149‑1 (MQ=255)
                                                   ccACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCatgatgat                                                                                                  <  6:205329/148‑1 (MQ=255)
                                                    cacaTAAACAAATTCACCAACGACATTCCAGTCGTTCGATTGGGCCGCACGTTGGATGTAATCACGCTCGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATAATCGAATGTAAGCAACAACGTTTGTCatgatgatg                                                                                                 >  3:342907/1‑148 (MQ=255)
                                                     acaTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGc                                                                                                                       <  7:63397/125‑1 (MQ=255)
                                                     acaTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGc                                                                                                                       >  8:63397/1‑125 (MQ=255)
                                                      caTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGa                                                                                                                                      >  7:314372/1‑109 (MQ=255)
                                                      caTCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGa                                                                                                                                      <  8:314372/109‑1 (MQ=255)
                                                          cacaAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGt                                                                                           <  2:329572/148‑1 (MQ=255)
                                                                 tCACCAACGACATTCCAGTCGTTCGCTAGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCATGGTAGATACTCGAATGTAAGCAACAACGGTTGTCAGGATGATGGCTCGTCTGTGAAt                                                                                   >  4:220174/1‑149 (MQ=255)
                                                                   aCCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTTACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGa                                                                                 <  5:300358/149‑1 (MQ=255)
                                                                   aCCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGa                                                                                 <  4:27477/149‑1 (MQ=255)
                                                                   aCCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGa                                                                                 <  6:99630/149‑1 (MQ=255)
                                                                      aaCGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGACTGTGTATGAGGc                                                                              >  3:136381/1‑149 (MQ=255)
                                                                         gACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCaa                                                                           >  1:93272/1‑149 (MQ=255)
                                                                               ccAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTAt                                                                      <  7:211718/148‑1 (MQ=255)
                                                                                     gTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGc                                                               <  8:82979/149‑1 (MQ=255)
                                                                                         gCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCAt                                                            >  6:256867/1‑148 (MQ=255)
                                                                                          cTTGGGCCTCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                          <  3:332506/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCAt                                                            >  4:9533/1‑147 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                          <  2:93272/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                          <  5:50579/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                          <  5:256867/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                          <  3:9533/149‑1 (MQ=255)
                                                                                               gCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTAtt                                                     >  2:98770/1‑149 (MQ=255)
                                                                                                 cGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTATGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTg                                                   >  5:49877/1‑149 (MQ=255)
                                                                                                    aCGATGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTTTGAGGCCAAATTATCAACTGCCATGAGTATTTGTc                                                 >  2:20123/1‑148 (MQ=255)
                                                                                                      gTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGACTGTGTATGAGGCCAAt                                                                          >  2:157494/1‑121 (MQ=255)
                                                                                                      gTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGACTGTGTATGAGGCCAAt                                                                          <  1:157494/121‑1 (MQ=255)
                                                                                                       ttGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAAtt                                             >  5:163581/1‑149 (MQ=255)
                                                                                                                     gCTGGGCTTCCAGACCTATGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATTATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATTAGTATTTGTCAATTTGGCAGActttcgt                               <  4:342907/149‑1 (MQ=255)
                                                                                                                      cTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTATGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGAcgttcgtt                              <  6:49877/149‑1 (MQ=255)
                                                                                                                          gCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGt                           <  5:163117/148‑1 (MQ=255)
                                                                                                                            ttCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATTAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGc                         >  7:54177/1‑148 (MQ=255)
                                                                                                                            ttCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGc                         >  7:55923/1‑148 (MQ=255)
                                                                                                                                 gttcTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTTTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATgg                   <  5:91087/146‑1 (MQ=255)
                                                                                                                                       ggCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTAAATTGGGCAGAAGTTCGTTAGTGCAAATgg                   >  6:91087/1‑143 (MQ=255)
                                                                                                                                         caagagTCAAGGTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGaaa                <  8:1262/142‑1 (MQ=255)
                                                                                                                                         caagagTCAAGGTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGaaa                >  7:1262/3‑144 (MQ=255)
                                                                                                                                         caagagTCAAGGTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTa           >  3:55753/3‑149 (MQ=255)
                                                                                                                                          cAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAg          <  7:35113/149‑1 (MQ=255)
                                                                                                                                             agTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTa                             <  6:223357/127‑1 (MQ=255)
                                                                                                                                             agTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTa                             >  5:223357/1‑127 (MQ=255)
                                                                                                                                                   cTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACATTCGTTAGTGCAAATGGAAATTCTAGATTTCCCa  >  3:363402/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCA  >  NC_002947/1277177‑1277471

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: