Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,197,901	A→G	Q62Q (CAA→CAG)	hemE →	uroporphyrinogen decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,197,901	0	A	G	100.0%	302.0 / NA	84	Q62Q (CAA→CAG)	hemE	uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/41);  total (43/41)

CACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAACCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGTTTTACCTCGCCAGTCACCTGCAAAGCTG  >  NC_000913/4197775‑4198037                                                                                                                               |                                                                                                                                         cACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCt                                                                                                                             <  1:165044/140‑1 (MQ=255) cACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCt                                                                                                                             <  2:696069/140‑1 (MQ=255) cACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCt                                                                                                                             >  1:525671/1‑140 (MQ=255)   cTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTAc                                                                                                                           >  2:109940/1‑140 (MQ=255)       aGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGc                                                                                                                       >  2:877805/1‑140 (MQ=255)        gTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCACCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCt                                                                                                                      <  1:203898/140‑1 (MQ=255)         tATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTg                                                                                                                     >  2:874444/1‑140 (MQ=255)         tATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTg                                                                                                                     >  1:381993/1‑140 (MQ=255)           tGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGa                                                                                                                   >  2:96548/1‑140 (MQ=255)             gatgatGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATg                                                                                                                 <  1:96548/140‑1 (MQ=255)              atgatgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgc                                                                                                                >  2:465521/1‑140 (MQ=255)                gatgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcgg                                                                                                              >  2:935033/1‑140 (MQ=255)                gatgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcgg                                                                                                              >  2:801/1‑140 (MQ=255)                 atgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggc                                                                                                             >  2:1028448/1‑140 (MQ=255)                 atgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggc                                                                                                             >  2:586154/1‑140 (MQ=255)                  tgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggcg                                                                                                            >  1:297828/1‑140 (MQ=255)                  tgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggcg                                                                                                            >  2:679554/1‑140 (MQ=255)                  tgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggcg                                                                                                            >  2:803883/1‑140 (MQ=255)                  tgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggcg                                                                                                            >  2:216546/1‑140 (MQ=255)                  tgCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATgcggcg                                                                                                            >  2:941134/1‑140 (MQ=255)                   gcgcCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGa                                                                                                           >  1:63392/1‑140 (MQ=255)                   gcgcCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGa                                                                                                           >  1:896033/1‑140 (MQ=255)                   gcgcCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGa                                                                                                           >  2:620378/1‑140 (MQ=255)                   gcgcCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGa                                                                                                           >  1:354962/1‑140 (MQ=255)                            gggTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCtttt                                                                                                  <  2:51185/140‑1 (MQ=255)                            gggTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCtttt                                                                                                  <  1:1028448/140‑1 (MQ=255)                              gTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTcc                                                                                                >  2:1022487/1‑140 (MQ=255)                              gTCGCTACCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTcc                                                                                                <  1:187416/140‑1 (MQ=255)                               tCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCg                                                                                               <  2:257958/140‑1 (MQ=255)                               tCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCg                                                                                               <  1:801/140‑1 (MQ=255)                                        aCCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTTACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCa                                                                                      >  1:428594/1‑140 (MQ=255)                                        aCCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCa                                                                                      <  1:935033/140‑1 (MQ=255)                                            gAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGt                                                                                  >  2:1037278/1‑140 (MQ=255)                                            gAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGt                                                                                  >  1:958544/1‑140 (MQ=255)                                             aaTATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTg                                                                                 >  1:144720/1‑140 (MQ=255)                                              atatAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGc                                                                                <  1:1022487/140‑1 (MQ=255)                                              atatAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGc                                                                                <  1:216546/140‑1 (MQ=255)                                               tataAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGcc                                                                               >  1:724370/1‑140 (MQ=255)                                                 taAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCgg                                                                             >  1:43618/1‑140 (MQ=255)                                                 taAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCgg                                                                             >  1:442420/1‑140 (MQ=255)                                                 taAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCgg                                                                             >  1:747946/1‑140 (MQ=255)                                                  aaaGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGa                                                                            <  1:941134/140‑1 (MQ=255)                                                  aaaGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGa                                                                            <  1:465521/140‑1 (MQ=255)                                                    aGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGAcg                                                                          <  2:63392/140‑1 (MQ=255)                                                    aGCCACGCGCGCCCAGGCAGGCGATTTTATGTAGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGAcg                                                                          <  2:958544/140‑1 (MQ=255)                                                     gCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGAcgc                                                                         <  1:874444/140‑1 (MQ=255)                                                     gCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGAcgc                                                                         <  2:354962/140‑1 (MQ=255)                                                      ccACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGa                                                                            >  2:910433/1‑136 (MQ=255)                                                      ccACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGa                                                                            <  1:910433/136‑1 (MQ=255)                                                       cACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGa                                                                       >  2:118815/1‑140 (MQ=255)                                                       cACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGa                                                                       <  2:137244/140‑1 (MQ=255)                                                        aCGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGAt                                                                      <  2:896033/140‑1 (MQ=255)                                                        aCGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGAt                                                                      <  1:679554/140‑1 (MQ=255)                                                        aCGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGAt                                                                      <  2:747946/140‑1 (MQ=255)                                                        aCGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGAt                                                                      <  1:1037278/140‑1 (MQ=255)                                                        aCGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGAt                                                                      <  1:109940/140‑1 (MQ=255)                                                         cgcgcgCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATg                                                                     >  1:1074398/1‑140 (MQ=255)                                                         cgcgcgCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATg                                                                     >  1:762263/1‑140 (MQ=255)                                                         cgcgcgCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATg                                                                     >  1:513028/1‑140 (MQ=255)                                                               cccAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTa                                                               >  2:100239/1‑140 (MQ=255)                                                                  aggcaggcGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAggg                                                            <  1:714036/140‑1 (MQ=255)                                                                    gcaggcGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGct                                                          <  2:236922/140‑1 (MQ=255)                                                                    gcaggcGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGct                                                          <  1:118815/140‑1 (MQ=255)                                                                           aTTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTAtttt                                                   <  2:1074398/140‑1 (MQ=255)                                                                            ttttATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTg                                                  <  2:542587/140‑1 (MQ=255)                                                                            ttttATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTg                                                  <  2:11902/140‑1 (MQ=255)                                                                               tATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAg                                               <  1:454767/140‑1 (MQ=255)                                                                                aTGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGTAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGc                                              <  1:456199/140‑1 (MQ=255)                                                                                aTGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGc                                              <  1:800969/140‑1 (MQ=255)                                                                                     gCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGgag                                         <  2:144720/140‑1 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGct                                                          >  2:748597/1‑122 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGct                                                          <  1:748597/122‑1 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGTCGgaga                                        >  2:684745/1‑140 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGgaga                                        >  2:720913/1‑140 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGgaga                                        >  1:160376/1‑140 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGgaga                                        >  2:1070078/1‑140 (MQ=255)                                                                                      cTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGgaga                                        >  1:752634/1‑140 (MQ=255)                                                                                                aCGCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGt                              <  1:803883/140‑1 (MQ=255)                                                                                                  gCCGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGttt                            <  1:720913/140‑1 (MQ=255)                                                                                                   ccGAGCTGGCGTGCGAAGTGAGTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGtttt                           >  2:769284/1‑140 (MQ=255)                                                                                                   ccGAGCTGGCGTGCGAAGTGACTTTGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGtttt                           >  2:625872/1‑140 (MQ=255)                                                                                                                           tGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGTTTTACCTCGCCAGTCACCTGCAAAGCt   <  1:625872/140‑1 (MQ=255)                                                                                                                           tGCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGTTTTACCTCGCCAGTCACCTGCAAAGCt   <  1:877805/140‑1 (MQ=255)                                                                                                                            gCAGCCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGTTTTACCTCGCCAGTCACCTGCAAAGCTg  <  2:843847/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                         CACTCCAGTATGGATGATGCGCCAGGCGGGTCGCTATCTACCGGAATATAAAGCCACGCGCGCCCAGGCAGGCGATTTTATGTCGCTGTGCAAAAACGCCGAGCTGGCGTGCGAAGTGACTTTGCAACCGCTGCGTCGCTACCCGCTGGATGCGGCGATCCTCTTTTCCGATATCCTCACCGTGCCGGACGCGATGGGGTTAGGGCTCTATTTTGAAGCCGGAGAAGGTCCGCGTTTTACCTCGCCAGTCACCTGCAAAGCTG  >  NC_000913/4197775‑4198037

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