Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,476,725 G→A 5.5% intergenic (‑13/+155) PP_3087 ← / ← PP_3088 excinuclease ABC subunit A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,476,7250GA 5.5% 225.6 / 2.5 74intergenic (‑13/+155)PP_3087/PP_3088excinuclease ABC subunit A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TTGCCCGAGCCGGACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAGGCTTTTA  >  NC_002947/3476579‑3476870
                                                                                                                                                  |                                                                                                                                                 
ttGCCCGAGCCGGACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTc                                                                                                                                                 >  7:197974/1‑149 (MQ=255)
           ggACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAg                                                                                                                                       >  4:323709/1‑148 (MQ=255)
            gACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAg                                                                                                                                       >  5:70920/1‑147 (MQ=255)
            gACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCa                                                                                                                                     <  3:28906/149‑1 (MQ=255)
               acGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTATCTGAAGCTGAGCAGTg                                                                                                                                  >  5:30295/1‑149 (MQ=255)
               acGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTATCTGAAGCTGAGCAGTg                                                                                                                                  >  5:24955/1‑149 (MQ=255)
                 gCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTgcg                                                                                                                                >  3:150592/1‑149 (MQ=255)
                    ggTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTc                                                                                                                             <  8:231914/149‑1 (MQ=255)
                     gTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCa                                                                                                                            <  8:290275/149‑1 (MQ=255)
                      tGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCaa                                                                                                                           <  1:269618/149‑1 (MQ=255)
                      tGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCaa                                                                                                                           <  5:65470/149‑1 (MQ=255)
                        tgcaccacCAAGGCGGCGCGGGGGATCTCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAATCCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGt                                                                                                                         <  4:267657/147‑1 (MQ=255)
                         acacCACCAATGCTGCGCGGGTGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGtt                                                                                                                        <  6:256744/149‑1 (MQ=255)
                          caccacCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTc                                                                                                                                                 <  5:263141/123‑1 (MQ=255)
                          caccacCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTc                                                                                                                                                 >  6:263141/1‑123 (MQ=255)
                           accaccAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTaa                                                                                                                      <  5:240026/149‑1 (MQ=255)
                             caccaAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGTGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTaa                                                                                                                      >  5:111455/1‑147 (MQ=255)
                             caccaAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACa                                                                                                                    <  2:235193/149‑1 (MQ=255)
                                   ggCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAAtt                                                                                                              >  6:73168/1‑149 (MQ=255)
                                    gCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTg                                                                                                             <  1:95157/149‑1 (MQ=255)
                                      gTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTg                                                                                                             >  6:60437/1‑147 (MQ=255)
                                      gTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGCGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCa                                                                                                           >  7:272282/1‑149 (MQ=255)
                                             gggATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCt                                                                                                    >  3:167567/1‑149 (MQ=255)
                                              ggATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCtt                                                                                                   <  5:105479/149‑1 (MQ=255)
                                                      aCATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAAt                                                                                                               <  8:239402/129‑1 (MQ=255)
                                                      aCATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAAt                                                                                                               >  7:239402/1‑129 (MQ=255)
                                                       cATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGgca                                                                                          <  2:94686/149‑1 (MQ=255)
                                                       cATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGgca                                                                                          <  4:167567/149‑1 (MQ=255)
                                                               ttcttcAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCaa                                                                                  <  3:323709/149‑1 (MQ=255)
                                                                 cttcAGGTTGTGCTCGCGTGCGCCAAGGACGTGGATGAATCGGCTTGGGGTGGTCGGGCTTTGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAg                                                                                 <  8:272282/148‑1 (MQ=255)
                                                                 cttcAGGTTGTGCTCGCGGGTGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAg                                                                                 <  6:111455/148‑1 (MQ=255)
                                                                 cttcAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAg                                                                                 <  6:70920/148‑1 (MQ=255)
                                                                   tcAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGaa                                                                              >  3:169027/1‑149 (MQ=255)
                                                                    cAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAAt                                                                             >  4:357970/1‑149 (MQ=255)
                                                                     aGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATc                                                                            >  8:15021/1‑149 (MQ=255)
                                                                     aGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATc                                                                            >  2:280343/1‑149 (MQ=255)
                                                                       gTTGTGCTCGCGGGCGCCACGGACGCGGATGATACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCat                                                                                                                <  7:239857/111‑2 (MQ=255)
                                                                       gTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCaa                                                                                                                >  8:239857/1‑111 (MQ=255)
                                                                       gTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTa                                                                          >  3:302042/1‑149 (MQ=255)
                                                                        ttGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTAc                                                                         <  5:60437/149‑1 (MQ=255)
                                                                          gtgCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTAcc                                                                        <  5:73168/148‑1 (MQ=255)
                                                                               cgtgggCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGa                                                                     >  7:127104/4‑146 (MQ=255)
                                                                                 cgGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATagag                                                                >  7:154655/1‑149 (MQ=255)
                                                                                      gcCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggt                                                           <  4:169027/149‑1 (MQ=255)
                                                                                       ccACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGATGCATGAATCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtg                                                          <  6:219873/149‑1 (MQ=255)
                                                                                       ccACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtg                                                          <  7:15021/149‑1 (MQ=255)
                                                                                        cACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtgg                                                         >  6:304951/1‑149 (MQ=255)
                                                                                          cGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtgggg                                                       <  1:280343/149‑1 (MQ=255)
                                                                                          cGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtgggg                                                       <  4:150592/149‑1 (MQ=255)
                                                                                              cgcgGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCggtggggtgg                                                    >  6:54612/1‑148 (MQ=255)
                                                                                                      aaaCCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGcc                                           >  2:152886/1‑149 (MQ=255)
                                                                                                        aCCGCTTGGGGTGGTCTGTCGTGGAGGCATGAAGCATGCCTTTTCTGAATCTGAGCAGTGCTCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGCGCGCCgg                                         <  8:154655/149‑1 (MQ=255)
                                                                                                               ggggTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTc                                    <  5:54612/147‑1 (MQ=255)
                                                                                                               ggggTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAAATGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCtt                                  >  2:357678/1‑149 (MQ=255)
                                                                                                                 ggtggtCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCg                                >  4:36787/1‑149 (MQ=255)
                                                                                                                  gtggtCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGc                               <  1:152886/149‑1 (MQ=255)
                                                                                                                  gtggtCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGc                               <  8:127104/149‑1 (MQ=255)
                                                                                                                        gggCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTaa                                                                                                                      >  4:311638/1‑56 (MQ=255)
                                                                                                                        gggCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTaa                                                                                                                      <  3:311638/56‑1 (MQ=255)
                                                                                                                        gggCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTAc                          <  3:36787/148‑1 (MQ=255)
                                                                                                                         ggCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCt                        <  4:302042/149‑1 (MQ=255)
                                                                                                                               gAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCt                        >  7:126305/1‑143 (MQ=255)
                                                                                                                               gAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCt                        <  8:126305/143‑1 (MQ=255)
                                                                                                                                aGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGc                               <  8:47210/135‑1 (MQ=255)
                                                                                                                                aGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGc                               >  7:47210/1‑135 (MQ=255)
                                                                                                                                aGGCATGAAGCAGTCCTTATCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTc                 <  6:30295/149‑1 (MQ=255)
                                                                                                                                aGGCATGAAGCAGTCCTTATCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTc                 <  6:24955/149‑1 (MQ=255)
                                                                                                                                    aTGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGtct             >  4:5407/1‑149 (MQ=255)
                                                                                                                                      gAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGtctc            >  4:27286/1‑148 (MQ=255)
                                                                                                                                         gCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAgg        >  2:194416/1‑149 (MQ=255)
                                                                                                                                           aGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCACCGCCCTTGCGTGGGCAGCAGTCAAGGAATCAACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGATCCTCTCAAGGCt      >  3:58193/1‑149 (MQ=255)
                                                                                                                                               cTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAGGCTTTTa  >  2:331804/1‑149 (MQ=255)
                                                                                                                                                ttGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAGGCtttt   >  4:227394/1‑147 (MQ=255)
                                                                                                                                                ttGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAGGCtttt   <  3:227394/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGCCCGAGCCGGACACGCCGGTGAACACCACCAAGGCGTCGCGGGGGATATCGACATCGATGTTCTTCAGGTTGTGCTCGCGGGCGCCACGGACGCGGATGAAACCGCTTGGGGTGGTCGGGCTTGGAGGCATGAAGCAGTCCTTGTCTGAAGCTGAGCAGTGCGCTCAAGTTAACAGCAATTGCAGCGCCCTTGCGTGGGCAGCAGTCAAGGAATCTACCCGATAGAGCCGGTGGGGTGGCAGGGCGCCGGCTGTCTTCGCCTTACCTGACGTTCGTCTCAAGGCTTTTA  >  NC_002947/3476579‑3476870

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