Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,733,263 T→A 10.2% intergenic (‑23/‑49) 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,2630TA10.2% 165.6 / 8.8 62intergenic (‑23/‑49)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 = 5.97e-01

CGCGACGTAATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTA  >  NC_002947/3733115‑3733409
                                                                                                                                                    |                                                                                                                                                  
cgcgACTTAATCATCGGCCAGCGCAACCTCTGGCTGTTAATAATCCGGGCCAACCATGCACCCTGCAGTGATACTGGCGCTAAAAAGCATAGAAATCGTTCGGTAGGGAAGCATTATGATTCTCATAGACCAAGAGGAGGTGTGACCat                                                                                                                                                    <  1:58095/149‑1 (MQ=255)
   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                                                                                                                <  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                                                                                >  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                                   >  4:172892/3‑125 (MQ=255)
                                                                                                                                         gagtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcagat                                   <  3:172892/123‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
CGCGACGTAATGATCGGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAATCTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTA  >  NC_002947/3733115‑3733409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: