Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,841	G→A	A450T (GCG→ACG)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,841	0	G	A	100.0%	126.1 / NA	43	A450T (GCG→ACG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/16);  total (27/16)

CGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTC  >  NC_000913/3220597‑3221082                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cGGTTTCTGTCAACTGGCGCGGGCATATCACGTTCTCGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGCGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGCGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGcc                                                                                                                                                                                                                                              >  2:34078/1‑250 (MQ=255)       ccGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCg                                                                                                                                                                                                                                         >  1:352167/1‑249 (MQ=255)                                     ggTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCg                                                                                                                                                                                                         >  1:186591/1‑251 (MQ=255)                                                     ggTAAAGGGATGTTTATGGCGAGTGAGTCTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGGTCCGCCAGCGCGTACTGGTGGCCGGATCGCTCAATAACGCCAAAACGATCCTCATTGAACCACCACTGACACTGCCCATTGAACAGTGTGAACTGGTGATCAAAGCGACGCTTAAGGCGCTGACGGCCATGCGAGTAAGTGGCGGAGAAGCGTAATATCCAGTCATATGGCGATacgacg                                                                                                                                                                                          <  2:253962/250‑1 (MQ=255)                                                            ggATGTTGTTGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCACATCGGATGGCGAACCGACGCCGCCATc                                                                                                                                                                                  >  2:203315/1‑251 (MQ=255)                                                              aTGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCg                                                                                                                                                                                >  1:37340/1‑251 (MQ=255)                                                                    aTGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCttt                                                                                                                                                                          >  1:353086/1‑251 (MQ=255)                                                                          aTTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATac                                                                                                                                                                    >  1:320642/1‑251 (MQ=255)                                                                               gTTTGTTGATAACGAACTCGGATATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAACCGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGAcgc                                                                                                                                                               >  2:262986/1‑251 (MQ=255)                                                                                     tgATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGc                                                                                                                                                         >  1:240919/1‑251 (MQ=255)                                                                                                 cGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGGGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGcatc                                                                                                                                              >  1:402100/1‑250 (MQ=255)                                                                                                  ggCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGcatcat                                                                                                                                            >  1:241295/1‑251 (MQ=255)                                                                                                          cTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCCTCATCCGCTCCg                                                                                                                                    >  2:350671/1‑251 (MQ=255)                                                                                                              gCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGt                                                                                                                                >  2:365660/1‑251 (MQ=255)                                                                                                                                cAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACt                                                                                                              >  1:181973/1‑251 (MQ=255)                                                                                                                                        aCTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACAATTGAACAGTGTGAACGGGGGAACAAAGCGGAGAGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTcct                                                                                                      >  2:66485/1‑249 (MQ=255)                                                                                                                                           ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCACCTGAAACTGACCATTGCACAGTGTGAACTGGTGATCAAACCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGCGTCGCAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCgg                                                                                                   >  1:31294/1‑251 (MQ=255)                                                                                                                                                       gCTCACTAACCCCAAAACGATCAGCATTGAACCGCCACTGACACTGACCATTGAACACTGTGAACTGGTGATCAAAGCGGCACGTAAGGCGCTAACGGCAATGCGAGTAAGCGTCGAAGAAGCGTAATATAAAATCGGATGGCGATACGACGTCGCCATCCGCTCTCTTTATACGACGCGTACGCCAGCAGGAATCATCCGCACCGGCGTTACCAGCCCGCTTTCCCTGCCGTAATCGGTTTCCGcaccc                                                                                        >  2:381535/1‑248 (MQ=255)                                                                                                                                                           aaacacgaaaaacggcaccgtatttaaCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACTCTTTCACTGCCGTCAGCGGTTTCCGCACACAACAg                                                                                   <  1:107407/226‑1 (MQ=255)                                                                                                                                                            aataagcaaaaccacactctattaaaCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAgc                                                                                  <  2:273804/227‑1 (MQ=255)                                                                                                                                                                 caaaatacaccccgtgttGACCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCCCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCTTTAACAGCACTCTTTCACTGCCGTCATCGTTTTCCGCCCACAACAgcatgca                                                                             <  2:338940/235‑1 (MQ=255)                                                                                                                                                                        cgatgcgactgcaacagCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATAGCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTAGACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGGCATCGGTTTCCGCACACAACAGCATGCATTCCGAc                                                                      <  1:367177/235‑1 (MQ=255)                                                                                                                                                                          acccataataacccGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGAc                                                                      <  1:155484/237‑1 (MQ=255)                                                                                                                                                                          aTCCGCATTGAACCGCCACTGACACTGACCATTGCACAGTGTGAACTGGTGATAAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTAGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGACTTCCGACGt                                                                    >  2:90647/1‑251 (MQ=255)                                                                                                                                                                           cccgaatgacaCCGCCACTGACCCTGACCATGTAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGtt                                                                   <  2:37340/241‑1 (MQ=255)                                                                                                                                                                            catctgccacccGCCACTGACACTGCCCATTGAACAGTGGGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAACTCGCATGGCGATACGACGTCGCCAGCCGATCTTTTTATACGACGCGTACGCCCTCAGGCATCCACCGCTCCGGCGTTAACCGCCCGCTTCCACTGCCGTCATCTGTTTCCGCAGACAGCTGCATGCTTTTCGACGtt                                                                   <  2:368144/240‑1 (MQ=255)                                                                                                                                                                                        gCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGTCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCAttg                                                      <  2:325061/251‑2 (MQ=255)                                                                                                                                                                                        gCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAAAAGCATGCATTCCGACGCTTCACCACGCAttt                                                      >  1:22469/1‑251 (MQ=255)                                                                                                                                                                                           actgacactgacCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCg                                                   >  1:188465/1‑251 (MQ=255)                                                                                                                                                                                             tgacactgacCATTTAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCACTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGAGATACGACGTCACCATCCGATAGTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCAGACGTTTCACCACGCATTTTCGCt                                                 >  2:109361/1‑251 (MQ=255)                                                                                                                                                                                              gacactgacCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCCCCACGCATTTTCCCtt                                                >  2:123559/1‑251 (MQ=255)                                                                                                                                                                                                    gacCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATGCCGACGTTTCACCACGCATTTTCGCTTTTTGCa                                          <  2:326198/251‑1 (MQ=255)                                                                                                                                                                                                     acCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAg                                         <  2:320642/251‑1 (MQ=255)                                                                                                                                                                                                      ccATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGa                                        <  1:157038/251‑1 (MQ=255)                                                                                                                                                                                                        aTTGAACAGTGTGAACTGGCGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGAtt                                      >  1:358165/1‑251 (MQ=255)                                                                                                                                                                                                           gAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGc                                    <  2:186591/250‑1 (MQ=255)                                                                                                                                                                                                             aCAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATa                                 <  1:408937/251‑1 (MQ=255)                                                                                                                                                                                                                   ggAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACATCAGCATGCATTCCGACGTTTCACCACGCATTGTCGCTTTTTGCAGATTGCATAAGacca                           <  2:180747/250‑1 (MQ=255)                                                                                                                                                                                                                                  aGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGAAGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGAGTTAAAAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTcccc             >  1:45520/1‑250 (MQ=255)                                                                                                                                                                                                                                  aGCGGCGCGTAAGGCGCTGACGGCCATGCGAGTAAGGGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTcccc             >  2:297499/1‑250 (MQ=255)                                                                                                                                                                                                                                     ggcgCGTAAGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCCAATCGGATGGCGATAAGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCCGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATc          >  2:380244/1‑250 (MQ=255)                                                                                                                                                                                                                                         cgTAAGGCGCTGACGGCCATGCGCGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCCTGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGctc     >  2:352277/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaGGCGCTGACGGCCATGCGAGTAAGTGTCGAAGAAGCACAATATCAAATCGGATGGCGATACGACGTCTCCATCCGATCTTTTTATACGACGCGTTCGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACGGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATGTTCGCGGTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCtcttc  <  2:13313/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTC  >  NC_000913/3220597‑3221082

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