Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,314,735:1	+T	10.2%	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	10.2%	215.1 / 18.6	59	coding (937/1209 nt)	aspC	aminotransferase
Reads supporting (aligned to +/- strand):  ref base . (21/32);  new base T (3/3);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATT  >  NC_002947/4314587‑4314872                                                                                                                                                      |                                                                                                                                          cAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg                                                                                                                                            <  2:74818/149‑1 (MQ=255)  aaGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg‑g                                                                                                                                          >  7:114594/1‑149 (MQ=255)    gCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑Gcc                                                                                                                                        >  6:34441/1‑149 (MQ=255)      tCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTg                                                                                                                                      >  4:225143/1‑149 (MQ=255)        aGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                    <  4:274977/149‑1 (MQ=255)           ctgctgCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGTCACATG‑GCCTGGAGgt                                                                                                                                 >  5:312580/1‑149 (MQ=255)           ctgctgCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgt                                                                                                                                 <  4:142706/149‑1 (MQ=255)              ctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGtt                                                                                                                              <  4:89526/149‑1 (MQ=255)              ctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAgaatat                                                                                                                              <  4:326968/149‑6 (MQ=255)                   tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACg                                                                                                                         <  3:238067/149‑1 (MQ=255)                   tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACg                                                                                                                         <  5:34441/149‑1 (MQ=255)                     aaaGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCt                                                                                                                       >  3:101403/1‑149 (MQ=255)                     aaaGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCt                                                                                                                       >  3:68588/1‑149 (MQ=255)                        gTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAg                                                                                                                                   <  7:92472/134‑1 (MQ=255)                        gTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAg                                                                                                                                   >  8:92472/1‑134 (MQ=255)                            gCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGGCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGTca                                                                                                                <  5:65080/149‑1 (MQ=255)                                   cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACattatt                                                                                                         <  6:312580/149‑1 (MQ=255)                                    gTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                        <  6:371656/149‑1 (MQ=255)                                    gTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                        <  3:225143/149‑1 (MQ=255)                                       ttttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGTTAATGGTCGATATTCTCCAGGGACATG‑GCCTGGAGGTGTTAAACTCTGACGGCACATTCTTTGt                                                                                                      <  5:275421/148‑1 (MQ=255)                                          tACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg‑g                                                                                                                                          <  7:50547/109‑1 (MQ=255)                                          tACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg‑g                                                                                                                                          >  8:50547/1‑109 (MQ=255)                                               gATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                              >  8:252086/1‑148 (MQ=255)                                                           tGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGc                                                                                 >  1:284971/1‑149 (MQ=255)                                                             cATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAATGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtg                                                                               <  8:114594/149‑1 (MQ=255)                                                              aTTCTCGACGACGGGGGCTTCGCCGAACGCCTTTGCGAGTGGATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtgt                                                                                                                               <  8:219422/100‑1 (MQ=255)                                                              aTTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                              <  7:338684/149‑1 (MQ=255)                                                              aTTCTCGACGACCCGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAAAGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtgt                                                                                                                               >  7:219422/1‑100 (MQ=255)                                                               ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGATGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgtg                                                                             <  4:68588/149‑1 (MQ=255)                                                                 ctcGACGACCGCGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGGGCGGTAATGCTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                            <  8:296764/148‑1 (MQ=255)                                                                  tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCATGGAGGTGTTAAACGCTCACGGCACATTCTTTGTTTACGTGCTATGTCCTGCGGATGTGAg                                                                           >  8:355865/1‑148 (MQ=255)                                                                    gacgacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCaa                                                                        >  4:91165/1‑149 (MQ=255)                                                                     acgacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAAc                                                                       <  7:252086/149‑1 (MQ=255)                                                                          cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                  <  3:101270/149‑1 (MQ=255)                                                                            gggCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGa                                                                >  6:128255/1‑149 (MQ=255)                                                                               cTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAt                                                             >  1:81660/1‑149 (MQ=255)                                                                                ttCGCCGAAGGCATTTGCCAGTGCATAGCTCTGCTAGGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCATCGGCAGGACAtt                                                            <  4:101403/149‑1 (MQ=255)                                                                                   gCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGc                                                         <  5:51659/149‑1 (MQ=255)                                                                                      gAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTg                                                                                                       <  1:297699/100‑1 (MQ=255)                                                                                      gAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTg                                                                                                       >  2:297699/1‑100 (MQ=255)                                                                                             tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAgg                                                                 >  1:200895/1‑131 (MQ=255)                                                                                             tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAgg                                                                 <  2:200895/131‑1 (MQ=255)                                                                                             tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAtt                                                            <  1:320497/136‑1 (MQ=255)                                                                                             tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAtt                                                            >  2:320497/1‑136 (MQ=255)                                                                                               tGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAATCGCCACTGCCAgg                                               <  3:91165/147‑1 (MQ=255)                                                                                                       cATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCg                                      <  2:81660/148‑1 (MQ=255)                                                                                                        aTAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGc                                     >  8:189242/1‑148 (MQ=255)                                                                                                               aGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTg                             <  2:284971/149‑1 (MQ=255)                                                                                                                  gACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGc                           >  7:289645/1‑148 (MQ=255)                                                                                                                   aCGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCa                         <  3:72477/149‑1 (MQ=255)                                                                                                                        cGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGa                     <  5:128255/148‑1 (MQ=255)                                                                                                                           tAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                 <  6:277060/149‑1 (MQ=255)                                                                                                                             aTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                <  6:248155/148‑1 (MQ=255)                                                                                                                                  cGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATcc          >  6:52354/1‑149 (MQ=255)                                                                                                                                   gATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                <  4:52549/142‑1 (MQ=255)                                                                                                                                   gATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                >  3:52549/1‑142 (MQ=255)                                                                                                                                     tatTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGa       <  5:137501/149‑1 (MQ=255)                                                                                                                                        tctcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGa    >  3:178208/1‑149 (MQ=255)                                                                                                                                         ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                  <  2:165817/86‑1 (MQ=255)                                                                                                                                         ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                  >  1:165817/1‑86 (MQ=255)                                                                                                                                         ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAt   >  1:59063/1‑149 (MQ=255)                                                                                                                                          tcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAgg                                               >  6:356232/1‑104 (MQ=255)                                                                                                                                            cAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt  <  7:189242/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                          CAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATT  >  NC_002947/4314587‑4314872

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