Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	816,792	C→A	20.0%	intergenic (‑145/‑252)	ybhK ← / → moaA	putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	816,792	0	C	A	20.0%	129.3 / 4.8	52	intergenic (‑145/‑252)	ybhK/moaA	putative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (32/9);  new base A (2/9);  total (34/18)
Fisher's exact test for biased strand distribution p-value = 4.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCC  >  NC_000913/816588‑817021                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      tCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGAAAAACCCCCATTTTCCCTTATTCCTTAGc                                                                                                                                                                                                     >  2:411947/1‑239 (MQ=255)  cccgcccagccgcggAGAAAAACGACACGATCAAGATCAGCCAGGGTCAGATGGCGCAAATAGTTTCCTGGAGTCAGATTATCCGCGTTACAGGAGCGCAAATGGCGGGT‑AAACAGCAATTAA‑CCTG‑CAAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGAGGAGGCACGCTATGCCTCCCTGAGG‑ATGTATATAAAAGAAAAACCGCCATTTTCCCTTATTCTGTa                                                                                                                                                                                                       <  1:501158/227‑1 (MQ=255)      gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGAACGCTATGCCTCCCTGATG‑ATGTATATCAAATTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAt                                                                                                                                                                                                >  2:553777/1‑239 (MQ=255)          tCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAGCCGCCCTTTTCCCTTATTCTGTAGCGAAATAGc                                                                                                                                                                                             >  1:803872/1‑238 (MQ=255)           ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGCTACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTAATTAAAAGAAAAACCGCCATTGTCCCTTATTCTCTAACGAAATAGCAc                                                                                                                                                                                           >  2:329348/1‑239 (MQ=255)           ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCCCGCTATGCTTCCCTGATG‑AGGTATATCAAAGAAAAACCGCACTTTTCCCTTATTCTGAAGGGaaaaaacca                                                                                                                                                                                           >  2:353492/1‑233 (MQ=255)           ccgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAAAGCAc                                                                                                                                                                                           >  1:648444/1‑239 (MQ=255)           cccgcgAGACCAACGCCCAGAACAAGACCATCCGGCGTACGTTGGCGCATAAAGTTCCCTGAAGAAAAATTACCCGCGCTACAGTAGCGAAAATGGCGGGG‑AAACAGAAATTAC‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                                                                                                                                                                                           <  1:440096/235‑1 (MQ=255)                    acaacaccacGAGCAAGATCAGCCGGCGTACGATTGCGAATATATTTTCCGGGATCAAGTTTACCCGCGCTACGGTAGCGAAAATTGCGGGA‑AAACAGCAATTAC‑CCTG‑CAAATATTGAGGATACAAACTATTTTCT‑TCTCCCTCAAAAGGACGATGCCCGCTATGCCTCCCTGATG‑AGGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGa                                                                                                                                                                                   <  1:911994/231‑1 (MQ=255)                     gttaaaagaaacagaGATCAGCCAGCTGACGATGGCGCATATATTTCCCGGAAGTAAGATTATCCGCGCTACAGGGGCGAAAGTGGCGGGA‑AAACAAAAATTAA‑CCTG‑CCAATATTGAGGATACAAACTTTTTTCT‑TCTACCTATAAAGGACGAGGCACGGTATCCCTCCTTGATG‑ATGTATATAAAAGAAAAACCGCCATGTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACg                                                                                                                                                                                 <  1:2321/225‑1 (MQ=255)                       aCGACACGATCAAGACCGCCCGGCTTACTATTGCGCATATTGTTCCCGGGAGTCAGATTACCCGCGCAACGGTAGCCAAAATGGCGGGA‑AAACAGAAATTAC‑CCTG‑CCAATTTTGAGGATACAAACTATTTTCT‑TCTCCCTCAAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATAAAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCt                                                                                                                                                                               <  1:289735/239‑1 (MQ=255)                               aTCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCGAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTTTAGCGAAAAAGCACGAtaaaaacgcaaaaaaaag                                                                                                                                                                       >  2:356061/1‑222 (MQ=255)                                          ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGc                                                                                                                                                            >  1:662084/1‑239 (MQ=255)                                            agcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGa                                                                                                                                                          >  1:212273/1‑239 (MQ=255)                                              cGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGCTG‑ATGTATTTCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATAATAACCCtaaaaaaattctgtcttttataca                                                                                                                                                        >  2:2321/1‑217 (MQ=255)                                              cGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACTTAGCGaaa                                                                                                                                                        >  1:240843/1‑239 (MQ=255)                                                    atggcgcATATAGTTTCCTGGAGTCAGATTATCCGAGCTACAGTAGCGCAAATTGCGCGG‑AAACAGCAATTAA‑ACTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCCGATG‑TTGTCTATAAAAGTAAAACCGCCCTTTTCCCTTATTCTGAAGCGAACTAGGACGAGCATGACGCAAAATACATAAttcaaaaacgaaatggggg                                                                                                                                                  >  2:621540/4‑222 (MQ=255)                                                        cgcATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAgtgggattgg                                                                                                                                              >  1:365634/1‑232 (MQ=255)                                                        cgcATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGAAAAACCGCCATTTTCCCTTATTCTGTGACGAAATAGCCCGATCATGACGCTATATACATGATTACATATCGAAAGCGTGGATgg                                                                                                                                              >  2:8165/1‑239 (MQ=255)                                                                  ttCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTAACTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGAAGCCATATACATGATTACATACCGAAAGTGGGGATGGTaaaaaaccc                                                                                                                                    >  2:67597/1‑235 (MQ=255)                                                                       ggAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATTATTACATAAGGAAAGTGTAGATGGTAAAAAACTAATTTa                                                                                                                               >  2:341389/1‑239 (MQ=255)                                                                            cataTTACCCGCATTACAGGACCAAAAATGGGGGGA‑AAACAGAAATTAA‑CCGG‑CAAATATTGAGAATAAAACCTATTTTCT‑TCTACCTAAAAAGGACGATGCCCGTTATGCCTCCGTGATG‑ATGTATATAAAAATAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGc                                                                                                                          <  1:356061/236‑1 (MQ=255)                                                                             aGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATAACGCAAAATACATGATTACATAGCGCAATGTTGGATGGGAAAAATCTCATTTACACGca                                                                                                                         >  1:285764/1‑238 (MQ=255)                                                                                 ttacccccattacatgggaaaaaTTGGGGGGCAAACA‑AAATTAACCCTG‑CAAATTGTGAGGAACAAAACTATTTTCT‑TCCCCCACAAAAGGAAGATCCATTCTATGCCCCCCTGAGT‑ATTAATATAAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGt                                                                                                                      <  2:212273/219‑1 (MQ=255)                                                                                   tCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTCTATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTCCACGCTAGTATc                                                                                                                   >  1:345504/1‑239 (MQ=255)                                                                                        cctaaatgcacaaaaaTTGCGGAG‑AAACAACATATAA‑CCGG‑ACAATTTAGAGGAAACAAATTTTTCTCT‑TACACCTCAAAAGGACGATGCACTCTATGCCTCTCGTAAG‑ATGTATATAAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCGCGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCAt                                                                                                              <  2:146264/226‑1 (MQ=255)                                                                                         cTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATa                                                                                                             >  1:219449/1‑239 (MQ=255)                                                                                          taaagtaccaaaaaTGGGGGGA‑AAACACCAATTAA‑CCTG‑CCATATTTGAGAAAACAAACTTTTTTCT‑TCTACCTCAAAAGGACGATGCCCGTTATGCCTCCCTGATG‑ATTTATATCAAGGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATaa                                                                                                            <  1:736341/228‑1 (MQ=255)                                                                                           acaggaccaaaaaTGGCGGGA‑AAACAGAAATTAC‑CCTG‑ACAATTGTGAGTATAAAAACTATTTTCT‑TTCACCTAAAAAGGACGATGCACGTTATGCCTCCTTGTGG‑TTTTAAATAAAATAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAAc                                                                                                           <  1:353492/229‑1 (MQ=255)                                                                                             agtggccaatatgggggga‑AAACAGAAATTAC‑CGTG‑CAAATATTGAGGATACAACCTATTTTCT‑CCTCCCTTAAAAGGACGATGCACCCTATGCCTCCCTGATG‑ATGTATATCAAAATAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCa                                                                                                         <  1:150914/224‑1 (MQ=255)                                                                                              agggcccgaaTTGGGGGGAAAACCACAATTAC‑CCTG‑CCAAATTGGAAAAAACAAACTATTTTCT‑TCAACCACAAAGTGGCGATCCACGCTACGCCTCCTGGATG‑ATTTAAAACAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCa                                                                                                         <  1:659655/231‑1 (MQ=255)                                                                                                  cgcAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAACCTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATTTCATTTACACGCTAGTATCGGCATAACAACTaa                                                                                                     >  1:582676/1‑238 (MQ=255)                                                                                                   caaaaattcgtggaaaaCCGAAATCAC‑CCCG‑CCAATTTGGAGCATAAAAAC‑ATTTTCTCTATACTTTAAAGGGACGGTCACCGTTATGCCCCCTTGATA‑ATGAATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTaa                                                                                                     <  1:341389/224‑1 (MQ=255)                                                                                                       aTTGGGGGG‑AACCACAATTTAA‑CCCGCCCAATATGGAGGTTAAAACTTTTTTTCT‑CATCCCACAAAAGGACAACCCGCGATATCCCCCCTTGAAG‑ATGTATAGCAAAGAAAACCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACCct                                                                                                <  2:582676/239‑1 (MQ=255)                                                                                                                    aagatatctacccTG‑CCAATATGGGGAAAACAACTTTTTTTCT‑TCTACCCAAAAAGGACGATGCACGCTATG‑CTACCCAATGTATTAAAAAAAAAATAAAGCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCa                                                                                    <  2:240843/228‑1 (MQ=255)                                                                                                                           tAA‑CCTG‑CAAATATTGGAGATACAAACTATTTTCT‑TCGACCTCAAAGGGACGTTGCACGATAGGCCCCCCGAAGA‑ATGAATAAAAAAGAAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTc                                                                            <  1:67597/239‑1 (MQ=255)                                                                                                                            tacccTG‑CAAAAATTGAGATAACAAACTATTTTCT‑TCTACCTCAAAAGGACGATGCAGTGTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTc                                                                            <  1:411947/236‑1 (MQ=255)                                                                                                                                cTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGAATCGCCATAACCACTAAACAATCTAGCCTCTGCACCGGGGTCAACt                                                                        >  1:120477/1‑239 (MQ=255)                                                                                                                                                  aCAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTaaa                                                                                                    >  1:711427/1‑194 (MQ=255)                                                                                                                                                  aCAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGcc                                                       >  1:777248/1‑239 (MQ=255)                                                                                                                                                  aCAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGGGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCCTAACCACTAAACACCCAAACCTCTTCACCTGGGGCAAATGTTACGGTGCTTTGGcc                                                       >  2:377427/1‑239 (MQ=255)                                                                                                                                                    aaaaTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGCGTGGATGGTAAAAATCTCTTTTACACGCTAGTATCGGCATAACCACTAAACACCCTAGCCTCTGCACCTGGGTCAACAGATAAGGTACttgacccgt                                                     >  2:273815/1‑232 (MQ=255)                                                                                                                                                         aTTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCCCTGCACCTGGGTCAACTGATACGGTGCTTTGTCCGTGACaa                                                >  1:770960/1‑239 (MQ=255)                                                                                                                                                          ttttCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTACTATCGGCAAAACCACTAAAAACCCTACCCTCTGCACCAGGGTCAACGGTTACGGGGCTTTAGCCGTGaaaag                                               >  2:750464/1‑238 (MQ=255)                                                                                                                                                            tttatactgcgtctAAAGGGCTACGCTCGCTATGCCTCCCTGATG‑ATGTATATCAAATTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCGCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATg                                              <  1:377427/228‑1 (MQ=255)                                                                                                                                                                   tACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCCCGTaa                                       >  1:341156/1‑239 (MQ=255)                                                                                                                                                                     ccTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGCAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAg                                     >  1:795395/1‑239 (MQ=255)                                                                                                                                                                          aaaGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACAATCTAGCATCTTCACCTGGGTCAACTGATACGGGACTTTTGCCGTCAAAAAGCCCGAAAAAATTGc                                >  2:710233/1‑239 (MQ=255)                                                                                                                                                                                 gATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCATCTGCACCTGGGTCCACTGATACGGTGCTTTGGCCGTGCCAATGCTCGTAAAAATTGcccccaag                         >  1:224537/1‑237 (MQ=255)                                                                                                                                                                                   tGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGc                                                                                           >  1:13328/1‑171 (MQ=255)                                                                                                                                                                                                  ccTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATg          >  1:524471/1‑237 (MQ=255)                                                                                                                                                                                                          atgTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAATTCTCATTTACACGCTAGGATCGGCATAACCACTAAACAATCCAGCCTCTGCACCTGGGTCACCTGATACGGGGCTTTGCCCGTGAAAAGGCTCGTAAAGATTGCCACCAGGGAGAAGGAAAAAAGAACTTCGcc  >  2:393510/1‑238 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      TCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGG‑AAACAGCAATTAA‑CCTG‑CCAATATTGAGGATACAAACTATTTTCT‑TCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATG‑ATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCC  >  NC_000913/816588‑817021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑