Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,277,325 T→G 5.5% Q15H (CAA→CAC 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,3250TG 5.5% 231.9 / 3.6 73Q15H (CAA→CACPP_1116resolvase family site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base G (2/2);  total (33/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 = 8.25e-01

TCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCAAA  >  NC_002947/1277177‑1277473
                                                                                                                                                    |                                                                                                                                                    
tCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCt                                                                                                                                                      >  1:32577/1‑149 (MQ=255)
tCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCt                                                                                                                                                      >  7:247046/1‑149 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                     >  5:316801/1‑148 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                     >  5:205329/1‑148 (MQ=255)
  agagCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCtt                                                                                                                                                     >  3:27477/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                                                                                                                                                   <  4:367514/140‑1 (MQ=255)
            cACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGa                                                                                                                                                   >  3:367514/1‑140 (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                                                                                                                           <  5:230896/143‑1 (MQ=255)
                                 gagcgagcTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTaa                                                                                                                           >  6:230896/1‑143 (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                                                            <  5:50579/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                            <  5:256867/149‑1 (MQ=255)
                                                                                          cTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGa                                                            <  2:93272/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                                                                            <  1:157494/121‑1 (MQ=255)
                                                                                                      gTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGACTGTGTATGAGGCCAAt                                                                            >  2:157494/1‑121 (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                               >  5:223357/1‑127 (MQ=255)
                                                                                                                                             agTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTa                               <  6:223357/127‑1 (MQ=255)
                                                                                                                                                   cTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACATTCGTTAGTGCAAATGGAAATTCTAGATTTCCCa    >  3:363402/1‑148 (MQ=255)
                                                                                                                                                    ttGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCaaa  >  4:176192/1‑149 (MQ=255)
                                                                                                                                                    ttGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTAGTCAATTGGGCAGACGATCGTTAGTGCATATGGAAATTCTAGATTTGCCaaa  >  4:53801/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGAGCTTTTGCACCTTCTGGACGATCTTCAAGAGCGAGCTTGCCCGATACCACATCCACAAATTCACCAACGACATTCCAGTCGTTCGCTTGGGCCGCACGTTGGATGTAATCACGCTGGGCTTCCAGACCTAGGCCAGAGTCACCTTGACTCTTGGTAGATACTCGAATGTAAGCAACAACGTTTGTCATGATGATGGCTCGTCTGTGTATGAGGCCAATTTATCAACTGCCATGAGTATTTGTCAATTGGGCAGACGTTCGTTAGTGCAAATGGAAATTCTAGATTTCCCAAA  >  NC_002947/1277177‑1277473

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: