Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,733,272	T→A	9.2%	intergenic (‑32/‑40)	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_002947	3,733,272	0	T	A	9.2%	183.9 / 6.9	65	intergenic (‑32/‑40)	PP_3299/PP_3300	lipoprotein/Transcriptional regulator, tetR family
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (3/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGGTA  >  NC_002947/3733130‑3733415                                                                                                                                               |                                                                                                                                                 ggCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTca                                                                                                                                            >  5:200633/1‑149 (MQ=255)   ccAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGACGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCAAAATATAATATGGTcaca                                                                                                                                          <  2:86488/149‑1 (MQ=255)         cAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGt                                                                                                                                    <  1:306474/149‑1 (MQ=255)          aaCCTCTGGCTGTTTATAATCCGAGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGt                                                                                                                                    >  1:63635/1‑148 (MQ=255)           aCCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTc                                                                                                                                   >  4:109262/1‑148 (MQ=255)             ctctGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                                <  8:117573/149‑1 (MQ=255)             ctctGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                                <  3:153505/149‑1 (MQ=255)             ctctGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGCTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                                <  2:101749/149‑1 (MQ=255)                tGGCTGTTTATAATCCGGCCCAACCATGCACCCTGCAGTGAGACTGTCGCTCAACATCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCa                                                                                                                              <  5:343092/148‑1 (MQ=255)                      tttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACActct                                                                                                                                      <  5:251884/134‑1 (MQ=255)                      tttATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACActct                                                                                                                                      >  6:251884/1‑134 (MQ=255)                             tacGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCtt                                                                                                                <  8:112347/147‑1 (MQ=255)                             tCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTAGGTCACACTCTGTCAATCAGTTCCAGGATGCtt                                                                                                                >  3:299780/1‑149 (MQ=255)                               cgtgCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCa                                                                                                              <  5:115232/146‑1 (MQ=255)                                   ccAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCa                                                                                                                        >  5:187361/1‑135 (MQ=255)                                   ccAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCa                                                                                                                        <  6:187361/135‑1 (MQ=255)                                   ccAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCaaaaa                                                                                                          >  4:10650/1‑149 (MQ=255)                                     aaCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAAt                                                                                                         >  6:15492/1‑148 (MQ=255)                                     aaCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAAt                                                                                                         >  6:70650/1‑148 (MQ=255)                                      aCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCGGGAt                                                                                                                    <  2:243032/136‑1 (MQ=255)                                      aCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCGGGAt                                                                                                                    >  1:243032/1‑136 (MQ=255)                                      aCCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGa                                                                                                       <  3:109262/149‑1 (MQ=255)                                        cATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGt                                                                                                     <  5:70650/149‑1 (MQ=255)                                          tGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCTAAAATGAGTa                                                                                                    >  7:125660/1‑148 (MQ=255)                                               ccTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTca                                                                                                                                            <  3:75979/103‑1 (MQ=255)                                               ccTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTca                                                                                                                                            >  4:75979/1‑103 (MQ=255)                                                cTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTa                                                                                                    <  6:221157/142‑1 (MQ=255)                                                cTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTa                                                                                                    >  5:221157/1‑142 (MQ=255)                                                   cAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTcc                                                                                         <  3:6438/149‑1 (MQ=255)                                                   cAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTAGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTcc                                                                                         <  4:299780/149‑1 (MQ=255)                                                                gcTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGgcg                                                                            >  5:239298/1‑149 (MQ=255)                                                                  tCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCGCCACCCTCGCCGgcgcg                                                                          >  6:82364/1‑149 (MQ=255)                                                                    aaCAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGcgcgcg                                                                        <  6:239298/149‑1 (MQ=255)                                                                    aaCAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGcgcgcg                                                                        >  8:160765/1‑149 (MQ=255)                                                                        gcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                                >  7:136008/1‑90 (MQ=255)                                                                        gcagAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAAt                                                                                                                                <  8:136008/90‑1 (MQ=255)                                                                                 tCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCATCGCCGGCGCGCGGGCCAGCCGACCa                                                           >  3:66314/1‑149 (MQ=255)                                                                                  cGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCac                                                          >  5:288498/1‑149 (MQ=255)                                                                                        gTATGGAATCATTATGAATCTCATAGAACAATACGATGTGTGACCATATTATATTTTGGTCAAAATCTGTCAATCATTTCCAGGATGCTTCAAAAATGATTAAATACTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGt                                                     <  1:65763/149‑1 (MQ=255)                                                                                           gggAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAgg                                                 <  3:10650/149‑1 (MQ=255)                                                                                               aGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatc                                             <  6:224363/149‑1 (MQ=255)                                                                                                  atgatgATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcaga                                          <  5:15492/149‑1 (MQ=255)                                                                                                      tgatTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGgatcagatca                                       <  6:200633/148‑1 (MQ=255)                                                                                                       gatTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGa                                                              >  4:24797/1‑124 (MQ=255)                                                                                                       gatTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGatcatc                                     >  6:319761/1‑149 (MQ=255)                                                                                                         ttCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGatcatc                                     >  2:323770/1‑147 (MQ=255)                                                                                                            tcATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGc                                >  7:13590/1‑149 (MQ=255)                                                                                                              aTAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTg                              >  5:251459/1‑149 (MQ=255)                                                                                                              aTAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTg                              >  6:179850/1‑149 (MQ=255)                                                                                                                 gACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCa                           >  8:118811/1‑149 (MQ=255)                                                                                                                     aaGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATTCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGa                       <  6:251459/149‑1 (MQ=255)                                                                                                                       gaggagGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGc                     <  2:63635/149‑1 (MQ=255)                                                                                                                            agtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACACTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCAAGGCTCCCATGCAGCACtac                >  8:258100/2‑147 (MQ=255)                                                                                                                             gtgtGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTtctc                                                                                          >  4:5245/1‑74 (MQ=255)                                                                                                                             gtgtGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTtctc                                                                                          <  3:5245/74‑1 (MQ=255)                                                                                                                             gtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGc                                >  7:172936/1‑132 (MQ=255)                                                                                                                             gtgtGACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGc                                <  8:172936/132‑1 (MQ=255)                                                                                                                              tgtgACCAAAATATAATATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGTTGCCATGGAGCACTTCGc              >  7:62737/1‑149 (MQ=255)                                                                                                                                  aCCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCAt          >  6:202317/1‑149 (MQ=255)                                                                                                                                  aCCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCAt          <  6:288498/149‑1 (MQ=255)                                                                                                                                    cATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTa        <  1:323770/149‑1 (MQ=255)                                                                                                                                      aaatataatATGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATg      <  7:258100/146‑1 (MQ=255)                                                                                                                                        ttatattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCGATTATggg    >  3:140690/1‑149 (MQ=255)                                                                                                                                        ttatattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATggg    >  3:318600/1‑149 (MQ=255)                                                                                                                                           tattTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGGTa  <  7:118811/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 GGCCAGCGCAACCTCTGGCTGTTTATAATCCGGGCCAACCATGCACCCTGCAGTGAGACTGGCGCTCAACAGCAGAGCAATCGTTCGGTAGGGAAGCATGATGATTCTCATAGACCAAGAGGAGGTGTGACCATATTATATTTTGGTCACACTCTGTCAATCAGTTCCAGGATGCTTCAAAAATGAGTAAAT‑CTTCTCCACCCTCGCCGGCGCGCGGGCCAGCCGACCACGACGTGAGGGATCAGATCATCCAGGCTGCCATGGAGCACTTCGCCCATTATGGGTA  >  NC_002947/3733130‑3733415

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