Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,045,934 C→A 8.0% Q262H (CAG→CAT) ‡ PP_3567 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,045,9340CA 8.0% 214.3 / 7.8 75Q262H (CAG→CAT) ‡PP_3567LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base A (3/3);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

TGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCA  >  NC_002947/4045798‑4046082
                                                                                                                                         |                                                                                                                                                    
tGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCggg                                                                                                                                           <  5:26574/148‑1 (MQ=255)
         cgcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                 >  6:80068/1‑149 (MQ=255)
           cAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTg                                                                                                                               >  2:249851/1‑149 (MQ=255)
            aaCACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGc                                                                                                                              <  1:249851/149‑1 (MQ=255)
                 tAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGc                                                                                                                                                   >  8:31321/1‑123 (MQ=255)
                 tAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGc                                                                                                                                                   <  7:31321/123‑1 (MQ=255)
                 tAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATTTTCGCTGCAGgcg                                                                                                                         >  7:120092/1‑149 (MQ=255)
                 tAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGgcg                                                                                                                         >  4:212485/1‑149 (MQ=255)
                   gTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGcc                                                                                                                       >  6:86561/1‑149 (MQ=255)
                      aCCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGTGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCt                                                                                                                     <  7:332596/148‑1 (MQ=255)
                      aCCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTg                                                                                                                    >  4:275896/1‑149 (MQ=255)
                            tATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAaca                                                                                                              >  4:8589/1‑149 (MQ=255)
                              tGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAAcaca                                                                                                            >  7:217068/1‑149 (MQ=255)
                              tGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcctggc                                                                                                                                         >  4:11566/1‑115 (MQ=255)
                              tGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcctggc                                                                                                                                         <  3:11566/120‑6 (MQ=255)
                              tGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGCCCTGCCAGTCGGGCAATATACGCTGCAGGTGCCCTGCCAacac                                                                                                             >  8:31423/1‑148 (MQ=255)
                                aGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCt                                                                                                                     >  7:308549/1‑138 (MQ=255)
                                aGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCt                                                                                                                     <  8:308549/138‑1 (MQ=255)
                                aGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACAt                                                                                                           <  1:91661/148‑1 (MQ=255)
                                aGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATc                                                                                                          >  4:294007/1‑149 (MQ=255)
                                  aaaaTTTTCGCAAAAGTCGATGAATGCGTTTACACGCCGTGAGCCCCGGTCGTTGGGCATATACAGGGCGTTGATTCTGCTATTGGCTGCGCCAGGACTGGCCTGCCACTCGGGGAAAATTCGCTCCAGGAGCGCTGCCAACACATcca                                                                                                        >  4:317348/2‑148 (MQ=255)
                                     aTTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGa                                                                                                                                         <  3:340324/113‑1 (MQ=255)
                                     aTTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGa                                                                                                                                         >  4:340324/1‑113 (MQ=255)
                                           gCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAcc                                                                                                <  5:86561/148‑1 (MQ=255)
                                           gCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcctg                                                                                                                                           <  4:299373/105‑4 (MQ=255)
                                           gCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcctg                                                                                                                                           >  3:299373/1‑102 (MQ=255)
                                             aaaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                             >  6:340348/1‑149 (MQ=255)
                                             aaaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                             <  3:212485/149‑1 (MQ=255)
                                             aaaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                             <  8:30174/149‑1 (MQ=255)
                                             aaaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGATGCGCCAGGACTGGCCTGCAAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCACGACATCCTCGCGTAccagc                                                                                             >  3:3799/1‑149 (MQ=255)
                                               atgTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGTGCCCTGCCAACACATCCTCGCGTAccagcca                                                                                           <  7:31423/147‑1 (MQ=255)
                                                    gATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                       <  1:142752/148‑1 (MQ=255)
                                                     aTGAATGCGTTTACCCGCCGTGAGTCCCCGGTGGTTGGGCATATACAGGGCGTTGATGATGCTATTGGCTGCGCCAGGACTGGCGTGCCAGTCGGGGATTATTCGATGCAGGCGCCGTGCCAACACATCATCGCGTACCAGCCAGTCgg                                                                                      >  7:52182/1‑149 (MQ=255)
                                                     aTGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAac                                                                                                               <  1:295070/123‑1 (MQ=255)
                                                     aTGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAac                                                                                                               >  2:295070/1‑123 (MQ=255)
                                                             gTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTg                                                                             >  1:110104/1‑149 (MQ=255)
                                                             gTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTg                                                                             >  1:190107/1‑149 (MQ=255)
                                                                  cccGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGAt                                                                         <  5:80068/148‑1 (MQ=255)
                                                                   ccgccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                             >  5:199339/1‑127 (MQ=255)
                                                                   ccgccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                             <  6:199339/127‑1 (MQ=255)
                                                                        gTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActg                                                                  <  2:227144/149‑1 (MQ=255)
                                                                               cccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcctgc                                                             >  7:293757/1‑148 (MQ=255)
                                                                                ccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcctgc                                                             <  4:319030/147‑1 (MQ=255)
                                                                                ccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcctgc                                                             >  7:151263/1‑147 (MQ=255)
                                                                                ccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcctgc                                                             >  3:319030/1‑147 (MQ=255)
                                                                                acggcggTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCa                                                           <  5:279607/148‑1 (MQ=255)
                                                                                 cGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAg                                                          >  1:327526/1‑149 (MQ=255)
                                                                                      gtTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCg                                                     >  7:212874/1‑149 (MQ=255)
                                                                                      gtTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCg                                                     >  6:7371/1‑149 (MQ=255)
                                                                                      gtTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTACCg                                                     >  6:307729/1‑149 (MQ=255)
                                                                                        tGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGtt                                                   >  5:270055/1‑149 (MQ=255)
                                                                                         gggCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGttt                                                  >  4:129896/1‑149 (MQ=255)
                                                                                          ggCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTc                                                 <  3:294007/149‑1 (MQ=255)
                                                                                           gCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCa                                                <  8:212874/149‑1 (MQ=255)
                                                                                           gCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCa                                                <  3:317348/149‑1 (MQ=255)
                                                                                            cATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCAc                                               <  2:110104/149‑1 (MQ=255)
                                                                                             atatACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACg                                              >  3:271909/1‑149 (MQ=255)
                                                                                                taCAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCaa                                           >  4:345531/1‑149 (MQ=255)
                                                                                                       cGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCg                                    <  3:8589/149‑1 (MQ=255)
                                                                                                              gctgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTg                              <  8:151263/148‑1 (MQ=255)
                                                                                                                  ctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGgcg                           <  8:217068/147‑1 (MQ=255)
                                                                                                                  ctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCt                         <  1:23445/149‑1 (MQ=255)
                                                                                                                   tATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTg                        >  2:147009/1‑149 (MQ=255)
                                                                                                                          tGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGc                                                                                     <  1:181468/81‑1 (MQ=255)
                                                                                                                          tGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGc                                                                                     >  2:181468/1‑81 (MQ=255)
                                                                                                                            cgcCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc               <  6:270055/149‑1 (MQ=255)
                                                                                                                               cccgACTGGCAGGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc               <  5:213385/143‑1 (MQ=255)
                                                                                                                               cccgACTGGCAGGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc               >  6:213385/4‑146 (MQ=255)
                                                                                                                               cAGGACTGGCCTGCCAGTCGGTGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTAGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCATTGCCGTTTCACTCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCa            <  8:293757/149‑1 (MQ=255)
                                                                                                                               cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTg                 <  6:38568/144‑1 (MQ=255)
                                                                                                                               cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTg                 >  5:38568/1‑144 (MQ=255)
                                                                                                                               cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCa            >  6:30246/1‑149 (MQ=255)
                                                                                                                               cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCa            >  6:280958/1‑149 (MQ=255)
                                                                                                                                   aCTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc               >  8:105100/1‑142 (MQ=255)
                                                                                                                                   aCTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc               <  7:105100/142‑1 (MQ=255)
                                                                                                                                       gacTGCGAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAAGACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGt     <  8:120092/146‑1 (MQ=255)
                                                                                                                                         cTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCa  <  3:345531/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCA  >  NC_002947/4045798‑4046082

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: