Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,912	A→C	K140T (AAA→ACA)	patA →	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,219,912	0	A	C	100.0%	275.7 / NA	72	K140T (AAA→ACA)	patA	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/34);  total (38/34)

TTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCT  >  NC_000913/3219764‑3220056                                                                                                                                                     |                                                                                                                                                 tttGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTacac                                                                                                                                                >  2:1072026/1‑151 (MQ=255)      aaTTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGaaa                                                                                                                                           <  2:623173/151‑1 (MQ=255)       aTTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAAt                                                                                                                                          <  2:314256/151‑1 (MQ=255)       aTTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAAt                                                                                                                                          >  2:282646/1‑151 (MQ=255)         tttCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATAc                                                                                                                                        >  1:339139/1‑151 (MQ=255)           tCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATAc                                                                                                                                        >  1:439019/1‑149 (MQ=255)           tCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACa                                                                                                                                       <  2:618030/150‑1 (MQ=255)            cAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGc                                                                                                                                     <  1:188848/151‑1 (MQ=255)               cGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGcttc                                                                                                                                  <  2:677783/151‑1 (MQ=255)                  ggggCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGcttctt                                                                                                                                <  2:1185960/150‑1 (MQ=255)                  ggggCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCttcttc                                                                                                                               <  1:1169598/151‑1 (MQ=255)                  gcggCCCCGTAGTCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGTGCGATGTTGGCGAAACCCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCttcttc                                                                                                                               <  2:52169/149‑1 (MQ=255)                    ggCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTg                                                                                                                             <  1:1579/151‑1 (MQ=255)                    ggCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTg                                                                                                                             <  2:1176671/151‑1 (MQ=255)                           tAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGc                                                                                                                      >  2:723630/1‑151 (MQ=255)                               ccggTTGTGTTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGc                                                                                                                   <  2:953923/147‑1 (MQ=255)                                cAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCAc                                                                                                                 >  1:507165/1‑151 (MQ=255)                                  gTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCg                                                                                                               >  2:361848/1‑151 (MQ=255)                                     gtgGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAg                                                                                                             <  1:417504/150‑1 (MQ=255)                                       ggTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTcc                                                                                                          <  2:139981/151‑1 (MQ=255)                                           tCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTg                                                                                                      >  2:423216/1‑151 (MQ=255)                                            ccgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTg                                                                                                      >  1:367204/1‑150 (MQ=255)                                            ccgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTg                                                                                                      >  1:405316/1‑150 (MQ=255)                                            ccgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCATTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTg                                                                                                      >  2:48726/1‑150 (MQ=255)                                             cgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCTAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGCCCGTTGa                                                                                                     >  2:273032/1‑150 (MQ=255)                                             cgccgTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGa                                                                                                     >  2:565284/1‑150 (MQ=255)                                                 gTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCGCCGAGTCCGTTGAagcag                                                                                                <  2:713467/151‑1 (MQ=255)                                                   aCAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAagcagc                                                                                               <  2:165676/150‑1 (MQ=255)                                                     aGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCt                                                                                            >  1:1143222/1‑151 (MQ=255)                                                       aaTCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTg                                                                                           >  1:834494/1‑150 (MQ=255)                                                           aaCTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGa                                                                                          >  2:7505/1‑147 (MQ=255)                                                               tGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTgg                                                                                    <  1:53527/149‑1 (MQ=255)                                                                        acCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTAc                                                                         >  1:320246/1‑151 (MQ=255)                                                                         ccGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGATGGCGAAAGCTTAcc                                                                        >  2:49540/1‑151 (MQ=255)                                                                         ccGCAGCAAAGCAAGGAGGTGGTCGATCCGGTACGGGCGATGTTGGCGAAAACCCTCGCTCAGCTAACGCCCGGTACACTGAAATACAGCTTTTTCTGTAATAGTCGCACCGAGTCCGTTGAACCACCGACGCAGCAGGCGACAGCTCAcc                                                                        >  2:860767/1‑151 (MQ=255)                                                                          cGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTAcc                                                                        >  2:801569/1‑150 (MQ=255)                                                                          cGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCa                                                                       <  1:440764/151‑1 (MQ=255)                                                                              gCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGt                                                                     <  1:940741/149‑1 (MQ=255)                                                                                   gCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcg                                                              <  2:163262/151‑1 (MQ=255)                                                                                    ccAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcgc                                                             >  1:991866/1‑151 (MQ=255)                                                                                      aGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCgg                                                           >  1:251155/1‑151 (MQ=255)                                                                                      aGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCgg                                                           >  2:789207/1‑151 (MQ=255)                                                                                       ggAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAACGCTTACCAGTCACCGCGCgg                                                           >  2:312039/1‑150 (MQ=255)                                                                                        gAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCa                                                         <  2:185012/151‑1 (MQ=255)                                                                                         aGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCa                                                         >  2:371204/1‑150 (MQ=255)                                                                                             gctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGttt                                                    >  2:986416/1‑151 (MQ=255)                                                                                              ctCGATCCGTTACGGGCGATGTTGGCAAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGttt                                                    >  1:344183/1‑150 (MQ=255)                                                                                                cGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACt                                                 >  1:685440/1‑151 (MQ=255)                                                                                                 gATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACtt                                                >  1:644275/1‑151 (MQ=255)                                                                                                  aTCCGTTACGGGCGATGTTGGCGAAACCCCTTGCTGCGCTAAGGGCCGGGACACTGAAATGCAGCTTCTTCTGTAATAGCGGCGCCGAGTCCGTTGATGCGGCGCTGAAGCTGGCGACATCTTACCAGTCACCGCGCGGCATGTTTACttt                                               <  2:620572/151‑1 (MQ=255)                                                                                                   tCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACtttt                                              >  1:1135169/1‑151 (MQ=255)                                                                                                    ccGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTa                                             <  1:730323/151‑1 (MQ=255)                                                                                                      gTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTAtt                                           <  2:1051436/151‑1 (MQ=255)                                                                                                        tACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGc                                         >  2:721554/1‑151 (MQ=255)                                                                                                           gggCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGccac                                      >  1:928056/1‑151 (MQ=255)                                                                                                            ggCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGccacc                                     <  2:425585/151‑1 (MQ=255)                                                                                                              cGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAg                                   <  2:593724/151‑1 (MQ=255)                                                                                                               gATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAgc                                  <  2:596396/151‑1 (MQ=255)                                                                                                                         aaaaCCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTcc                        >  1:146201/1‑151 (MQ=255)                                                                                                                           aaCCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCAc                      >  2:1113677/1‑151 (MQ=255)                                                                                                                            aCCCTTGCTGCGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACg                     <  2:281456/151‑1 (MQ=255)                                                                                                                                 tgctgcGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTa                  <  2:150174/149‑1 (MQ=255)                                                                                                                                   ctgcGCTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAAt               <  2:607733/150‑1 (MQ=255)                                                                                                                                     gcgcTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCa             <  2:1180597/150‑1 (MQ=255)                                                                                                                                      cgcTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACt           >  1:300974/1‑151 (MQ=255)                                                                                                                                      cgcTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACt           >  1:237063/1‑151 (MQ=255)                                                                                                                                        cTAACGCCCGGTGCACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACtt          <  2:808727/150‑1 (MQ=255)                                                                                                                                        cTAACGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTg         >  1:46899/1‑151 (MQ=255)                                                                                                                                          aaCGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgc       >  1:1097255/1‑151 (MQ=255)                                                                                                                                           aCGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACt           <  1:84003/146‑1 (MQ=255)                                                                                                                                           aCGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgcg      <  1:850166/151‑1 (MQ=255)                                                                                                                                            cGCCCGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgcgc     <  2:1063371/151‑1 (MQ=255)                                                                                                                                               ccGGTACACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCt  >  2:1174662/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTTGGAATTTTCAACGTGGGGCACCGTAATCCAGTTGTGGTTTCCGCCGTACAGAATCAACTTGCGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCT  >  NC_000913/3219764‑3220056

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