Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,441,765	C→T	intergenic (‑22/+89)	hslJ ← / ← ldhA	lipoprotein implicated in Novobiocin resistance/D‑lactate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,441,765	0	C	T	97.6%	133.2 / ‑2.9	42	intergenic (‑22/+89)	hslJ/ldhA	lipoprotein implicated in Novobiocin resistance/D‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/14);  minor base . (0/1);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

GCCTTCACCGCTAAAGCGGTTACACATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCA  >  NC_000913/1441522‑1442010                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       gCCTTCACCGTTAAAGCGGTTACTCATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTGTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCa                                                                                                                                                                                                                                                <  2:335976/251‑1 (MQ=255)   cttaaCCGCTAAAGCGGTTACACATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAACGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGACATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTATTCATTGTGGTTCTAGATTACAGTTTTTG‑AATCAgg                                                                                                                                                                                                                                              >  1:273007/5‑251 (MQ=255)   cTTCACCGCTAAAGCGGTTACACATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGGTTTCTGGCGGCTTTTTACCGCTGGTCACGGGCTTACCGTTTACGGTTTCCTGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTATCGCAACACACGCGGCTACTTTCTTCAGTGTGGTTCTCAATTACAGTTTTTG‑TCTCAgg                                                                                                                                                                                                                                              >  1:272985/1‑251 (MQ=255)    ttCACCGCTAAAGCGGTTACACATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAgg                                                                                                                                                                                                                                              >  2:218453/1‑250 (MQ=255)                     tatacaTGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCGTGCGGATTTTTATCGCTGGTTACGGGCTTACCGTTTTCGCTTTCCAGCACAAAGCGTTGATGCTGTAGCTGTTCTGGCGTTACAGCAATTTTGTCATTACTTACACATCCCGCCAGCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATaga                                                                                                                                                                                                                             <  1:24558/247‑1 (MQ=255)                          aTGCTGCCGGAAATCATCATTTTATCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTCCCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCACTTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTAAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATg                                                                                                                                                                                                                        >  2:536876/1‑250 (MQ=255)                            gctgcCTGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGTATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTGCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTGTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGaa                                                                                                                                                                                                                      <  2:346080/250‑1 (MQ=255)                            gctgcCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGaa                                                                                                                                                                                                                      >  1:398106/1‑250 (MQ=255)                             ctgcCGGAAATCATCATTTTTTCACCAAAGCAGGTTTCTGGAGGACTTTAACCGCTGGTCACGGGCTTACCGTTTACGTTTTCCCGCACAAAGCGACGATGCTGTAGCTATTCTGGCCTAACAGCAATTTTGTCATTACTTACACATCCAGCCATCAGCAGGCTTAGCGCAACAAACGCAGCTAATTTCTTCCTAGTGGTACTCCATTACAGTTTTTG‑ACCAAGGATTATTTTACGAATAGAGGATGaa                                                                                                                                                                                                                      >  1:44204/1‑249 (MQ=255)                                ccGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGCTTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTGCACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTTAGGACTATTTTAAGAATAGAGGATGAAAGGTc                                                                                                                                                                                                                 >  2:83484/1‑251 (MQ=255)                                 cGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGTTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCa                                                                                                                                                                                                                <  1:235323/251‑1 (MQ=255)                                 cGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTc                                                                                                                                                                                                                 >  2:265981/1‑250 (MQ=255)                                    aaaTCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTg                                                                                                                                                                                                             <  2:198841/251‑1 (MQ=255)                                             ttttttCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATc                                                                                                                                                                                                    >  2:192988/1‑251 (MQ=255)                                                             gATTTCTTGCGGATTTTTAGCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCTGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGTCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAAGCAGCTCCCCt                                                                                                                                                                                     <  2:398106/250‑1 (MQ=255)                                                             gATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCt                                                                                                                                                                                     <  2:357725/250‑1 (MQ=255)                                                                    ggcggttatttagctctggtaaatggCTTACCGTTTACGCTTTCCAGCACAATGCGATGATTCTGTAGCTGTTCTGTCGTAACAGCAATTTTGTCTTTACTTACACAGCCCGTCATCAGCAGTCTTATCTCAACAAATGCTGCTACTTTCTTCATGGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAGGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAAt                                                                                                                                                                                <  1:367544/224‑1 (MQ=255)                                                                           ttttATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTATGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTTAATCAGCTCCCCTGGAATGCAGGGGAGc                                                                                                                                                                      >  1:184965/1‑251 (MQ=255)                                                                                         aCGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGCTGATGCTGTAGCTGTTCTGGCTTAACATCAATTTTATCATTACTTCCACATCCCGCCATCAGCAGGCTTATCTCAACAAACGCGGCTACTTTCTTCATTGTTGTTCTCAATTACAGTTTTTG‑ACTCAGTCTTATTTTATGACTCGAGTCTGAAAGGTCATTGGGTATTATCTGAATCCTCTCCCCTGGACTGCAGGCGCGCGTCAAGATTAAAcc                                                                                                                                                        >  2:497066/1‑251 (MQ=255)                                                                                         aCGGGCTTACCGTTTACGATTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAAcc                                                                                                                                                        >  1:179935/1‑251 (MQ=255)                                                                                           gggCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGCTTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAAcc                                                                                                                                                        >  2:329293/1‑249 (MQ=255)                                                                                              cTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACACACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGTATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAgttc                                                                                                                                                   >  2:204323/1‑251 (MQ=255)                                                                                                                  gCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCtt                                                                                                                               >  1:415338/1‑251 (MQ=255)                                                                                                                                   tGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACCGTTTTTG‑ACTCAGGACTATTTTAAGCATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTa                                                                                                                >  1:283409/1‑249 (MQ=255)                                                                                                                                                 cGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGt                                                                                                >  1:109116/1‑251 (MQ=255)                                                                                                                                                    aaCAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTg                                                                                             >  1:409145/1‑251 (MQ=255)                                                                                                                                                           cttttGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGCGGTTCTCAATTACAGTTTTTG‑ACTCAGGCCTATTTTAAGCATAGAGGATGAACGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTACAGATTGCTTAAGTTTTGCAGCGAAGTCTGAGAAAt                                                                                       >  2:117464/2‑250 (MQ=255)                                                                                                                                                           aTTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGGTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATa                                                                                      <  1:417744/251‑1 (MQ=255)                                                                                                                                                             tttGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATAc                                                                                     <  1:525338/250‑1 (MQ=255)                                                                                                                                                               tGTCATTCCTTACACATCCCGCCATCAGCAGGTTTAGCGCACCAAACGCGGATACTTTCTTCCTTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTAATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGATTAAGTTTTGCAGCGTAGTCTGAGCAATACTgg                                                                                  >  2:585916/1‑251 (MQ=255)                                                                                                                                                                  cATTACTTACCCATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCATGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTcc                                                                               >  1:94176/1‑250 (MQ=255)                                                                                                                                                                                                        acaaACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAAc                                         >  2:50245/1‑251 (MQ=255)                                                                                                                                                                                                           aaCGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCCGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTgcat                                                                                                                                                                                               >  1:337581/1‑95 (MQ=255)                                                                                                                                                                                                                 gcttcTTTCATCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCATGACTATTTTAAGGATAGATGATGAAAGGTCATTGTGAATTATCTGAATCTGCTCCCCTTGAATGCGGGGGCGCTGCAAGATTAAACCAGTTCCTTCGGGCAGGTTTAACCGTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGCCTGTTAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGATCCCTGGAAACAGCACGTtg                                <  2:430365/247‑1 (MQ=255)                                                                                                                                                                                                                    actgtcttcATTGTTGTGCGCAATTACAGGTT‑TGTACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCATGGGAGCGGCAAGATTAAACCAGTTCGTGCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCTTAGTCTGAGATATGCTGGTCAGAGCTTCTGCTGTCAGGAATTCCTGGTGCCCGGCAAACCGCACGTtgtg                              <  2:360394/246‑1 (MQ=255)                                                                                                                                                                                                                    aCTTTCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTtgtg                              <  2:409145/250‑1 (MQ=255)                                                                                                                                                                                                                      tttCTTCATTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTggc                            >  1:438586/1‑250 (MQ=255)                                                                                                                                                                                                                             aTTGTGGTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGAc                    >  1:290083/1‑251 (MQ=255)                                                                                                                                                                                                                                   gTTCTCAATTACAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCg               >  1:70034/1‑250 (MQ=255)                                                                                                                                                                                                                                        cAATTACAGTTTTTG‑ACTCAGGACTAGTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGTTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGaa         <  1:559882/251‑1 (MQ=255)                                                                                                                                                                                                                                             aCAGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGGCAGGCGACGGACTACGt    >  1:338274/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGTTTTTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCa  <  2:265901/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GCCTTCACCGCTAAAGCGGTTACACATGCTGCCGGAAATCATCATTTTTTCACCAAAGCTGATTTCTGGCGGATTTTTATCGCTGGTCACGGGCTTACCGTTTACGCTTTCCAGCACAAAGCGATGATGCTGTAGCTGTTCTGGCGTAACAGCAATTTTGTCATTACTTACACATCCCGCCATCAGCAGGCTTAGCGCAACAAACGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTG‑ACTCAGGACTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCA  >  NC_000913/1441522‑1442010

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