Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	782,930	A→G	A282A (GCA→GCG)	nadA →	quinolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	782,930	0	A	G	98.6%	238.4 / ‑4.8	71	A282A (GCA→GCG)	nadA	quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (33/37);  total (33/38)
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.96e-01

ATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCT  >  NC_000913/782797‑783064                                                                                                                                      |                                                                                                                                       aTGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAAcc                                                                                                                                  <  1:1400173/140‑1 (MQ=255)  tGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCg                                                                                                                                 <  1:941482/140‑1 (MQ=255)       tcggtTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGt                                                                                                                            >  2:519914/1‑140 (MQ=255)        cggtTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTg                                                                                                                           <  2:311876/140‑1 (MQ=255)         ggtTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGTGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGa                                                                                                                          <  2:122359/140‑1 (MQ=255)            tCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAggg                                                                                                                       >  2:247258/1‑140 (MQ=255)                 caGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCaa                                                                                                                  <  1:940966/140‑1 (MQ=255)                   gTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCa                                                                                                                   <  2:1107105/137‑1 (MQ=255)                   gTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCa                                                                                                                   >  1:1107105/1‑137 (MQ=255)                   gTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAAcc                                                                                                                >  2:917387/1‑140 (MQ=255)                   gTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAAcc                                                                                                                >  2:482103/1‑140 (MQ=255)                      aaCTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGc                                                                                                             >  2:441270/1‑140 (MQ=255)                              gctgctGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTg                                                                                                     <  1:1016649/140‑1 (MQ=255)                                  ctgcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTgcgcg                                                                                                 >  2:709951/1‑140 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  2:50565/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  1:1229605/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  1:917387/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  1:788946/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  2:588759/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCa                                                                                               <  2:552778/140‑1 (MQ=255)                                    gcGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTCCAACCTGCCGCAGCTGCGCGCa                                                                                               <  2:444300/140‑1 (MQ=255)                                            aTTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGt                                                                                       >  1:1291630/1‑140 (MQ=255)                                              tGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTgg                                                                                     <  1:1189204/140‑1 (MQ=255)                                               gCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGa                                                                                    <  1:666419/140‑1 (MQ=255)                                               gCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGa                                                                                    <  1:1321242/140‑1 (MQ=255)                                               gCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGa                                                                                    <  2:129316/140‑1 (MQ=255)                                                 cacaTCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATg                                                                                  <  1:519914/140‑1 (MQ=255)                                                              tATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGcc                                                                     >  2:1383186/1‑140 (MQ=255)                                                                ttGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCtt                                                                   >  1:818635/1‑140 (MQ=255)                                                                ttGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCtt                                                                   <  1:804286/140‑1 (MQ=255)                                                                   tgGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAg                                                                >  1:734810/1‑140 (MQ=255)                                                                     tcAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAGGGGTTACTGGAAGCGCCAACCGCAGGTGAGGGGGCAACCTGCCGGAGGTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGc                                                              >  1:939811/2‑140 (MQ=255)                                                                     gCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGc                                                              >  2:672874/1‑140 (MQ=255)                                                                        aCCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCAt                                                           >  1:313800/1‑140 (MQ=255)                                                                         ccGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATc                                                          >  1:1098793/1‑140 (MQ=255)                                                                         ccGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGGTACTGGAAGCGCCAACCGCAGGTGAGGGCGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAAGCCATc                                                          >  2:686348/1‑140 (MQ=255)                                                                          cGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCg                                                         >  2:779415/1‑140 (MQ=255)                                                                             tCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCag                                                      >  2:674057/1‑140 (MQ=255)                                                                              cGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCaga                                                     <  1:441270/140‑1 (MQ=255)                                                                              cGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCaga                                                     <  2:673123/140‑1 (MQ=255)                                                                                   tATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCAt                                                <  1:247258/140‑1 (MQ=255)                                                                                      tttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAg                                             >  1:736093/1‑140 (MQ=255)                                                                                        tCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGaa                                           <  2:182708/140‑1 (MQ=255)                                                                                         cTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAAc                                          <  1:482103/140‑1 (MQ=255)                                                                                         cTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAAc                                          <  1:674057/140‑1 (MQ=255)                                                                                         cTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAAc                                          <  1:1266933/140‑1 (MQ=255)                                                                                           aCAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACag                                        >  2:381192/1‑140 (MQ=255)                                                                                               aaTGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaag                                    >  2:1164006/1‑140 (MQ=255)                                                                                                         gCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCa                          <  2:456683/140‑1 (MQ=255)                                                                                                         gCGGTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCa                          >  2:1359592/1‑140 (MQ=255)                                                                                                           ggTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACg                        <  2:1098793/140‑1 (MQ=255)                                                                                                           ggTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACg                        <  2:818635/140‑1 (MQ=255)                                                                                                           ggTGCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCAGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACg                        >  2:1259913/1‑140 (MQ=255)                                                                                                              gCCAGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAgg                     <  1:381192/140‑1 (MQ=255)                                                                                                               ccAGATAAAGGGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGt                    <  1:779415/140‑1 (MQ=255)                                                                                                                 aGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTc                  <  1:709951/140‑1 (MQ=255)                                                                                                                 aGATAAAGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTc                  <  1:937714/140‑1 (MQ=255)                                                                                                                     aaaGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGt              <  2:126889/140‑1 (MQ=255)                                                                                                                     aaaGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGt              <  1:965852/140‑1 (MQ=255)                                                                                                                     aaaGAGTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGt              <  1:812315/140‑1 (MQ=255)                                                                                                                       agagTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtg            >  2:79779/1‑140 (MQ=255)                                                                                                                       agagTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtg            >  1:172257/1‑140 (MQ=255)                                                                                                                       agagTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtg            >  1:134060/1‑140 (MQ=255)                                                                                                                          gTTACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTtgatg         >  2:1128274/1‑140 (MQ=255)                                                                                                                            tACTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCAACGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaa       >  2:367072/1‑140 (MQ=255)                                                                                                                             aCTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      >  2:997639/1‑140 (MQ=255)                                                                                                                             aCTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      >  2:45757/1‑140 (MQ=255)                                                                                                                             aCTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      >  2:43985/1‑140 (MQ=255)                                                                                                                             aCTGGAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      >  2:304863/1‑140 (MQ=255)                                                                                                                                 gAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      >  2:1053370/1‑136 (MQ=255)                                                                                                                                 gAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGaaa      <  1:1053370/136‑1 (MQ=255)                                                                                                                                 gAAGCGCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCt  >  2:1335715/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       ATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCT  >  NC_000913/782797‑783064

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑