Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,220,673	G→C	5.4%	E394Q (GAG→CAG)	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,220,673	0	G	C	5.4%	247.1 / 2.0	70	E394Q (GAG→CAG)	patA	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base C (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCA  >  NC_000913/3220526‑3220820                                                                                                                                                    |                                                                                                                                                    gCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGa                                                                                                                                                    <  1:2161/149‑1 (MQ=255) gCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGa                                                                                                                                                    <  1:547617/149‑1 (MQ=255)   gACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGCTATGTTGCTGGCCGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGt                                                                                                                                                  <  2:193890/149‑1 (MQ=255)   gACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGt                                                                                                                                                  <  3:316777/149‑1 (MQ=255)   gACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGt                                                                                                                                                  >  3:461225/1‑149 (MQ=255)      cATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTttg                                                                                                                                               >  2:135221/1‑149 (MQ=255)        tCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTttgtt                                                                                                                                             <  1:135221/149‑1 (MQ=255)        tCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTttgtt                                                                                                                                             <  3:535257/149‑1 (MQ=255)          gaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTTCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGCAAAGGGATGTTGATGGCGATTCAGTTTGTTGa                                                                                                                                           >  3:479057/2‑149 (MQ=255)          aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGa                                                                                                                                           <  4:463328/149‑1 (MQ=255)          aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGa                                                                                                                                           >  2:165516/1‑149 (MQ=255)          aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGa                                                                                                                                           >  1:440226/1‑149 (MQ=255)          aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGa                                                                                                                                           >  1:80198/1‑149 (MQ=255)          aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTCAGTTTGTTGa                                                                                                                                           >  3:396077/1‑149 (MQ=255)                ttGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGa                                                                                                                                     >  1:107470/1‑149 (MQ=255)                   cTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGGCGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAAt                                                                                                                                  >  3:527245/1‑149 (MQ=255)                    tGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATc                                                                                                                                 <  4:381803/149‑1 (MQ=255)                      gAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCgg                                                                                                                               <  1:528696/149‑1 (MQ=255)                       aGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTCAGTTTGTTGATAACGAAATCGGc                                                                                                                              <  4:479057/149‑1 (MQ=255)                       aGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTCAGTTTGTTGATAACGAAATCGGc                                                                                                                              <  4:396077/149‑1 (MQ=255)                           gAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCtata                                                                                                                          >  1:599575/1‑149 (MQ=255)                           gAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCtata                                                                                                                          >  1:670470/1‑149 (MQ=255)                            aaCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATaa                                                                                                                         >  4:232299/1‑149 (MQ=255)                                 aCCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTg                                                                                                                    >  3:395262/1‑149 (MQ=255)                                  ccGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGc                                                                                                                   <  4:168258/149‑1 (MQ=255)                                    ggCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCa                                                                                                                 <  1:165516/149‑1 (MQ=255)                                        cAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCga                                                                                                             <  4:123114/149‑1 (MQ=255)                                        cAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCAATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCga                                                                                                             <  4:301434/149‑1 (MQ=255)                                                gCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTcc                                                                                                     <  3:548343/149‑1 (MQ=255)                                                    aaaGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCa                                                                                                 <  4:333839/149‑1 (MQ=255)                                                    aaaGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCa                                                                                                 >  1:383325/1‑149 (MQ=255)                                                          gATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGt                                                                                           >  3:1095/1‑149 (MQ=255)                                                           atatGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTa                                                                                          <  4:461225/149‑1 (MQ=255)                                                            tatGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCTTTCTGGTACAGGAAGCGCGTGGTAAAGAGATGTTGATGGAGATTGAGTGTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGCGATGTTCCGCCAGCGCGTAc                                                                                         <  2:534399/149‑1 (MQ=255)                                                               gTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACtgg                                                                                      >  1:326869/1‑149 (MQ=255)                                                                         gtttcccGCCAACTGGCGCGGGAATATCCCGATTAGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAAc                                                                            <  4:564718/144‑1 (MQ=255)                                                                         ggTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAAc                                                                            >  2:70209/1‑149 (MQ=255)                                                                         ggTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAAc                                                                            <  4:552018/149‑1 (MQ=255)                                                                           tttCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGc                                                                          <  2:107470/149‑1 (MQ=255)                                                                             tCCGTCAACTGGCGCGGGAATATCCCGATATGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGAGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTc                                                                        <  3:378879/149‑1 (MQ=255)                                                                                gTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                     >  1:587074/1‑149 (MQ=255)                                                                                gTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                     >  2:173530/1‑149 (MQ=255)                                                                                   aaCTGGTGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                  <  2:242961/149‑1 (MQ=255)                                                                                   aaCTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAgcgc                                                                                             <  3:243372/122‑1 (MQ=255)                                                                                   aaCTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAgcgc                                                                                             >  4:243372/1‑122 (MQ=255)                                                                                   aaCTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                  >  1:629617/1‑149 (MQ=255)                                                                                   aaCTGGCGCGGGAATAGCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                  <  3:506985/149‑1 (MQ=255)                                                                                              gAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGAt                                                       <  4:658769/149‑1 (MQ=255)                                                                                              gAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGAt                                                       <  3:116051/149‑1 (MQ=255)                                                                                                 tatCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCg                                                    <  2:143425/149‑1 (MQ=255)                                                                                                        aTCTGGTACAGGAAGCGCGTGGAAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGaa                                             >  4:550093/1‑149 (MQ=255)                                                                                                          cTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                           <  1:229840/149‑1 (MQ=255)                                                                                                              tACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCca                                       >  1:139604/1‑149 (MQ=255)                                                                                                               aCAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcac                                      >  1:370762/1‑149 (MQ=255)                                                                                                               aCAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcac                                      >  4:437114/1‑149 (MQ=255)                                                                                                                  ggAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactga                                   >  2:569440/1‑149 (MQ=255)                                                                                                                  ggAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactga                                   <  3:437114/149‑1 (MQ=255)                                                                                                                   gAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGTTGGCCGGAACGCTCAACAACGCCAAAACGATCCGCATTGAACCGCcactgac                                  >  1:432783/1‑149 (MQ=255)                                                                                                                      gcgcGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacact                               >  2:212306/1‑149 (MQ=255)                                                                                                                        gcgTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactga                             <  3:232299/149‑1 (MQ=255)                                                                                                                           tGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCa                          <  2:80198/149‑1 (MQ=255)                                                                                                                            ggTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCAt                         >  1:661976/1‑149 (MQ=255)                                                                                                                               aaaGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTg                       <  4:380494/148‑1 (MQ=255)                                                                                                                               aaaGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTg                       >  3:380494/1‑148 (MQ=255)                                                                                                                                 aGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                    <  1:569440/149‑1 (MQ=255)                                                                                                                                    gATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAgt                 <  2:164492/149‑1 (MQ=255)                                                                                                                                             ggCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTgg        <  2:349424/149‑1 (MQ=255)                                                                                                                                              gCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGt       <  4:388356/149‑1 (MQ=255)                                                                                                                                                 aTTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGAt    >  1:39064/1‑149 (MQ=255)                                                                                                                                                   tGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCa  >  4:547725/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCA  >  NC_000913/3220526‑3220820

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