Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,616	G→A	L52F (CTT→TTT)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,616	0	G	A	93.7%	174.0 / 2.2	65	L52F (CTT→TTT)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (21/38);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTG  >  NC_000913/2823374‑2823860                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ggCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGatat                                                                                                                                                                                                                                                <  1:37632/251‑1 (MQ=255)                         cTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGAGGTTTTACCGGACGATTCCGGTCCGTAGGTTTCGACGATACGGCCCATCGGCAGACCACCTGACCCACACGCGATATACAGTGAAAGCGAGCCGGTagagc                                                                                                                                                                                                                        >  1:661597/1‑250 (MQ=255)                          ccGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGAt                                                                                                                                                                                                                       <  1:67669/250‑1 (MQ=255)                              cGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGtt                                                                                                                                                                                                                   <  1:22899/251‑1 (MQ=255)                                             aGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGaa                                                                                                                                                                                                   <  1:321869/251‑1 (MQ=255)                                              gCACAGCAGGGTGTCGAAATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAAc                                                                                                                                                                                                  <  2:261879/251‑1 (MQ=255)                                               cacaGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACg                                                                                                                                                                                                 >  1:514882/1‑251 (MQ=255)                                               cacaGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACg                                                                                                                                                                                                 <  1:27297/251‑1 (MQ=255)                                                  agcagGTTGTCGATATCGCCGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGTTAAACGCACAGGTTTTACCTTCCCGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTCTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTc                                                                                                                                                                                              <  2:563646/251‑1 (MQ=255)                                                            cgagatCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCAGCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAgg                                                                                                                                                                                    <  1:70855/247‑1 (MQ=255)                                                            cGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCTCAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGTTCTTCACCCAgg                                                                                                                                                                                    <  1:550730/251‑1 (MQ=255)                                                                   gACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCatg                                                                                                                                                                              >  1:346802/1‑250 (MQ=255)                                                                     cGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCatgatt                                                                                                                                                                           <  1:498764/251‑2 (MQ=255)                                                                      gCCCAGTTTACGTGCGGAGATTGGGCCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTGGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCAAGTGAAAGCGCACCGATAAAGATGGTTTCC‑CAGTCCAAGGAACGGTCTTCACCCAGGCGCATGATgg                                                                                                                                                                          >  2:237494/1‑251 (MQ=255)                                                                        ccAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAg                                                                                                                                                                        <  1:258357/251‑1 (MQ=255)                                                                         cAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGc                                                                                                                                                                       <  1:620884/251‑1 (MQ=255)                                                                           gTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCt                                                                                                                                                                     <  1:322426/251‑1 (MQ=255)                                                                             ttACGTGCGTTGAGTGGGTCCAGCGCGGGTTCAGACTCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCGGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCATA‑TCCATGGAACGGTCTGCACCCAGGCGCATGATGGAGCCtt                                                                                                                                                                    <  2:607639/250‑1 (MQ=255)                                                                                 gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTAAAAGCGAACCGGTAGAGATGGTTTACACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCa                                                                                                                                                              <  1:240365/251‑1 (MQ=255)                                                                                 gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACATTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGACCGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑Acc                                                                                                                                                               >  1:588020/1‑250 (MQ=255)                                                                                        gATTTGGTCCAGCTCGTGTTCAGCTCCGATAAACGCACTGGTTTTTCCTTCACGCGGCGCTGCGG‑CGCTCACCTGCAGCGCCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCTAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCTCATGATGGAGCCTTT‑ACCAAATTGtt                                                                                                                                                       <  1:582868/251‑1 (MQ=255)                                                                                         aTTGGGTCCTGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑AGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATAACCAGTGAAAACGAACCGGTAGAGATGGTTTCCACA‑TCCATTGAACGGTCCTCACCCAGGCGCATGATGGCGCCTTT‑ACCaacacgttt                                                                                                                                                      >  2:361063/1‑244 (MQ=255)                                                                                           tGGGTCCAGCCCGTTGTCAGCATCGATAAATGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑GGATAACCTGCAGCGTCAGCGTGGTTTGACCGGAAGATTCCGGGCCGTATATTTCGTCGATACGGCCCATCGGCAGACCCCCTGCCCCAAACGCGATATCCAGTGAAAGCGATCCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTGT‑ACCAAATTGTTtc                                                                                                                                                     <  2:380190/250‑1 (MQ=255)                                                                                                ccAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAAt                                                                                                                                                <  1:669648/250‑1 (MQ=255)                                                                                                ccAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAAt                                                                                                                                                <  1:18931/250‑1 (MQ=255)                                                                                                 cagcgctggtTCAGCACCGATAAACGCACAGGTTTTTCCTTCACGCTGCGCTGCGG‑CGAGTACCTGCAGCGGCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCt                                                                                                                                              <  2:553252/243‑1 (MQ=255)                                                                                                 cAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTCCCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCt                                                                                                                                              >  1:207357/1‑251 (MQ=255)                                                                                                      cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGccc                                                                                                                                         >  1:358648/1‑251 (MQ=255)                                                                                                            cAGCATCGATAAACGCACAGGTTTTACCTTAACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGtgct                                                                                                                                   >  2:122446/1‑251 (MQ=255)                                                                                                            aagcttCGAAAAACGCCCAGGTGTTTCCTCCACGCGGCGCT‑CGGTCGACCACCTGCTGGGGCAGCTTGCTGTGACCTGAAGATTCCTGTCCGTAGATTTTGACGATACGGCGCCTCGGCAGACCACCTGCCCCACACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAGCGGTCTTCACCCTGGCGCATGATGGAGCCTTT‑ACCAAATAGTTTCTCAAGCTGGcttctttct                                                                                                                                   <  2:518527/246‑9 (MQ=255)                                                                                                                aTCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTgccg                                                                                                                               >  1:578345/1‑251 (MQ=255)                                                                                                                     aaaaCGCACAGGTTTTACCTTCCCGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAAc                                                                                                                          <  1:127785/250‑1 (MQ=255)                                                                                                                            acaGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCt                                                                                                                   >  1:238792/1‑251 (MQ=255)                                                                                                                              aGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCTGGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCtgt                                                                                                                 <  2:140853/251‑1 (MQ=255)                                                                                                                              aGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCtgt                                                                                                                 <  1:601495/251‑1 (MQ=255)                                                                                                                                                  gctgcgG‑CAATAACCTGCAGCGTCACCGTGGTTTTACCGGAAGCTTCCGGTCCGTAGCTTTCGTCGATACGGCCCATCGCCAGGCCACCTGCCCCGACCGCGAACTCCAGTGAAAGCGAACCGGAAGAGAAGATTTCCACA‑TCCTTGGAACGGTCTTCACCCAGTCGCGTGATGGGCCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTACCCCCCACGCATTCTGTTTGTGTTCGTCGATATCCAt                                                                                             >  2:352440/1‑251 (MQ=255)                                                                                                                                                   ctgcgT‑CGATCACCTGCAGCGTCAGCGGGGTTTTACCGGAAGATTCCGGTCAGTAGATTTCGACGCTACGGCCCATCGGCAGACCACCTGACCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGATTCCACA‑TCCATGGCACGGTCTTTACCCCGGCGCCGGATGGAGCCTTTCACCAA‑TTGTTTCTCAATATGGCCCCGTTCGGCCGCAACCTCTTTCTGGTTATTTTCGTCGCTTGCCAtt                                                                                            >  2:560142/1‑251 (MQ=255)                                                                                                                                                     gcgg‑tgATCACCTGCAGCGTCAGGTTGGTTTTACCGGAATCTTCCTGCACGTAGATTTCGTTGATACGACCCATCGTCAGACCCCCTGCCCCAAACGCGAGATCCAGTGAAAGCGAACCGGTAGATATTGTTTCCACA‑TCCATTGAACGGTCTTCCCCCAGTCGCATGATGGAGCTTTT‑ACCAAATTGTTTCTCAATCTGTCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCTTCGATAGCCAtttt                                                                                          <  2:199454/251‑1 (MQ=255)                                                                                                                                                     gcgg‑cgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAtttt                                                                                          <  1:563994/251‑1 (MQ=255)                                                                                                                                                       cggtgctcacccgtggcgtccgcgtGGTTTTCCCGGAAGATTCCGTTCCGTAGAGTCCTGCGAGAGGTCCAACCGGCAGACCACCTGCCCCAGCCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑GCCATGGAACGGTCTTCCCCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTTCTGCCTCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAttttt                                                                                         <  2:439721/236‑1 (MQ=255)                                                                                                                                                         gcgacccccTGCAGCGTAAGTGTGGTTTTCCCTGAAGATGCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCGTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTAc                                                                                       <  2:413950/244‑1 (MQ=255)                                                                                                                                                          cgATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAAAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGAATGATGGCGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAtttttcct                                                                                      >  1:271082/1‑248 (MQ=255)                                                                                                                                                           gATCACCTGAAGCGTCAGCGTGGCTTTACCGGAAGGTTCCGGTCCGTAGCTTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGCTATCCAGTGAAAGCGAACCGGACGAGATGGTTTCCACA‑TCAATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCAGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTc                                                                                     >  2:283958/1‑251 (MQ=255)                                                                                                                                                             tCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCt                                                                                   >  1:499001/1‑251 (MQ=255)                                                                                                                                                                            gcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCCCCCAGGCGCATGAGGGAGCCTTT‑ACCAAATTGTGTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAAt                                                                    <  1:283958/251‑1 (MQ=255)                                                                                                                                                                            gcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAAt                                                                    <  1:122446/251‑1 (MQ=255)                                                                                                                                                                               tGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGgc                                                                                                                                                                                   >  1:20667/1‑138 (MQ=255)                                                                                                                                                                                      aCCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAAGTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAg                                                           <  1:394544/250‑1 (MQ=255)                                                                                                                                                                                          gaagtgtCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCTCCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTGGGGTTGTTTGCGCCGATAGCCACTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAat                                                      <  1:392846/245‑1 (MQ=255)                                                                                                                                                                                          gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCCCCTGCCCCAAAAGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCCGTGATGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCgg                                                               >  2:625255/1‑242 (MQ=255)                                                                                                                                                                                          gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGTACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCGGGCCCAGGGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATGGCCAATTTTACTGCTGTCATGCCGGTTATTTCCGTATAGTccat                                                      >  2:491342/1‑248 (MQ=255)                                                                                                                                                                                          gAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCgg                                                               <  1:625255/242‑1 (MQ=255)                                                                                                                                                                                             gagtCCGGTCCGTAGATTGCGACGATACGGACCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCTGGAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATATCCATGTTTACTCCTTTCATGCCGGGTAATACCGGATAGTCAATATg                                                   <  1:560142/248‑1 (MQ=255)                                                                                                                                                                                              aTTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGt                                                  >  1:238596/1‑251 (MQ=255)                                                                                                                                                                                                gccgttccgtcaaTTTCTACTATACGGCCCATCGGCAGAACACTCGCCTCAAACGCGATATCCATCGAAAGCGAACCGGTAGAGATTGTGTCCACA‑TTCATGGAACGGTCTACACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTTTTTCTCAATCAGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTGTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCTGATAGTCAATATGtt                                                 <  2:19745/238‑1 (MQ=255)                                                                                                                                                                                                                cgacgaTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAt                                                                                             <  1:65268/190‑1 (MQ=255)                                                                                                                                                                                                                cgacgaTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAt                                                                                             >  2:65268/1‑190 (MQ=255)                                                                                                                                                                                                                 gacgaTACGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATAc                               <  1:315139/251‑1 (MQ=255)                                                                                                                                                                                                                        cGGCCCATCGGCAGACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGGTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGc                        <  1:162839/251‑1 (MQ=255)                                                                                                                                                                                                                                     gACCACCTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATc            >  1:134386/1‑250 (MQ=255)                                                                                                                                                                                                                                      accaccTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGCTGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGACCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGAAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTAACCTGTATGCTAATACAgaaaca           >  2:643246/1‑245 (MQ=255)                                                                                                                                                                                                                                       ccaccTGCCCCAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTGCTGTTGAAGCAATTATACTGTATTCTCATACAGTATCAAg         <  2:667111/251‑1 (MQ=255)                                                                                                                                                                                                                                             gCCCCAAACGCGATCACCAGTGAAAGCGACCCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgtggtt    <  2:569050/250‑4 (MQ=255)                                                                                                                                                                                                                                             gCCCCAAACGCGATATCCAGTGAAAGCGAACCGGGAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGttt    <  2:174746/250‑1 (MQ=255)                                                                                                                                                                                                                                              ccccAAACGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGGTTCCACA‑TCCATGGAACGGTCTTCACACAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTg  >  1:603796/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG‑CGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACA‑TCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTT‑ACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTG  >  NC_000913/2823374‑2823860

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