Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,442,673	A→C	8.2%	K604T (AAG→ACG)	topA →	DNA topoisomerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,442,673	0	A	C	8.2%	217.7 / 10.2	73	K604T (AAG→ACG)	topA	DNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base A (30/37);  new base C (3/3);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

AACCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGCTGACGACGAGGGC  >  NC_002947/2442525‑2442818                                                                                                                                                     |                                                                                                                                                   aaCCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCaa                                                                                                                                                    <  8:111934/149‑1 (MQ=255)    cTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAgg                                                                                                                                                  >  4:298403/1‑148 (MQ=255)      cgacaaCGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAt                                                                                                                                               >  8:305607/1‑149 (MQ=255)         cgacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACC‑ACATTCCTTGCACGGAAtg‑ttg                                                                                                                                          <  7:14521/149‑3 (MQ=255)          gacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCACGGAAtg‑ttg                                                                                                                                          <  7:115327/149‑3 (MQ=255)          gacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑Tgg                                                                                                                                          <  8:299247/149‑1 (MQ=255)            cgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCACGGAAtg‑ttggt                                                                                                                                        <  2:176463/149‑5 (MQ=255)              tCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCgg                                                                                                                                      <  4:792/149‑1 (MQ=255)                  ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAtg‑tg                                                                                                                                           <  1:301002/140‑1 (MQ=255)                  ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAtg‑tg                                                                                                                                           >  2:301002/1‑140 (MQ=255)                  ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCga                                                                                                                                  >  8:268119/1‑149 (MQ=255)                     gggCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgatga                                                                                                                               <  1:374263/149‑1 (MQ=255)                       gCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgatgatg                                                                                                                             <  2:157374/149‑1 (MQ=255)                        cGAGCGTGACTGGAAGAACGTGCGCGACGAGTTCTACGGCGTATTCAGCAAAATGCTGCAGTCTGGCGAGTCCAGCGAAGACGGCATGCTTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgatgatga                                                                                                                            <  6:112094/149‑1 (MQ=255)                          aGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCAAGGAATG‑TGGTCGGCCGATGATGAtt                                                                                                                          <  1:225541/149‑1 (MQ=255)                          aGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGAtt                                                                                                                          <  3:298403/149‑1 (MQ=255)                           gCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTc                                                                                                                         <  3:133884/149‑1 (MQ=255)                              tGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGAtt                                                                                                                          <  5:58100/145‑1 (MQ=255)                              tGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGAtt                                                                                                                          >  6:58100/1‑145 (MQ=255)                                aCTGGAAGAACGTGCTCGACGAGTTCTAGGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCAc                                                                                                                     <  8:353758/148‑1 (MQ=255)                                aCTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGGCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACGAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCAcc                                                                                                                    <  5:121497/149‑1 (MQ=255)                                   ggAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCgg                                                                                                                                      <  8:296278/128‑1 (MQ=255)                                   ggAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCgg                                                                                                                                      >  7:296278/1‑128 (MQ=255)                                    gaagaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCt                                                                                                                 >  1:357071/1‑148 (MQ=255)                                    gaagaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCAGATGATGCTTCGCACCGCtt                                                                                                                >  2:186499/1‑149 (MQ=255)                                       gaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTa                                                                                                             >  3:57793/1‑149 (MQ=255)                                          cGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACgg                                                                                                          <  4:57793/149‑1 (MQ=255)                                          cGTGCTCGACGAGTTATACCGCGATTTCAGCAAAAACCTGCAGACTACCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCCTGCAAGGAATG‑TGGTCGCCCGATGATGATTCGCACAGCTTCTACgg                                                                                                          >  8:114918/1‑149 (MQ=255)                                             gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgat                                                                                                                                 <  3:282120/123‑1 (MQ=255)                                             gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgat                                                                                                                                 >  4:282120/1‑123 (MQ=255)                                             gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgtg                                                                                                       >  4:87784/1‑149 (MQ=255)                                                 gacgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                                        <  1:353547/144‑1 (MQ=255)                                                 gacgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                                        >  2:353547/1‑144 (MQ=255)                                                  acgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTc                                                                                                  >  5:266742/1‑149 (MQ=255)                                                    gaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGc                                                                                                                  >  8:76441/1‑131 (MQ=255)                                                    gaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGc                                                                                                                  <  7:76441/131‑1 (MQ=255)                                                          tACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCt                                                                                                >  6:183392/1‑143 (MQ=255)                                                          tACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCt                                                                                                <  5:183392/143‑1 (MQ=255)                                                              gCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                                      <  6:86324/149‑1 (MQ=255)                                                                gATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                                    <  3:87784/149‑1 (MQ=255)                                                                        cAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                                    <  3:38275/141‑1 (MQ=255)                                                                        cAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                                    >  4:38275/1‑141 (MQ=255)                                                                                  cAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                                      >  6:85405/1‑129 (MQ=255)                                                                                  cAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                                      <  5:85405/129‑1 (MQ=255)                                                                                         ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                                        <  4:127246/104‑1 (MQ=255)                                                                                         ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                                        >  3:127246/1‑104 (MQ=255)                                                                                         ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTg                                                                                          <  5:159992/118‑1 (MQ=255)                                                                                         ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTg                                                                                          >  6:159992/1‑118 (MQ=255)                                                                                              tCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCaa                                                      <  2:357071/149‑1 (MQ=255)                                                                                                cAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAgg                                                    >  2:359118/1‑149 (MQ=255)                                                                                                cAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAgg                                                    <  6:266742/149‑1 (MQ=255)                                                                                                  gCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGc                                                   >  8:273129/1‑148 (MQ=255)                                                                                                      aGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCg                                              <  7:273129/149‑1 (MQ=255)                                                                                                        aCGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCt                                                  <  2:80588/143‑1 (MQ=255)                                                                                                        aCGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCt                                                  >  1:80588/1‑143 (MQ=255)                                                                                                            cATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgatgat                                                                                                                              <  6:76019/63‑1 (MQ=255)                                                                                                            cATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCgatgat                                                                                                                              >  5:76019/1‑63 (MQ=255)                                                                                                                    ccAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGcc                                 <  5:333909/148‑1 (MQ=255)                                                                                                                    ccAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCGGCAAGGCTACCGTCAACCTGGTGCCg                                >  3:54905/1‑149 (MQ=255)                                                                                                                      aaCCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                                      >  7:89144/1‑93 (MQ=255)                                                                                                                      aaCCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                                      <  8:89144/93‑1 (MQ=255)                                                                                                                       aCCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGc                             >  1:186301/1‑149 (MQ=255)                                                                                                                          aGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTCCAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGcgac                          <  7:305607/149‑1 (MQ=255)                                                                                                                           gCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTg                                                                                     <  6:45067/89‑1 (MQ=255)                                                                                                                            ccGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGcgacga                        <  7:268119/149‑1 (MQ=255)                                                                                                                              gagcATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCGTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGTCgaga                      <  7:114918/146‑1 (MQ=255)                                                                                                                                     gaccaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCg                >  2:157665/5‑149 (MQ=255)                                                                                                                                      accaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGc               >  5:107811/4‑149 (MQ=255)                                                                                                                                      accaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGc               >  4:276545/4‑149 (MQ=255)                                                                                                                                            ttCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTAc                                                <  4:88307/109‑1 (MQ=255)                                                                                                                                            ttCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTAc                                                >  3:88307/1‑109 (MQ=255)                                                                                                                                            ttCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGt                                             >  2:218647/1‑112 (MQ=255)                                                                                                                                            ttCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGt                                             <  1:218647/112‑1 (MQ=255)                                                                                                                                                   cAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGCTGACGACGAGGGc  >  5:115290/1‑148 (MQ=255)                                                                                                                                                    aaGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGCTGACGACGAGGGc  >  1:48056/1‑147 (MQ=255)                                                                                                                                                    aaGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGCTGACGACGAGGGc  >  1:156972/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AACCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATG‑TGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGCTGACGACGAGGGC  >  NC_002947/2442525‑2442818

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