Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,194,176	A→G	6.0%	L38P (CTG→CCG)	nudJ ←	bifunctional thiamine pyrimidine pyrophosphate hydrolase and thiamine pyrophosphate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,194,176	0	A	G	6.0%	206.7 / 3.1	66	L38P (CTG→CCG)	nudJ	bifunctional thiamine pyrimidine pyrophosphate hydrolase and thiamine pyrophosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (38/24);  new base G (2/2);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGCCATTCAC  >  NC_000913/1194031‑1194323                                                                                                                                                  |                                                                                                                                                    aaTGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAg                                                                                                                                                    <  8:421584/147‑1 (MQ=255)        aaGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCgg                                                                                                                                            >  6:355565/1‑148 (MQ=255)         agagGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCCCGGGCGGCGGCTTCCACTAAGGTTGCATCGGCTTCCAGATGCCCgg                                                                                                                                            >  7:147655/1‑147 (MQ=255)            ggAAACGCAAAAACGCCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAgg                                                                                                                                        <  2:230998/148‑1 (MQ=255)              aaaCGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAggtt                                                                                                                                      <  2:23401/148‑1 (MQ=255)              aaaCGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAggtt                                                                                                                                      <  3:148657/148‑1 (MQ=255)                 cgcgAAAACGGCGGTTTATCTGGCGCAAACCACTGATGGATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCGCACAGTTCACGGGCGGCGGCTTCCTCTAAGGTTTCATCGGCTTGCAGATGCCCGGCAggttggt                                                                                                                                   >  5:430715/1‑148 (MQ=255)                    aaaaaCGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTc                                                                                                                                 >  4:44832/1‑147 (MQ=255)                              ttttATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCttt                                                                                                                      <  3:384801/148‑1 (MQ=255)                               tttATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAg                                                                                                                                                    >  3:11902/1‑117 (MQ=255)                               tttATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAg                                                                                                                                                    <  4:11902/117‑1 (MQ=255)                               tttATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTa                                                                                                                     >  5:316529/1‑148 (MQ=255)                                             tCCACTGATGCATCCGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGtctc                                                                                                       <  6:172974/148‑1 (MQ=255)                                              cctcTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCGGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAAt                                                                                                             >  8:33378/4‑141 (MQ=255)                                              ccACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCGGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAAt                                                                                                             <  7:33378/141‑1 (MQ=255)                                              ccACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCGGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAAt                                                                                                             >  8:40087/1‑141 (MQ=255)                                              ccACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCGGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAAt                                                                                                             <  7:40087/141‑1 (MQ=255)                                              ccACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCaa                                                                                                                                         <  6:92575/113‑2 (MQ=255)                                              ccACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCaa                                                                                                                                         >  5:92575/1‑112 (MQ=255)                                                aCTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCtt                                                                                                     >  4:364927/1‑147 (MQ=255)                                                 cTGATGCATACGAATAATGTGTTGCGGCTGCGCGCTGATGCCGGGTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCACATGCCCGGCAGGTTGGTTCCAAAACGCTTTACCATTAATCGTCTCTTCa                                                                                                   <  6:430715/148‑1 (MQ=255)                                                 cTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTc                                                                                                    >  8:200486/1‑147 (MQ=255)                                                          aCGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTaa                                                                                           >  8:234771/1‑147 (MQ=255)                                                             aaTAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAAt                                                                                       >  5:303426/1‑148 (MQ=255)                                                              aTAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAAtt                                                                                      >  6:142643/1‑148 (MQ=255)                                                              aTAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAAtt                                                                                      <  3:44832/148‑1 (MQ=255)                                                                   gtgtTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTctct                                                                                                      >  8:420240/1‑127 (MQ=255)                                                                   gtgtTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTctct                                                                                                      <  7:420240/127‑1 (MQ=255)                                                                     gtTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTc                                                                               >  4:314894/1‑148 (MQ=255)                                                                      ttGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCt                                                                              >  7:163627/1‑148 (MQ=255)                                                                             tGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGca                                                                       >  7:76977/1‑148 (MQ=255)                                                                             tGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACAATTAATCGTCTCTTCAACGACTAAAAATTAGCCTACTGCGTGca                                                                       >  8:208820/1‑148 (MQ=255)                                                                              gcgcgcTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGcac                                                                      >  8:360979/1‑148 (MQ=255)                                                                                      aTGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACgca                                                                <  5:355565/146‑1 (MQ=255)                                                                                        gCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAg                                                                                                                                                    >  5:218457/1‑60 (MQ=255)                                                                                        gCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAg                                                                                                                                                    <  6:218457/60‑1 (MQ=255)                                                                                            gTTTCTTCCCACAGTTCACGGGCGGCGGCTTACACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCAGCGTGCACCACGCAAGAAACgg                                                        >  4:229709/1‑148 (MQ=255)                                                                                                  tcCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGt                                                  >  5:228162/1‑148 (MQ=255)                                                                                                    ccACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTg                                                 <  4:306239/147‑1 (MQ=255)                                                                                                     cacaGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGtc                                                                                                         <  7:316425/90‑1 (MQ=255)                                                                                                     cacaGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCg                                               >  1:236309/1‑148 (MQ=255)                                                                                                     cacaGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCg                                               >  5:13726/1‑148 (MQ=255)                                                                                                     cacaGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGGTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCg                                               >  2:211088/1‑148 (MQ=255)                                                                                                      acaGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGATCGTTCCATAACGCATTACCATTAATCGtct                                                                                                        >  8:316425/1‑90 (MQ=255)                                                                                                                 gcggcggcTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACAttt                                   <  3:364927/148‑1 (MQ=255)                                                                                                                     cggcTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCa                               <  7:360979/148‑1 (MQ=255)                                                                                                                      ggcTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCAt                              >  4:127446/1‑148 (MQ=255)                                                                                                                          tCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCt                           >  1:846/1‑147 (MQ=255)                                                                                                                               tAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCt                     >  6:198723/1‑148 (MQ=255)                                                                                                                               tAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCt                     >  4:102312/1‑148 (MQ=255)                                                                                                                                aaGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTg                    >  4:174806/1‑148 (MQ=255)                                                                                                                                  ggTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGt                   >  5:378023/1‑147 (MQ=255)                                                                                                                                     ttCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTg                                 >  1:220721/1‑130 (MQ=255)                                                                                                                                     ttCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTg                                 <  2:220721/130‑1 (MQ=255)                                                                                                                                       cATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTaa                        <  7:196195/137‑1 (MQ=255)                                                                                                                                       cATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTaa                        >  8:196195/1‑137 (MQ=255)                                                                                                                                       cATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTtc             >  4:354405/1‑148 (MQ=255)                                                                                                                                          cGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCg          >  4:88607/1‑148 (MQ=255)                                                                                                                                           ggCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAAt                       >  6:274308/1‑134 (MQ=255)                                                                                                                                           ggCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAAt                       >  6:268099/1‑134 (MQ=255)                                                                                                                                           ggCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAAt                       <  5:274308/134‑1 (MQ=255)                                                                                                                                           ggCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAAt                       <  5:268099/134‑1 (MQ=255)                                                                                                                                           ggCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGc         >  4:22417/1‑148 (MQ=255)                                                                                                                                            gCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGcc        <  3:174806/148‑1 (MQ=255)                                                                                                                                            gCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGcc        <  8:76977/148‑1 (MQ=255)                                                                                                                                                 cAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGCCATTCa   <  8:163627/148‑1 (MQ=255)                                                                                                                                                  aGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGCCATTCAc  <  7:208820/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCCGGTTTCTTCCCACAGTTCACGGGCGGCGGCTTCCACTAAGGTTTCATCGGCTTCCAGATGCCCGGCAGGTTGGTTCCATAACGCTTTACCATTAATCGTCTCTTCAACGACTAAAAATTTGCCTTCTGCGTGCACCACGCAAGCAACGGTAACGTGCGGTTTAAACATTTTGCATCCTTAATCTGTCACTTCTCGCCATTCAC  >  NC_000913/1194031‑1194323

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