Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,024,816	A→C	16.6%	V3G (GTG→GGG)	dsrB ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,024,816	0	A	C	16.6%	107.4 / 2.5	53	V3G (GTG→GGG)	dsrB	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base C (5/5);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCAAAAGAGTACAGCGATACCGCAAAACGTGAGG  >  NC_000913/2024596‑2025028                                                                                                                                                                                                                              |                                                                                                                                                                                                                     ttCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTc                                                                                                                                                                                                                 >  1:100901/1‑226 (MQ=255)          gtggtTTCACACCAGTCAGGCTTAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTCATGCCCTGCTTCATTAAAGAACCAGATGCCGATCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGCACCACGCCAGGCCGCCCCGGAACGCCAACCGATTTGACGGTTCACCCATCACTCCACTTACTTGTTTtcct                                                                                                                                                                                                      >  2:601422/1‑228 (MQ=255)            ggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAAATGCCAGGACCACGCCAGGACGAGCCGGACCACCATCCGGTTTTACTGTTACCCGATCATTCACCTTCATCGTTTtcctcc                                                                                                                                                                                                    >  2:518558/1‑228 (MQ=255)                tCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCGGCGACACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGAACCACGCCAGGACCACGCGGACCCCCATCCGCATTGACTCTTACCCAATCATTCACCCTCAACCGTTtcctccctgg                                                                                                                                                                                                >  2:338087/1‑224 (MQ=255)                tCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCCCCTCCATCGTTTTCCTCCtgtg                                                                                                                                                                                                >  1:526812/1‑228 (MQ=255)                        gTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCACATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGAACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACCCGGACCCCCATCCGTGTTGGCTGTTACCCGATCATTCCCCTTCACCGTTTTCCTCCTGTGGCtttttt                                                                                                                                                                                        >  2:381766/1‑225 (MQ=255)                                  aTAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCCGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCTTTCACCTTCATCGTTTTCCTCCTGTGGCTTTTTGCCAggtcata                                                                                                                                                                              >  2:26118/1‑222 (MQ=255)                                           ctctGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTc                                                                                                                                                                      >  1:512081/1‑227 (MQ=255)                                              tGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAACCTCCTCAACTGCCAGTACCACGCCAGGCCGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGATTTCCCACTGTGGCTTTGTGCCAGCGGTGAACAATTTCGtt                                                                                                                                                                   >  1:424266/1‑227 (MQ=255)                                                cTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGtttt                                                                                                                                                                 >  1:328251/1‑227 (MQ=255)                                                cTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGCCTTTGTGCCAGGGTAGACAAATTTCGttttt                                                                                                                                                                >  2:152028/1‑228 (MQ=255)                                                cTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCCCCTGTGGCTTTTTGCCCGTGTAGAACAATTTCGttttt                                                                                                                                                                >  2:3946/1‑228 (MQ=255)                                                   tctcCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGCTCTGACTGTTACCCGCTCATGCCCCTTAATCGATTTCCTCCTGTGGCTTTGTGCCCGTGTAGAACAATTTCGTTTTTCt                                                                                                                                                              >  1:308487/1‑227 (MQ=255)                                                     tcCACAAAGAATCCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAACCTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGc                                                                                                                                                           <  1:514382/228‑1 (MQ=255)                                                                    tCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGAAGAACAATTTCGTTTTTCTGGCAGCGCCAGccgcccg                                                                                                                                            >  1:16950/1‑228 (MQ=255)                                                                              cTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGCAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCCCCTTCATCGTTTTCCTCCTGTGGCTTTTGTCCAGTTTAGAACACTTTCGTTTTTCTGGCAACCCCCGgcccgccagaggcg                                                                                                                                     >  2:548717/1‑217 (MQ=255)                                                                                 cTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGGACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGGTTTCCTCCTGTGGCTTTGTGCCAGTGTCGAACAATTTCGGTTTTCTGGCAGCGCCGGGcgcgcgcgtgggcgagt                                                                                                                                 >  1:421687/1‑217 (MQ=255)                                                                                    cATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTtct                                                                                                                            >  1:447088/1‑228 (MQ=255)                                                                                      ttAAAGAACCAGATGCCGAGCGGGTTGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCg                                                                                                                          <  1:197296/228‑1 (MQ=255)                                                                                      ttAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGGCGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCAGCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCT‑GCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGa                                                                                                                         <  1:158975/228‑1 (MQ=255)                                                                                      tgaaaGAACAAGATGCGGAGCGGGGAATCTTCCAACGAAGCCAGGTACAGTGGTC‑CTCCTCTAACCTCCTCACTTGCCAGTACCACGCCAGGACGCCGCGGACCGCCATCCGTTTTGGCTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTcgcg                                                                                                                          <  2:141934/226‑1 (MQ=255)                                                                                        aaaGCACCAGATGCCTGGCGGGTAGCCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGAc                                                                                                                        <  1:439439/228‑1 (MQ=255)                                                                                        aaaGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTCCATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGa                                                                                                                         <  1:186255/227‑1 (MQ=255)                                                                                           gAACCAGAGGCCGAGCGGGGAGTCTTCCAGCGAAACCAGGTACATTGGGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGt                                                                                                                     <  1:339703/228‑1 (MQ=255)                                                                                           aaaCCAGGTGCCGCGCGGGTATTCTTCCAGCCAAACCAGGTCCATTGGCC‑CTTCCCTAAACTCCTCAACTGCCCGTCCCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGt                                                                                                                     <  2:16950/227‑1 (MQ=255)                                                                                           aaaCCAGCTGCCGGGCGGGTAGTCTTCCAGCGAAACCAGGTACATGTGGC‑CTTCACAAAACTCCCCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGt                                                                                                                     <  2:328251/227‑1 (MQ=255)                                                                                              cccgATGCCGAGGGGGTTGTCTTCCCGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCta                                                                                                                  <  2:264963/225‑1 (MQ=255)                                                                                              ccAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTCTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCATGCGCGCGCGGGTGCTGATTTTCTAGCAGGTCta                                                                                                                  >  1:595888/1‑228 (MQ=255)                                                                                              ccAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACCCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCta                                                                                                                  <  1:437276/228‑1 (MQ=255)                                                                                                            ggTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCTCCCGGCGCGCGCGAGTGCTGACTTTCTCGTCGGTTTACACtataaagtgtc                                                                                                    >  2:352781/1‑218 (MQ=255)                                                                                                                 ttttCCAGCAAAACCAGGGCCATT‑TGCGCTTCCCTAAACCCCCCAACGGCCCGTCCCCAGCCAGGACGCCGCGGACCGCCATCCGGTTTGACTGTTTCCCGATCATTCCCCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGc                                                                                               <  2:447088/226‑1 (MQ=255)                                                                                                                                      aTTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTtgtg                                                                                                                                                                                        >  1:282269/1‑118 (MQ=255)                                                                                                                                       ttGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGc                                                                                                                                                                                       <  2:282269/118‑1 (MQ=255)                                                                                                                                       ttGTCC‑CTCCACTAACCGCCCCAGCTGCCCGTACCACGCCCGGACGACGCGGACTGTCATTCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCAGGCAGCGCCAGGCGCTCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGcaa                                                                          <  2:45771/227‑1 (MQ=255)                                                                                                                                            c‑cTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGCGTGCTGATTTTCTCGACGGTCTATAATTAAGAGATGCCAGGCGGCCTTAACGACTGGGGGCAACAACag                                                                    >  1:563924/1‑228 (MQ=255)                                                                                                                                                     aaaCTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATCCTTAAGAGAGGCCAGCCGGCCTCAACGACTGGGGGCAACACCAGGGTacagg                                                            >  2:107663/1‑224 (MQ=255)                                                                                                                                                          ccccaccTGCTAGTACCCCGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTGCCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCg                                                       <  2:182739/222‑1 (MQ=255)                                                                                                                                                            tcAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGGGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACTACTGGCGGCAACACCAGAGTAACGGTTGCGAg                                                     >  1:182739/1‑228 (MQ=255)                                                                                                                                                                 cgCCAGTCCCACGGCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTTCCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAg                                                <  1:370235/227‑1 (MQ=255)                                                                                                                                                                         ccaagccaagACGACTCTGACCGCGATCCGTTTTGACTGTTACCCGATCATTCCCCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGaa                                        <  1:352781/228‑1 (MQ=255)                                                                                                                                                                                 ggACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGGGGAAAGATGATGAAAACCCCaa                                >  1:115970/1‑228 (MQ=255)                                                                                                                                                                                   gcggcgcgGACCGCCATCCGTTTTGACTGTTACCCTATCATTCACCTTCATCGTTTTCCGCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCaaaa                              <  2:115970/224‑1 (MQ=255)                                                                                                                                                                                   cacgcgcgGACCGCCATCCGTTTTGACGGTTACCCGATCATTCACCTTCATCGTTTTCCGCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGAAGCGCCAGGCGCGCGCGAGTGCTGATTTGCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCaaaa                              <  1:601422/224‑1 (MQ=255)                                                                                                                                                                                       cgcgGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCa                                                                                                                                                                                     >  1:66043/1‑73 (MQ=255)                                                                                                                                                                                       cgcgGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACACCACAGGAACGGTTGCGAGGAAAGATGATGAAACCACCaaaaagga                          >  1:424728/1‑224 (MQ=255)                                                                                                                                                                                       cgcgGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAACGATGCTGAAAACCGCAAAagagg                          >  1:270732/1‑227 (MQ=255)                                                                                                                                                                                        gcgGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAg                                                                                                                                                                                    <  2:66043/73‑1 (MQ=255)                                                                                                                                                                                              cGCCATCCGTTTTGACTGGTACCCGATCTTTCCCCTTCATCGTTTTCCTCCTGTGGCTTTGGGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCAAAAGAGTACAg                      <  1:26118/225‑1 (MQ=255)                                                                                                                                                                                                 cATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGCACAATTTCGTTTTTCTGGCAGCGCAAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGCGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTCACGGCTGCTGGGAAAGCTGATGAAAACCGCAAAAGAGTAAAGCGATAc                >  1:552968/1‑228 (MQ=255)                                                                                                                                                                                                           ggcTGTTACCCGATCATTCCCCTTCATCGTTGTCCCCCGGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCAAAAGAGTACAGCGATACCGCAAAACGt      <  2:512081/226‑1 (MQ=255)                                                                                                                                                                                                           gACTTTTACCCGATCATTCCCCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCAAAAGAGTACAGCGATACCGCAAAACGt      <  1:3946/228‑1 (MQ=255)                                                                                                                                                                                                              tGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGGAAGATTATGAAAACCCCAAAAAAGAAAAACGGAACCCCAAAAccggag   >  2:76538/1‑223 (MQ=255)                                                                                                                                                                                                               gTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCCGCCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAAATGATGAAAACGGCAAAAGGGTAACGCGATACCGCAAAACGTCAgg  >  1:109552/1‑228 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                     TTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGC‑CTTCACTAAACTCCTCAACTGCCAGTACCACGCCAGGACGACGCGGACCGCCATCCGTTTTGACTGTTACCCGATCATTCACCTTCATCGTTTTCCTCCTGTGGCTTTGTGCCAGTGTAGAACAATTTCGTTTTTCTGGCAGCGCCAGGCGCGCGCGAGTGCTGATTTTCTCGACGGTCTATACTTAAGAGATGCCAGGCGGACTTAACGACTGGCGGCAACAACAGAGTAACGGTTGCGAGGAAAGATGATGAAAACCGCAAAAGAGTACAGCGATACCGCAAAACGTGAGG  >  NC_000913/2024596‑2025028

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