Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,770,911	T→A	S73C (AGT→TGT)	hemD ←	uroporphyrinogen‑III synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,770,911	0	T	A	100.0%	324.8 / NA	92	S73C (AGT→TGT)	hemD	uroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (45/47);  total (45/47)

CGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGCTTT  >  CP000730/1770777‑1771044                                                                                                                                       |                                                                                                                                      cGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGCTTTTTTCCATCATCCAAATGAAAAGTCGTTTGACATTAAATCATCTCGAATTTCAAGTGAGTCACAATATTGCGCTGTCTTACATCCTAt                                                                                                                                 <  1:435108/141‑1 (MQ=255)  gCCAATTCACTCGAAGGCAAAAGTATTTTTTTTTTAGTTTTATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATc                                                                                                                                >  2:364733/1‑141 (MQ=255)    cAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATcac                                                                                                                              >  2:634234/1‑141 (MQ=255)     aaTTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATcaca                                                                                                                             <  1:195000/141‑1 (MQ=255)         gacaattaaacAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCaa                                                                                                                         <  2:2693/131‑1 (MQ=255)         cACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCaa                                                                                                                         >  2:489373/1‑141 (MQ=255)         cACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCaa                                                                                                                         >  2:213873/1‑141 (MQ=255)           cTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATg                                                                                                                       <  1:296183/141‑1 (MQ=255)           cTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATg                                                                                                                       >  1:682274/1‑141 (MQ=255)            tCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGt                                                                                                                      >  2:491028/1‑141 (MQ=255)            tCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGt                                                                                                                      >  1:350604/1‑141 (MQ=255)             cGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGtt                                                                                                                     <  2:111578/141‑1 (MQ=255)             cGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATGTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAAGATTGCGCTGTCTTACAGCCTATCACAGCAATGtt                                                                                                                     <  2:256930/141‑1 (MQ=255)                aGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATc                                                                                                                  <  1:937444/141‑1 (MQ=255)                  gCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCaa                                                                                                                <  2:570535/141‑1 (MQ=255)                   cAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAAc                                                                                                               >  2:539223/1‑141 (MQ=255)                   cAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAAc                                                                                                               >  2:217083/1‑141 (MQ=255)                     aaaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACAt                                                                                                             >  2:928535/1‑141 (MQ=255)                      aaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACAtt                                                                                                            <  2:198690/141‑1 (MQ=255)                      aaGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACAtt                                                                                                            <  1:210927/141‑1 (MQ=255)                       aGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTa                                                                                                           <  2:30750/141‑1 (MQ=255)                        gTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTaa                                                                                                          <  2:900479/141‑1 (MQ=255)                        gTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTaa                                                                                                          >  2:669674/1‑141 (MQ=255)                           ttttttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTc                                                                                                       <  2:865327/141‑1 (MQ=255)                             ttttGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCt                                                                                                     >  1:716187/1‑141 (MQ=255)                                ttgTTATTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCtttt                                                                                                  >  2:575530/3‑141 (MQ=255)                                  gTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAg                                                                                                >  1:797297/1‑141 (MQ=255)                                   ttAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGa                                                                                               <  1:284933/141‑1 (MQ=255)                                           aTTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTtata                                                                                       <  1:547630/141‑1 (MQ=255)                                             tAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATaaa                                                                                     >  2:803927/1‑141 (MQ=255)                                             tAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCAATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATaaa                                                                                     <  2:923173/141‑1 (MQ=255)                                              aaaTGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATaaaa                                                                                    >  1:293822/1‑141 (MQ=255)                                              aaaTGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATaaaa                                                                                    >  1:105523/1‑141 (MQ=255)                                               aaTGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAg                                                                                   <  1:948339/141‑1 (MQ=255)                                                aTGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGa                                                                                  >  2:158000/1‑141 (MQ=255)                                                   aTTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAAtt                                                                               >  2:892630/1‑141 (MQ=255)                                                   aTTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAAtt                                                                               <  1:314590/141‑1 (MQ=255)                                                   aTTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAAtt                                                                               <  2:763638/141‑1 (MQ=255)                                                      ttAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTca                                                                            <  1:588071/141‑1 (MQ=255)                                                        aaaaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTcaca                                                                          <  2:556787/141‑1 (MQ=255)                                                        aaaaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTcaca                                                                          >  2:344362/1‑141 (MQ=255)                                                         aaaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAg                                                                         <  2:874712/141‑1 (MQ=255)                                                          aaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATCTTGCGCTGTCTTACATCCTATCACAGCAATGTCATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGc                                                                        <  1:725803/141‑1 (MQ=255)                                                          aaaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGc                                                                        <  1:531778/141‑1 (MQ=255)                                                           aaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACATCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCa                                                                       <  1:746190/141‑1 (MQ=255)                                                           aaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCa                                                                       >  1:782045/1‑141 (MQ=255)                                                           aaTCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCa                                                                       <  2:741553/141‑1 (MQ=255)                                                              ccTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAttt                                                                    >  1:250633/1‑141 (MQ=255)                                                               cttcttGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAtttt                                                                   >  2:302273/1‑141 (MQ=255)                                                                ttcttGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAttttt                                                                  >  1:747421/1‑141 (MQ=255)                                                                     gagaAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATg                                                             <  2:219583/141‑1 (MQ=255)                                                                     gagaAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATg                                                             <  1:344362/141‑1 (MQ=255)                                                                       gaAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGaa                                                           <  2:496731/141‑1 (MQ=255)                                                                            gTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAAt                                                      >  1:756333/1‑141 (MQ=255)                                                                             tCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATa                                                     >  2:570730/1‑141 (MQ=255)                                                                              cGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATaa                                                    <  1:599369/141‑1 (MQ=255)                                                                              cGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATaa                                                    >  2:418308/1‑141 (MQ=255)                                                                              cGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATACTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATaa                                                    >  2:838852/1‑141 (MQ=255)                                                                               gTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAg                                                   >  2:739184/1‑141 (MQ=255)                                                                               gTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAg                                                   <  1:397393/141‑1 (MQ=255)                                                                                tttGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACACAATTTTTAGATGAAAAAATAAGc                                                  >  2:44804/1‑141 (MQ=255)                                                                                 ttGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACCGCATTTTTAGATGAAAAAATAAGcc                                                 <  2:252981/141‑1 (MQ=255)                                                                                    gCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGt                                              >  1:581924/1‑141 (MQ=255)                                                                                    gCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAATTTTAGAAGAAAAAATAAGCCAGt                                              <  2:883336/141‑1 (MQ=255)                                                                                      aTAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCa                                            <  2:746213/141‑1 (MQ=255)                                                                                      aTAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCa                                            <  1:11366/141‑1 (MQ=255)                                                                                      aTAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCa                                            <  2:858006/141‑1 (MQ=255)                                                                                        aaaaTCAGCTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATa                                          <  1:622937/141‑1 (MQ=255)                                                                                            tCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGt                                      >  1:312983/1‑141 (MQ=255)                                                                                             cAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGtt                                     >  1:247164/1‑141 (MQ=255)                                                                                             aaaatCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATAACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCAAAAAAAAATAAAAAAATAAGCCAGTCATAGCGtt                                     <  2:198949/137‑1 (MQ=255)                                                                                               aCTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCATTGa                                   >  1:245637/1‑141 (MQ=255)                                                                                                cTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGAt                                  >  2:427417/1‑141 (MQ=255)                                                                                                cTCGAATGCCAAGTGATTCACAATATTGCGCTGTCGTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGAt                                  <  1:770058/141‑1 (MQ=255)                                                                                                  cGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGAttt                                >  1:493104/1‑141 (MQ=255)                                                                                                  cGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGAttt                                >  1:571327/1‑141 (MQ=255)                                                                                                     aTGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAAc                             >  2:963671/1‑141 (MQ=255)                                                                                                      tGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACa                            <  1:254691/141‑1 (MQ=255)                                                                                                      tGCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACa                            <  2:29833/141‑1 (MQ=255)                                                                                                       gCCAAGTGATTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGAAGAAAAAATAAGCCAGTCATAGCGTTGATTTAACaa                           <  2:170057/141‑1 (MQ=255)                                                                                                               aTTCACAATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAatat                   >  2:356524/1‑141 (MQ=255)                                                                                                                   acaATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCaaa               >  1:776879/1‑141 (MQ=255)                                                                                                                    caATATTGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCaaaa              <  1:623867/141‑1 (MQ=255)                                                                                                                         ttGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAg         <  1:372540/141‑1 (MQ=255)                                                                                                                          tGCGCTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGt        >  1:426616/1‑141 (MQ=255)                                                                                                                            cgcTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTgg      <  2:376174/141‑1 (MQ=255)                                                                                                                            cgcTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTgg      >  2:252355/1‑141 (MQ=255)                                                                                                                             gcTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGc     <  1:494179/141‑1 (MQ=255)                                                                                                                              cTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGCt    >  1:775918/1‑141 (MQ=255)                                                                                                                              cTGTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAAGTTAGTGGCt    <  1:853345/141‑1 (MQ=255)                                                                                                                                gTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGCttt  >  1:183673/1‑141 (MQ=255)                                                                                                                                gTCTTACATCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGCttt  >  1:183732/1‑141 (MQ=255)                                                                                                                                       |                                                                                                                                      CGCCAATTCACTCGAAGGCAAAAGTATTTTTTGGTTAGTTTGATTAAATGATTTTAAAAATCCTTCTTGAGAAAAGTCGTTTGGCATAAAATCAACTCGAATGCCAAGTGATTCACAATATTGCGCTGTCTTACTTCCTATCACAGCAATGTTATCAACATTAATTCCTTTTAGATATTTATAAAAGAATTTCACAGCATTTTTAGATGAAAAAATAAGCCAGTCATAGCGTTGATTTAACAAATGAATATCAAAATTTAGTGGCTTT  >  CP000730/1770777‑1771044

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