Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,733,265 T→A 9.7% intergenic (‑25/‑47) PP_3299 ← / → PP_3300 lipoprotein/Transcriptional regulator, tetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,733,2650TA 9.7% 167.5 / 8.1 62intergenic (‑25/‑47)PP_3299/PP_3300lipoprotein/Transcriptional regulator, tetR family
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (3/3);  total (33/29)
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.20e-01

GACGTAATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGG  >  NC_002947/3733118‑3733413
                                                                                                                                                   |                                                                                                                                                    
gACGTAATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCGGGGCCAACCATGCACCCTGCAGTGAGTCTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCatatt                                                                                                                                                     <  5:17853/149‑1 (MQ=255)
      aTGATCGGCCATCGCAACCTCTGGCTGTTTATAATGCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATAttt                                                                                                                                               <  4:342736/149‑1 (MQ=255)
         aTCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTg                                                                                                                                             <  7:373732/148‑1 (MQ=255)
             gCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTca                                                                                                                                         >  1:325343/1‑148 (MQ=255)
               cAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACGCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTGGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCacac                                                                                                                                      <  2:325343/149‑1 (MQ=255)
                  cgcAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCAAAATATAATATGGTCacacctc                                                                                                                                   <  7:205033/149‑4 (MQ=255)
                    cAACCTCTTGCTGTTTATAATCCGGTCCAACCATGCACCATTCATTGAGACTGTCGCTCAACAGCAGATCAATCGTTCTGTAGGGAAGCATGATGATTATCATAGACCACGAGGAGGTGTGACCATATTATATTTTGTTCACACTCTGt                                                                                                                                 <  5:43585/149‑1 (MQ=255)
                        ctctGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                             <  2:91871/149‑1 (MQ=255)
                                 tttATAATCCTGGCCAACCATGCACCCTTCAGTGAGACTGGCGCTAAACAGCAGAGCAATCGTTCGGTAGGGAATCATGATGATTCTCATAGACCAAGAGTAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCa                                                                                                                     <  7:135798/148‑1 (MQ=255)
                                 tttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAg                                                                                                                    >  5:105515/1‑149 (MQ=255)
                                 tttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAg                                                                                                                    <  6:161724/149‑1 (MQ=255)
                                 tttATAATCCGGCCCAACCATGCTCCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGTGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCa                                                                                                                     <  7:179645/148‑1 (MQ=255)
                                  ttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAg                                                                                                                    <  6:105515/148‑1 (MQ=255)
                                   tataATCGGGGCCCACCCTGCACCCTGCCGTGCGACTGGCGCTCAACAGCAGAGCAATCGTTCGGGAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAgg                                                                                                                   <  4:224751/148‑1 (MQ=255)
                                   tataATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGa                                                                                                                  >  5:277822/1‑149 (MQ=255)
                                   tataATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTAGTTTTTGGTCACACTCTGTCAATCAGTTCGAGGa                                                                                                                  >  6:232509/1‑149 (MQ=255)
                                     taATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATg                                                                                                                <  3:363489/149‑1 (MQ=255)
                                        tCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACActct                                                                                                                                   >  7:173351/1‑127 (MQ=255)
                                        tCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACActct                                                                                                                                   <  8:173351/127‑1 (MQ=255)
                                            ggCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCaaa                                                                                                         <  1:230762/149‑1 (MQ=255)
                                               cAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAAt                                                                                                      >  5:171038/1‑149 (MQ=255)
                                                 aCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGa                                                                                                    <  5:110866/149‑1 (MQ=255)
                                                  ccATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAg                                                                                                   >  3:52069/1‑149 (MQ=255)
                                                   cATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAg                                                                                                   >  7:164944/1‑148 (MQ=255)
                                                          ccTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAtctt                                                                                           >  2:74941/1‑149 (MQ=255)
                                                              cAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGAATCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTtctc                                                                                        <  4:102550/148‑1 (MQ=255)
                                                               aGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCa                                                                                      >  7:279534/1‑149 (MQ=255)
                                                               aGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTATGGTCACCCTCTGTCAGTCAGTTCCAGGATGCTTCAAAAATGCGTAAATGTTCTCCa                                                                                      >  8:300810/1‑149 (MQ=255)
                                                                gTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGATTCAAAAATGAGTAAATCTTCTCCAc                                                                                     >  4:221434/1‑149 (MQ=255)
                                                                 tGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCAcc                                                                                    <  6:171038/149‑1 (MQ=255)
                                                                   agaCTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTAGATTTTGGTCCCACTCTGTCAATCAGGTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCt                                                                                  >  6:351517/1‑149 (MQ=255)
                                                                               aaCAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGcgcgcg                                                                      >  7:13903/1‑149 (MQ=255)
                                                                                  agcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCAcc                                                                                    >  3:151291/1‑132 (MQ=255)
                                                                                  agcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCAcc                                                                                    <  4:151291/132‑1 (MQ=255)
                                                                                  agcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGc                                                                   >  5:410444/1‑149 (MQ=255)
                                                                                    cagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCa                                                                 >  5:2391/1‑149 (MQ=255)
                                                                                      gagCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGc                                                               >  6:128430/1‑149 (MQ=255)
                                                                                      gagCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCCGCGCGCCGGCGAgg                                                               <  6:277822/149‑2 (MQ=255)
                                                                                       agCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGcc                                                              <  6:84809/149‑1 (MQ=255)
                                                                                                   gTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGt                                                   >  6:93321/1‑148 (MQ=255)
                                                                                                       ggAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAgg                                               >  4:291767/1‑148 (MQ=255)
                                                                                                        gAATCATGATGATTCTCATAGACGAAGAGGAGGTGTGAGCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGga                                             <  5:232509/149‑1 (MQ=255)
                                                                                                         aaGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgat                                            >  4:360699/1‑149 (MQ=255)
                                                                                                            cATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcag                                         <  5:350743/149‑1 (MQ=255)
                                                                                                             atgatgATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGCCCACGACGTGAGGgat                                            >  7:79160/1‑145 (MQ=255)
                                                                                                             atgatgATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgat                                            >  3:417143/1‑145 (MQ=255)
                                                                                                               gatgatTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcagat                                       >  8:257442/1‑148 (MQ=255)
                                                                                                                  gatTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGatcatc                                   >  8:294149/1‑149 (MQ=255)
                                                                                                                      ctcATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAgg                               >  8:36060/1‑149 (MQ=255)
                                                                                                                         aTAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAgg                               <  8:130536/146‑1 (MQ=255)
                                                                                                                         aTAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTg                            >  5:312825/1‑149 (MQ=255)
                                                                                                                          tAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGc                           <  8:13903/149‑1 (MQ=255)
                                                                                                                            gACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCa                         >  1:35389/1‑149 (MQ=255)
                                                                                                                                aaGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGa                     <  2:35389/149‑1 (MQ=255)
                                                                                                                                  gaggagGGGTCACCACAGTCTATTTTGGTCACACTCGGTCACTCCGTTCCAGGATGCTGCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGc                   <  4:260178/149‑1 (MQ=255)
                                                                                                                                    cagagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcag                                         >  7:361527/5‑125 (MQ=39)
                                                                                                                                    cagagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcag                                         <  8:361527/121‑1 (MQ=39)
                                                                                                                                      gagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcagat                                       <  3:172892/123‑1 (MQ=255)
                                                                                                                                      gagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcagat                                       >  4:172892/3‑125 (MQ=255)
                                                                                                                                       ggTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTc              <  1:74941/149‑1 (MQ=255)
                                                                                                                                       agtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTc              >  2:209757/2‑149 (MQ=255)
                                                                                                                                             aCCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCAt        >  5:13265/1‑149 (MQ=255)
                                                                                                                                               cATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGACGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCAGCCAGGCTGCCATGGAGCACTTCGCCCATTa      >  8:250938/1‑149 (MQ=255)
                                                                                                                                                  attatattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATgg   >  4:27721/1‑149 (MQ=255)
                                                                                                                                                   ttatattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATggg  >  6:416236/1‑149 (MQ=255)
                                                                                                                                                   ttatattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATgg   <  8:279534/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACGTAATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGG  >  NC_002947/3733118‑3733413

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: