Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,515,898	T→G	V107V (GTT→GTG)	patD →	gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,515,898	0	T	G	100.0%	241.5 / NA	63	V107V (GTT→GTG)	patD	gamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/32);  total (31/32)

CCACGCCGAAAGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCG  >  NC_000913/1515750‑1516043                                                                                                                                                     |                                                                                                                                                  ccACGCCGAAAGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtt                                                                                                                                                 >  2:454103/1‑151 (MQ=255)   aCGCCGAAAGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttc                                                                                                                                               >  1:321023/1‑151 (MQ=255)       cGAAAGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcgctt                                                                                                                                           >  1:425933/1‑151 (MQ=255)         aaaGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcgctttt                                                                                                                                         <  2:693641/151‑1 (MQ=255)         aaaGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcgctttt                                                                                                                                         <  1:851799/151‑1 (MQ=255)          aaGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcgcttttt                                                                                                                                        >  1:531240/1‑151 (MQ=255)              gcgtgcgGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcgctttttcgc                                                                                                                                     <  2:900674/150‑1 (MQ=255)               cgtgcgGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTcgcg                                                                                                                                    <  2:829463/150‑1 (MQ=255)                gtgcgGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCggg                                                                                                                                  <  2:617342/151‑1 (MQ=255)                   cgGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTgc                                                                                                                               <  2:814565/151‑1 (MQ=255)                   cgGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTgc                                                                                                                               >  2:810050/1‑151 (MQ=255)                      aaTGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTgcggc                                                                                                                            >  2:660216/1‑151 (MQ=255)                      aaTGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTgcggc                                                                                                                            >  1:31461/1‑151 (MQ=255)                            tgctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTg                                                                                                                      >  2:514491/1‑151 (MQ=255)                             gctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGt                                                                                                                     >  1:368922/1‑151 (MQ=255)                              ctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTc                                                                                                                    >  1:145594/1‑151 (MQ=255)                              ctgAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCAGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGt                                                                                                                     >  2:854979/1‑150 (MQ=255)                                    cTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcg                                                                                                                                              <  2:585133/119‑1 (MQ=255)                                    cTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGtgtttcg                                                                                                                                              >  1:585148/1‑119 (MQ=255)                                     tGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATgg                                                                                                             >  2:826917/1‑151 (MQ=255)                                       gCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTc                                                                                                           <  2:839265/151‑1 (MQ=255)                                       gCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTc                                                                                                           <  1:772923/151‑1 (MQ=255)                                         tGATGTTATCGCAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTg                                                                                                         >  2:219573/1‑151 (MQ=255)                                          gATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCt                                                                                                          <  2:509333/149‑1 (MQ=255)                                          gATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTgg                                                                                                        <  1:32314/151‑1 (MQ=255)                                            tGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggcg                                                                                                      >  1:60061/1‑151 (MQ=255)                                             gTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTggcgg                                                                                                     >  2:394116/1‑151 (MQ=255)                                                aTCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCa                                                                                                   <  1:213099/150‑1 (MQ=255)                                                      gaaAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAat                                                                                            >  1:381664/1‑151 (MQ=255)                                                       aaaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAata                                                                                           <  1:563691/151‑1 (MQ=255)                                                       aaaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAat                                                                                            >  1:404255/1‑150 (MQ=255)                                                        aaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAatat                                                                                          <  1:838434/151‑1 (MQ=255)                                                        aaaTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAatat                                                                                          <  1:330744/151‑1 (MQ=255)                                                          aTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCt                                                                                        >  2:24116/1‑151 (MQ=255)                                                          aTGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCt                                                                                        >  2:277064/1‑151 (MQ=255)                                                           tGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCtt                                                                                       >  1:829705/1‑151 (MQ=255)                                                            ggTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCtt                                                                                       <  1:350468/150‑1 (MQ=255)                                                               cAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGaa                                                                                    <  1:896034/150‑1 (MQ=255)                                                                    ttttGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCa                                                                               <  2:683616/150‑1 (MQ=255)                                                                    ttttGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCa                                                                               <  2:179407/150‑1 (MQ=255)                                                                      ttGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATAc                                                                            >  2:731155/1‑151 (MQ=255)                                                                         ccGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTc                                                                         >  2:548336/1‑151 (MQ=255)                                                                           gAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCCTACTTCg                                                                        >  1:583918/1‑150 (MQ=255)                                                                            aaCTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCgat                                                                      >  1:861845/1‑151 (MQ=255)                                                                                  aGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCg                                                                >  2:92259/1‑151 (MQ=255)                                                                                   gTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGt                                                               <  2:166585/151‑1 (MQ=255)                                                                                    tCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTc                                                              >  1:708499/1‑151 (MQ=255)                                                                                     cccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTc                                                              <  2:139942/150‑1 (MQ=255)                                                                                     cccGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTc                                                              <  2:633960/150‑1 (MQ=255)                                                                                         tAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGAt                                                         >  1:46283/1‑151 (MQ=255)                                                                                               tgGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTg                                                   >  1:131246/1‑151 (MQ=255)                                                                                                    aaCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGt                                              <  2:709465/151‑1 (MQ=255)                                                                                                     aCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTc                                             <  2:491415/151‑1 (MQ=255)                                                                                                     aCCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTc                                             <  2:78756/151‑1 (MQ=255)                                                                                                      ccGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCg                                            <  2:60061/151‑1 (MQ=255)                                                                                                      ccGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCg                                            <  2:388679/151‑1 (MQ=255)                                                                                                        gCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTg                                          <  1:912773/151‑1 (MQ=255)                                                                                                           gCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCt                                       <  1:756230/151‑1 (MQ=255)                                                                                                            cATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCtt                                      <  1:855582/151‑1 (MQ=255)                                                                                                               aGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCt                                    <  2:389796/150‑1 (MQ=255)                                                                                                                        aaTGATGAAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGt                          <  1:162932/151‑1 (MQ=255)                                                                                                                          tgatgaAATCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTgg                        >  1:462858/1‑151 (MQ=255)                                                                                                                                 aTCCCGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATc                 >  2:855061/1‑151 (MQ=255)                                                                                                                                   cccGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCg               <  1:882267/151‑1 (MQ=255)                                                                                                                                   cccGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCg               >  2:80390/1‑151 (MQ=255)                                                                                                                                     cGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCt             <  1:548351/151‑1 (MQ=255)                                                                                                                                     cGGCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCt             >  1:539008/1‑151 (MQ=255)                                                                                                                                      ggCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtg            <  2:505001/151‑1 (MQ=255)                                                                                                                                       gCGATTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtga           <  2:702077/151‑1 (MQ=255)                                                                                                                                          aTTGTCGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgatg         >  2:392040/1‑150 (MQ=255)                                                                                                                                               cGATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGcc   >  1:408733/1‑151 (MQ=255)                                                                                                                                                gATGTGTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGGGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCcg  >  1:34185/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCACGCCGAAAGTGCGTGCGGAATGTCTGCTGAAACTGGCTGATGTTATCGAAGAAAATGGTCAGGTTTTTGCCGAACTGGAGTCCCGTAATTGTGGCAAACCGCTGCATAGTGCGTTCAATGATGAAATCCCGGCGATTGTCGATGTTTTTCGCTTTTTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCG  >  NC_000913/1515750‑1516043

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