Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,314,735:1	+T	6.5%	coding (937/1209 nt)	aspC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,314,735	1	.	T	6.5%	307.2 / 13.9	77	coding (937/1209 nt)	aspC	aminotransferase
Reads supporting (aligned to +/- strand):  ref base . (36/36);  new base T (2/3);  total (38/39)
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.37e-01

AAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATC  >  NC_002947/4314588‑4314879                                                                                                                                                     |                                                                                                                                                 aaGCTCAGCCTGCTGCATGAAGGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg‑g                                                                                                                                                 <  4:192389/149‑1 (MQ=255)  aGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑Gc                                                                                                                                                <  7:65582/149‑1 (MQ=255)       aGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTTCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                           <  1:89016/149‑1 (MQ=255)       aGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGa                                                                                                                                           <  8:21085/149‑1 (MQ=255)       aGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGa                                                                                                                                           <  1:352824/149‑1 (MQ=255)       aGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATCGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                           <  7:192979/149‑1 (MQ=255)        gCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAg                                                                                                                                          <  7:25368/149‑1 (MQ=255)              tgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGtt                                                                                                                                     >  4:100231/1‑148 (MQ=255)              tgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGtt                                                                                                                                     >  2:146472/1‑148 (MQ=255)                cATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTaa                                                                                                                                   <  3:100231/148‑1 (MQ=255)                 aTGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAAc                                                                                                                                 <  6:27565/149‑1 (MQ=255)                   gAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                               >  7:243071/1‑149 (MQ=255)                   gAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                               >  2:293414/1‑149 (MQ=255)                   gAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACAGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                               >  3:70731/1‑149 (MQ=255)                    aaaGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCCACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                               <  1:281605/148‑1 (MQ=255)                     aaGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCt                                                                                                                              >  7:42068/1‑148 (MQ=255)                      aGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGa                                                                                                                            <  1:392336/149‑1 (MQ=255)                          ggCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACgg                                                                                                                         <  7:160667/148‑1 (MQ=255)                               ccGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACAtt                                                                                                                   >  2:104877/1‑149 (MQ=255)                                cgcgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTc                                                                                                                  >  8:160577/1‑149 (MQ=255)                                 gcgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCt                                                                                                                 >  1:36940/1‑149 (MQ=255)                                   gTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                               <  6:341968/149‑1 (MQ=255)                                   gTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                               <  5:270837/149‑1 (MQ=255)                                      ttttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGt                                                                                                             <  7:209340/148‑1 (MQ=255)                                       tttACTCGATGCCTGGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGGCTTTGCGAGTGCTTAGCTCATCGTCGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACAATCTTTGtt                                                                                                            <  4:70731/148‑1 (MQ=255)                                                      gACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTg                                                                                            >  3:228633/1‑149 (MQ=255)                                                         cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTAc                                                                                                         >  2:241772/1‑133 (MQ=255)                                                         cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTAc                                                                                                         <  1:241772/133‑1 (MQ=255)                                                         cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCgg                                                                                         >  6:89833/1‑149 (MQ=255)                                                             aTTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCATCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                     <  4:10166/149‑1 (MQ=255)                                                             aTTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAACCGCTGACGGCGCATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                     >  2:392658/1‑149 (MQ=255)                                                              ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgtg                                                                                    <  8:42068/149‑1 (MQ=255)                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                     <  6:235223/145‑1 (MQ=255)                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                     >  5:235223/1‑145 (MQ=255)                                                                 tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAg                                                                                  >  1:21085/1‑148 (MQ=255)                                                                        ccGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCCCATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACgg                                                                            >  2:357995/1‑147 (MQ=255)                                                                         cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                         <  7:160577/149‑1 (MQ=255)                                                                           gggCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                       >  6:256916/1‑149 (MQ=255)                                                                                 cGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTc                                                                  <  2:21085/148‑1 (MQ=255)                                                                                   ccGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGc                                                                >  3:286082/1‑148 (MQ=255)                                                                                           cTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAAAGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAg                                                       >  1:255917/1‑149 (MQ=255)                                                                                              tGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAgg                                                      <  3:317796/147‑1 (MQ=255)                                                                                              tGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAgg                                                      <  4:383640/147‑1 (MQ=255)                                                                                                ccAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACt                                                  <  8:243071/149‑1 (MQ=255)                                                                                                 cAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACtt                                                 >  1:275368/1‑149 (MQ=255)                                                                                                   gTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACtttt                                               >  4:147620/1‑149 (MQ=255)                                                                                                    tGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTc                                              <  1:146472/149‑1 (MQ=255)                                                                                                       aTAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCga                                            >  3:253490/1‑147 (MQ=255)                                                                                                       aTAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGc                                            >  5:35080/1‑148 (MQ=255)                                                                                                            tCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATAcc                                      >  6:14612/1‑149 (MQ=255)                                                                                                             cAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCgg                                                                                         >  1:193085/1‑97 (MQ=255)                                                                                                             cAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCgg                                                                                         <  2:193085/97‑1 (MQ=255)                                                                                                               gCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTgg                                   <  5:237204/149‑1 (MQ=255)                                                                                                                cGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGc                                  <  2:36940/149‑1 (MQ=255)                                                                                                                    gTGCGGTAATGGTCGATATTCTCCAGGCACATA‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACt                              >  1:247386/1‑149 (MQ=255)                                                                                                                     tGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACt                              >  4:356904/1‑148 (MQ=255)                                                                                                                      gCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACt                              >  3:277402/1‑147 (MQ=255)                                                                                                                          tAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                        <  5:256916/149‑1 (MQ=255)                                                                                                                           aaTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                       >  4:292092/1‑149 (MQ=255)                                                                                                                            aTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                       <  6:35080/148‑1 (MQ=255)                                                                                                                            aTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                       <  3:147620/148‑1 (MQ=255)                                                                                                                            aTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                       <  3:186184/148‑1 (MQ=255)                                                                                                                            aTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                       <  4:286082/148‑1 (MQ=255)                                                                                                                              ggTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTa                    >  3:359362/1‑149 (MQ=255)                                                                                                                               gTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTAt                   >  5:59297/1‑149 (MQ=255)                                                                                                                                tCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCTCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATc                  >  2:380038/1‑149 (MQ=255)                                                                                                                                tCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATc                  >  2:385947/1‑149 (MQ=255)                                                                                                                                tCGATATTCTCCAGGCACAGG‑GCCTGGAGGTGTTAAACGCTGACGGCACAATCTTTGTTTACGTGCGATGTCCTGCGGCTGGGAGCAACGGCAGGACAATCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCCCTGAGCTGGGTATc                  >  6:279227/1‑149 (MQ=255)                                                                                                                                      ttCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGTTGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGAc             <  4:228633/148‑1 (MQ=255)                                                                                                                                        ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCaa                                          >  6:25741/1‑117 (MQ=255)                                                                                                                                        ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCaa                                          <  5:25741/117‑1 (MQ=255)                                                                                                                                        ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAt          >  5:287762/1‑149 (MQ=255)                                                                                                                                         tcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACt                              >  6:294479/1‑128 (MQ=255)                                                                                                                                         tcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACt                              <  5:294479/128‑1 (MQ=255)                                                                                                                                           cAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTcc       <  5:89833/149‑1 (MQ=255)                                                                                                                                            aGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTccc      >  3:206560/1‑149 (MQ=255)                                                                                                                                            aGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTccc      <  5:14612/149‑1 (MQ=255)                                                                                                                                                acaTG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATc  <  4:277402/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATC  >  NC_002947/4314588‑4314879

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