Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,000,741 G→A 13.0% intergenic (+211/‑1475) 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,7410GA13.0% 150.6 / 11.7 63intergenic (+211/‑1475)PP_4407/PP_4409gamma‑aminobutyrate transporter/site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (3/5);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCC  >  NC_002947/5000600‑5000888
                                                                                                                                             |                                                                                                                                                   
gAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCt                                                                                                                                              >  5:89355/1‑149 (MQ=255)
 aaTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTg                                                                                                                                             >  8:193070/1‑149 (MQ=255)
  aTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGc                                                                                                                                            <  2:225188/149‑1 (MQ=255)
  aTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGc                                                                                                                                            <  2:230249/149‑1 (MQ=255)
  aTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCCAGACTCCTGc                                                                                                                                            <  6:308067/149‑1 (MQ=255)
  aTGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCCAGACTCCTGc                                                                                                                                            <  8:162401/149‑1 (MQ=255)
         gaATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCt                                                                                                                                     >  8:239294/1‑149 (MQ=255)
           aTGGCGTCCCAGTCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTg                                                                                                                                    <  6:57241/148‑1 (MQ=255)
           aTGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTACGTAGAGTTGTTTTCTGCAGGAGTCCAGACTCCTGCGGaaaacaa                                                                                                                                   <  1:255129/149‑6 (MQ=255)
            tGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGcc                                                                                                                                  <  6:89355/149‑1 (MQ=255)
              gCGTCCCAGGCAGGAATTTAACCATTATCTGGGGTGCCTCCTGCCACATCAAGGTGCCGTTGTGCGGTATTCATCTCAATTTCCGAAGAGAAGTACTTCTACGTATAGTAGATTTCCGCAGCAGTCTCGACACCTGCAGAAGCAGCCa                                                                                                                                 >  8:63867/1‑148 (MQ=255)
               cGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCAtc                                                                                                                               <  7:193070/149‑1 (MQ=255)
               cGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCAtc                                                                                                                               <  6:248087/149‑1 (MQ=255)
                       gcaggAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGa                                                                                                                       >  3:80694/1‑149 (MQ=255)
                        caggTATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGAc                                                                                                                      <  7:63867/149‑1 (MQ=255)
                           gAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGc                                                                                                                   >  1:276376/1‑149 (MQ=255)
                            aaTTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGc                                                                                                                   <  7:135846/148‑1 (MQ=255)
                                 aaCCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGctc                                                                                                             >  2:201638/1‑149 (MQ=255)
                                             gggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTgg                                                                                                          >  5:123319/1‑140 (MQ=255)
                                             gggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTgg                                                                                                          <  6:123319/140‑1 (MQ=255)
                                             gggTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTaa                                                                                                 <  7:239294/149‑1 (MQ=255)
                                                   ctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGTCTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGGTGTGTTAACAGCTc                                                                                           <  2:276376/149‑1 (MQ=255)
                                                   ctcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTc                                                                                           >  4:158894/1‑149 (MQ=255)
                                                    tcctGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          >  6:370803/1‑149 (MQ=255)
                                                     cctgcctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTcc                                                                                          >  5:20034/1‑148 (MQ=255)
                                                      ctgcctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTgtgt                                                                                                    <  3:18550/137‑1 (MQ=255)
                                                      ctgcctgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTgtgt                                                                                                    >  4:18550/1‑137 (MQ=255)
                                                           tgATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCtt                                                                                   >  1:282170/1‑149 (MQ=255)
                                                             ttCATGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTGTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCTTGCTCTGGCTGTGTTAACAGCTCCCATCCTTgc                                                                                 <  1:200503/148‑1 (MQ=255)
                                                             aTCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCtt                                                                                   <  4:80694/147‑1 (MQ=255)
                                                             aTCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTgc                                                                                 >  1:83866/1‑149 (MQ=255)
                                                              tCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTc                                                                                           <  8:250039/138‑1 (MQ=255)
                                                                            tgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGtt                                                                  <  2:282170/149‑1 (MQ=255)
                                                                            tgCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGtt                                                                  <  6:20034/149‑1 (MQ=255)
                                                                              cGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCt                                                                >  5:39573/1‑149 (MQ=255)
                                                                                      aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACg                                                           <  1:86512/146‑1 (MQ=255)
                                                                                      aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACg                                                           >  2:86512/1‑146 (MQ=255)
                                                                                      aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGaaa                                                        >  3:118578/1‑149 (MQ=255)
                                                                                      aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGaaa                                                        >  7:150690/1‑149 (MQ=255)
                                                                                      aTCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGAAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAAAAGCTCACATCCTTGCCCCGCTCCATAGGTTCTGTAAGaaa                                                        >  8:274028/1‑149 (MQ=255)
                                                                                       tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAAt                                                       >  7:205233/1‑149 (MQ=255)
                                                                                       tCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGATTTCCGCAGGAGTATGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTTTGTTAAAAGATACCATCCTTGCGACGCTCACTAGGTTCTGTAAGAAAt                                                       >  7:199899/1‑149 (MQ=255)
                                                                                          gATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATccc                                                    >  8:15341/1‑149 (MQ=255)
                                                                                                     gagTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGaa                                         <  6:201895/149‑1 (MQ=255)
                                                                                                       gtagtaCTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAAt                                        <  8:150690/148‑1 (MQ=255)
                                                                                                                    gTAGTGTTGTTTTCCGCAGGATTCTTGACTCCTTCAGAAGCTGCCATCTCAATCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGATACTGAATCGACGCCCCTCtgg                          <  8:199899/149‑1 (MQ=255)
                                                                                                                    gTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGATTCGACGCCCCTCggg                          <  8:205233/149‑1 (MQ=255)
                                                                                                                     tAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTcc                                                                         <  4:307664/101‑1 (MQ=255)
                                                                                                                     tAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGa                         <  2:83866/149‑1 (MQ=255)
                                                                                                                         gttgttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCt                      >  4:18408/1‑148 (MQ=255)
                                                                                                                           cagcttCTGCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGa                    >  3:25047/5‑148 (MQ=38)
                                                                                                                            gttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAg                                                                     <  8:89622/98‑1 (MQ=255)
                                                                                                                            gttTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAg                                                                     >  7:89622/1‑98 (MQ=255)
                                                                                                                              tttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAAtt                >  1:142127/1‑149 (MQ=255)
                                                                                                                               ttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAAtt                <  3:18408/148‑1 (MQ=255)
                                                                                                                                 ctgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCt              <  4:25047/146‑1 (MQ=255)
                                                                                                                                  tgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGa                         <  3:280517/135‑1 (MQ=255)
                                                                                                                                  tgCAGGAGTCCAGACTCCTGCGGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGa                         >  4:280517/2‑136 (MQ=255)
                                                                                                                                  tgCAGGAGTCCAGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTg             >  7:60612/2‑148 (MQ=255)
                                                                                                                                   gCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTa           <  7:15341/149‑1 (MQ=255)
                                                                                                                                     aGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTCCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAgg         >  1:255303/1‑149 (MQ=255)
                                                                                                                                            tGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTcc  >  1:278362/1‑149 (MQ=255)
                                                                                                                                            tGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTcc  >  8:6103/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GAATGACGAGAATGGCGTCCCAGGCAGGAATTGAACCCTTATCTGGGGTGCCTCCTGCCTGATCAAGGTGCCGTGGTGCGGTATGCATCTGATTTTCCGTAGAGTAGTACTTCTTCGTAGAGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCC  >  NC_002947/5000600‑5000888

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