Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,000,740 T→C 5.5% intergenic (+210/‑1476) PP_4407 → / → PP_4409 gamma‑aminobutyrate transporter/site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,000,7400TC 5.5% 220.1 / 2.2 73intergenic (+210/‑1476)PP_4407/PP_4409gamma‑aminobutyrate transporter/site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (2/2);  total (36/37)
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

AAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCC  >  NC_002947/5000593‑5000888
                                                                                                                                                   |                                                                                                                                                    
aaGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTg                                                                                                                                                     >  7:69746/1‑149 (MQ=255)
aaGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTg                                                                                                                                                     <  6:92177/149‑1 (MQ=255)
    tGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGAc                                                                                                                                                  <  8:69746/148‑1 (MQ=255)
      ggAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTc                                                                                                                                                <  8:64079/148‑1 (MQ=255)
           gACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTg                                                                                                                                             >  8:249522/1‑146 (MQ=255)
              gagaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCa                                                                                                                                           <  8:119911/145‑1 (MQ=255)
              gagaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCa                                                                                                                                           >  7:119911/1‑145 (MQ=255)
                gaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCTGCAGGAGTCCAGACTCCTGCGGaaaac                                                                                                                                     <  1:178628/149‑4 (MQ=255)
                    ggCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCa                                                                                                                                 >  5:311399/1‑149 (MQ=255)
                      cGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCAtc                                                                                                                               >  5:289549/1‑149 (MQ=255)
                              gcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGa                                                                                                                       >  3:118260/1‑149 (MQ=255)
                              gcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGa                                                                                                                       <  2:210368/149‑1 (MQ=255)
                                aggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTATAGTTGTTTTCCGAAGGCTTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGagt                                                                                                                     >  7:244447/1‑147 (MQ=255)
                                    aTTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGt                                                                                                                 <  6:289549/149‑1 (MQ=255)
                                      tGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGc                                                                                                               >  2:10198/1‑149 (MQ=255)
                                       gAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGc                                                                                                               <  5:291879/148‑1 (MQ=255)
                                        aaCCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGctc                                                                                                             <  2:99168/149‑1 (MQ=255)
                                           ccTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTgg                                                                                                          <  1:204555/149‑1 (MQ=255)
                                           ccTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTgg                                                                                                          <  7:249686/149‑1 (MQ=255)
                                            cTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGc                                                                                                         <  4:280127/149‑1 (MQ=255)
                                                     ggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAAc                                                                                                <  1:10198/149‑1 (MQ=255)
                                                      gTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACt                                                                                                                     <  7:190165/127‑1 (MQ=255)
                                                      gTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACt                                                                                                                     >  8:190165/1‑127 (MQ=255)
                                                         cctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGAATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTg                                                                                                                                                     >  7:190029/1‑92 (MQ=255)
                                                         catcctGACTGATAAAGGTTCCGTGGTGAGGTATGAATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCtt                                                                                                                                                     <  8:190029/90‑2 (MQ=255)
                                                           tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCAt                                                                                                                                <  2:138576/111‑1 (MQ=255)
                                                           tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCAt                                                                                                                                >  1:138576/1‑111 (MQ=255)
                                                           tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          >  2:65739/1‑149 (MQ=255)
                                                           tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          >  1:340573/1‑149 (MQ=255)
                                                                 ctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCt                                                                                    >  8:186862/1‑149 (MQ=255)
                                                                   gATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTg                                                                                  >  5:207706/1‑149 (MQ=255)
                                                                    aTCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          >  3:159466/1‑140 (MQ=255)
                                                                    aTCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          <  4:159466/140‑1 (MQ=255)
                                                                    aTCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCtt                                                                                   <  8:295499/147‑1 (MQ=255)
                                                                     tCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTgc                                                                                 >  1:122691/1‑148 (MQ=255)
                                                                        aGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCtgt                                                                                                      <  2:177574/124‑1 (MQ=255)
                                                                         ggTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACATCTCCCATcc                                                                                     >  3:94977/1‑140 (MQ=255)
                                                                         ggTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACATCTCCCATcc                                                                                     <  4:94977/140‑1 (MQ=255)
                                                                         ggTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGcgccg                                                                             >  3:173058/1‑148 (MQ=255)
                                                                          gTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGTCTGTGTTAACAGCTCCCATCCTTGCGCCGCt                                                                           >  4:230870/1‑149 (MQ=255)
                                                                                tggtgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAg                                                                     <  1:65739/149‑1 (MQ=255)
                                                                                  gtgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAgg                                                                    <  6:207706/148‑1 (MQ=255)
                                                                                     cGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                <  8:150521/149‑1 (MQ=255)
                                                                                     cGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                <  6:311399/149‑1 (MQ=255)
                                                                                       gTATGCATCTGATTTTCCTTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGt                                                              <  8:298574/149‑1 (MQ=255)
                                                                                        tATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCt                                                                                                        <  2:70520/106‑1 (MQ=255)
                                                                                           gCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACg                                                           >  4:47833/1‑148 (MQ=255)
                                                                                            cATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGTGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGa                                                          >  1:282246/1‑148 (MQ=255)
                                                                                             aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCATGCAGAAGCTGCCATCTCAAGCGACAGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGCTTCTGTACGaaa                                                        >  7:190041/1‑149 (MQ=255)
                                                                                              tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAAt                                                       >  3:282209/1‑149 (MQ=255)
                                                                                                  aTTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCCGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAAt                                                       <  4:208939/145‑1 (MQ=255)
                                                                                                  aTTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGTGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATccc                                                    <  2:282246/148‑1 (MQ=255)
                                                                                                   ttttCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCga                                                  <  4:118260/149‑1 (MQ=255)
                                                                                                         gTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACt                                            <  7:186862/149‑1 (MQ=255)
                                                                                                               tagtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCGTCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATc                                       >  6:104696/1‑148 (MQ=255)
                                                                                                               tagtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATc                                       >  4:122409/1‑148 (MQ=255)
                                                                                                                 gtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCt                                                                                    >  5:241693/1‑101 (MQ=255)
                                                                                                                 gtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCt                                                                                    <  6:241693/101‑1 (MQ=255)
                                                                                                                  taCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACg                                   >  2:90788/1‑149 (MQ=255)
                                                                                                                  taCTTCTTCGTAGAGTTGTTTTCCCCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGAGACTGCGTGCTCTGGCTTTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTATGTACGAAATCACGAGAAACTGAATCGACg                                   >  1:327443/1‑149 (MQ=255)
                                                                                                                    cttcttCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGcc                                 >  8:154455/1‑149 (MQ=255)
                                                                                                                           gTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCgg                           <  2:122691/148‑1 (MQ=255)
                                                                                                                           gTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCgg                           <  7:249522/148‑1 (MQ=255)
                                                                                                                            tAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCGCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCgg                           >  1:247072/1‑147 (MQ=255)
                                                                                                                                gttgttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTg                     <  2:271621/149‑1 (MQ=255)
                                                                                                                                   agcttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCt                      <  3:310397/142‑1 (MQ=38)
                                                                                                                                   agcttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTGCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCt                      >  4:310397/4‑145 (MQ=38)
                                                                                                                                    ttttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTACCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTAGGGAGCTGAAAt                 >  5:136613/1‑149 (MQ=255)
                                                                                                                                        ccGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTg             >  6:271066/1‑149 (MQ=255)
                                                                                                                                         tgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTg             >  5:285552/2‑148 (MQ=255)
                                                                                                                                          gCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTa           <  4:173058/149‑1 (MQ=255)
                                                                                                                                             cgcgccTGGACCCCTGCCGAAGCTCCCCTCGCAAGCGACTGCGTGCTCTGGCTGTGTACACAGCTCCCATCCTTGCGCCGCTCCCTCGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAggg        <  2:340573/144‑1 (MQ=25)
                                                                                                                                              gAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGAGGCCCCTCGGGAGCTGAAATTCTGTAgg         <  5:176641/147‑1 (MQ=255)
                                                                                                                                              gAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAgg         >  6:176641/1‑147 (MQ=255)
                                                                                                                                                   tGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTcc  >  7:220602/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGATGGGAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCC  >  NC_002947/5000593‑5000888

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: