Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,523,466	T→A	6.9%	N210I (AAT→ATT)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,523,466	0	T	A	6.9%	201.1 / 4.9	58	N210I (AAT→ATT)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (34/20);  new base A (2/2);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

GCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAC  >  NC_000913/4523319‑4523613                                                                                                                                                    |                                                                                                                                                    gCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCAtt                                                                                                                                                    <  3:265083/149‑1 (MQ=255)  ccAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTg                                                                                                                                                   <  3:84810/149‑1 (MQ=255)       gCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGaa                                                                                                                                              <  2:564485/149‑1 (MQ=255)         tGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGc                                                                                                                                            >  4:468365/1‑149 (MQ=255)               tgGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTgcgc                                                                                                                                      >  1:4440/1‑149 (MQ=255)                 ggATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCa                                                                                                                                    >  3:612947/1‑149 (MQ=255)                        cAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCAATGGCGaattatc                                                                                                                                         <  1:420022/137‑6 (MQ=255)                        cAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCAATGGCGaattatc                                                                                                                                         >  2:420022/1‑132 (MQ=255)                          gTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGtt                                                                                                                              >  3:495853/1‑146 (MQ=255)                          gTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGtt                                                                                                                              <  4:495853/146‑1 (MQ=255)                                aTCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTg                                                                                                                     >  4:562006/1‑149 (MQ=255)                                  ccTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGcc                                                                                                                   <  3:468365/149‑1 (MQ=255)                                   cTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCg                                                                                                                  >  1:297664/1‑149 (MQ=255)                                     tCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGtt                                                                                                                <  4:67652/149‑1 (MQ=255)                                          caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACTGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtctt                                                                                                           >  4:615424/1‑149 (MQ=255)                                          caccacCTGAGATGTCCCGGCAGTGCCCAAAAATTGACTGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtctt                                                                                                           >  4:69922/1‑149 (MQ=255)                                          caccacCTGAGAGGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtctt                                                                                                           >  3:50773/1‑149 (MQ=255)                                                 tGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATc                                                                                                    <  3:562006/149‑1 (MQ=255)                                                                 tGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTc                                                                                    <  3:194777/149‑1 (MQ=255)                                                                          aTTGACAGCCGCCGCCGGATGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGccgccg                                                                           >  3:407139/1‑149 (MQ=255)                                                                             gACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCACGCCGCCgggg                                                                        <  3:393623/149‑1 (MQ=255)                                                                              aCAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGa                                                                       >  1:630933/1‑149 (MQ=255)                                                                                aGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGAcc                                                                     <  2:664962/149‑1 (MQ=255)                                                                                  ccgccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtcttct                                                                                                         <  4:172155/111‑1 (MQ=255)                                                                                  ccgccgccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtcttct                                                                                                         >  3:172155/1‑111 (MQ=255)                                                                                       gccgGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACgg                                                              <  1:352730/149‑1 (MQ=255)                                                                                           gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                                                          <  2:4440/149‑1 (MQ=255)                                                                                           gAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGc                                                          >  3:354995/1‑149 (MQ=255)                                                                                             agagGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGa                                                        <  2:270018/149‑1 (MQ=255)                                                                                                 gCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                                                    >  4:334100/1‑149 (MQ=255)                                                                                                 gCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgtt                                                    >  3:55042/1‑149 (MQ=255)                                                                                                  cacaGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgttg                                                   >  4:687163/1‑149 (MQ=255)                                                                                                   acaGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATTGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGAtgttgt                                                  >  2:330170/1‑149 (MQ=255)                                                                                                      ggCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGtgct                                               <  2:601054/149‑1 (MQ=255)                                                                                                       gCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtct                                                                                                            <  3:169392/87‑1 (MQ=255)                                                                                                       gCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGtct                                                                                                            >  4:169392/1‑87 (MQ=255)                                                                                                       gCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGtgctg                                              >  4:9641/1‑149 (MQ=255)                                                                                                        cTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGtgctgc                                             <  3:9641/149‑1 (MQ=255)                                                                                                         tttACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                            >  2:570374/1‑149 (MQ=255)                                                                                                         tttACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                            >  1:640492/1‑149 (MQ=255)                                                                                                         tctACAGGACGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                            >  2:244833/3‑149 (MQ=255)                                                                                                               ggacGCACGCCGTGCGGGCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGa                                      >  4:61731/4‑149 (MQ=255)                                                                                                                cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                            <  3:13118/142‑1 (MQ=255)                                                                                                                cccGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTgcg                                            >  4:13118/1‑142 (MQ=255)                                                                                                                  cGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTg                                   >  1:368276/1‑149 (MQ=255)                                                                                                                       ggCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCcatc                              <  2:368276/149‑1 (MQ=255)                                                                                                                        gcgtgcgtCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCcatca                             >  3:335476/1‑149 (MQ=255)                                                                                                                                 cTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCg                    >  3:538661/1‑149 (MQ=255)                                                                                                                                        gagaaTTGGCCGTTTGCGAAGCGTGTGCCAATGGTTGGCACCTTTCTGTTGTCTTCTCCATCGTTGGCGGGCAGCGTCGCGCCGCCGGGGACCAGCTAGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGccc             >  2:593601/4‑149 (MQ=255)                                                                                                                                        gaTAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGt                                    <  3:522681/126‑1 (MQ=255)                                                                                                                                        gaTAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGt                                    >  4:522681/1‑126 (MQ=255)                                                                                                                                         acgcttCGCCAATGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGcccc            >  1:252068/5‑149 (MQ=255)                                                                                                                                            aTTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAg                                     >  1:539542/1‑121 (MQ=255)                                                                                                                                            aTTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAg                                     <  2:539542/121‑1 (MQ=255)                                                                                                                                              tCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATc       >  2:391792/1‑149 (MQ=255)                                                                                                                                              tCGCCAATGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATc       <  2:252068/149‑1 (MQ=255)                                                                                                                                                 ccATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCgaa    >  4:480209/1‑147 (MQ=255)                                                                                                                                                   aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAc  >  2:240460/1‑149 (MQ=255)                                                                                                                                                   aTTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAc  >  4:395825/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCCAGGGCTGAATGTGGGATTGCCAGTCCCAATCCTTCGGCCACCACCTGAGATGTCCCGGCAGTGCCCAAAAATTGACAGCCGCCGCCGGAAGAGGCACAGGCTTTACAGCCCGCACGGCGTGCGTCCTGTAGAGATAATTCGCCATTGGCGAAGCGTGCGCCAATGGTTTGCACCTTGCCGTTGTCTTCTCCATCCTTTGCGGGCAGCGTCGCGCCGCCGGGGACCAGCACGGTTGCGATGTTGTGCTGCGCGGCGAGTGCCATCATGGTGGCCGGAAGCCCCTTATCGCAAC  >  NC_000913/4523319‑4523613

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