Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,733,270 T→A 9.8% intergenic (‑30/‑42) 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,2700TA 9.8% 166.9 / 8.1 61intergenic (‑30/‑42)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/28)
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.44e-01

ATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGG  >  NC_002947/3733124‑3733413
                                                                                                                                                  |                                                                                                                                               
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                                                                                                                    <  6:161724/149‑1 (MQ=255)
                           tttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAg                                                                                                                    >  5:105515/1‑149 (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                                                                                    <  4:151291/132‑1 (MQ=255)
                                                                            agcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCAcc                                                                                    >  3:151291/1‑132 (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                                         <  8:361527/121‑1 (MQ=39)
                                                                                                                              cagagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcag                                         >  7:361527/5‑125 (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)
                                                                                                                                               atattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGATTAAAt                                                                                              <  7:274708/55‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGG  >  NC_002947/3733124‑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: