Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,570,524	C→T	6.4%	A410T (GCT→ACT)	eutE ←	aldehyde oxidoreductase, ethanolamine utilization protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,570,524	0	C	T	6.4%	203.8 / 2.8	65	A410T (GCT→ACT)	eutE	aldehyde oxidoreductase, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CAGCGACGCAGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCATCAGTTCA  >  NC_000913/2570379‑2570669                                                                                                                                                  |                                                                                                                                                  caGCGACGCAGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCAt                                                                                                                                                 >  3:288533/1‑148 (MQ=255) caGCGACGCAGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCAt                                                                                                                                                 >  1:78194/1‑148 (MQ=255)       cGCAGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCa                                                                                                                                           <  3:79503/148‑1 (MQ=255)          aGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCt                                                                                                                                        >  1:57763/1‑148 (MQ=255)                   aaCGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGGAATGCGCGGTCCGGTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCAttt                                                                                                                               >  6:91756/1‑146 (MQ=255)                         cgcgcgCTGGTTACCCCTTCACCGGTTGGCGTGGGGATGGTCATGGTGGTCCAGCCTTCCCCGCCAAGACCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCg                                                                                                                          <  1:101030/147‑1 (MQ=255)                            gcgcTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGt                                                                                                                      <  3:280900/148‑1 (MQ=255)                                     aCCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAg                                                                                                              <  1:365570/147‑1 (MQ=255)                                          ttCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCACCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCAtt                                                                                                        >  3:18164/1‑148 (MQ=255)                                           tCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTg                                                                                                       >  5:177713/1‑148 (MQ=255)                                            cACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGc                                                                                                                                                   >  7:260845/1‑103 (MQ=255)                                            cACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGc                                                                                                                                                   <  8:260845/103‑1 (MQ=255)                                            cACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGTATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTgc                                                                                                      >  3:61852/1‑148 (MQ=255)                                            cACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGTATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTgc                                                                                                      >  3:66823/1‑148 (MQ=255)                                                     gCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGtgtg                                                                                             <  8:207620/148‑1 (MQ=255)                                                     gCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGAGCGAGTGCATTGCCGCCgtgt                                                                                              <  2:13115/147‑1 (MQ=255)                                                      cGTGGTGAAGCTCAAGGTGGTCCACCATTCCCATAAAAACCCCAGCCCGGCAATGCCCGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCTTGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTgtgg                                                                                            >  8:67425/1‑148 (MQ=255)                                                           tgATGGTCATGGTGGTCCAGCCTTCCCCGCCCCGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGc                                                                                        >  7:303128/1‑147 (MQ=255)                                                           tgATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCCCGGTCCGGTCTTTACGAAAATGCTGGTATCCAGAGCCTTCGGCATCGGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGGCGTGtggtggg                                                                                        >  2:134294/1‑146 (MQ=255)                                                           tgATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCa                                                                                       <  4:107185/148‑1 (MQ=255)                                                                    tggtggTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGcgcg                                                                                                                >  5:318421/1‑114 (MQ=255)                                                                    tggtggTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGcgcg                                                                                                                <  6:318421/114‑1 (MQ=255)                                                                     ggtggtCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGGCCTCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCCCCCGCCGTcc                                                                             >  7:354442/1‑148 (MQ=255)                                                                     ggtggtCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAAAGAAAATGCTGGTAACAATAGGATTCGCCATCTGGGTCATGTTTTAGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTc                                                                              >  7:313444/1‑147 (MQ=255)                                                                      gtggtCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACTGTCCTTTCTTATCGAAAAAGCTGGTTTCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGAAACCGCCTTcc                                                                             <  6:284353/147‑1 (MQ=255)                                                                       tggtCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAg                                                                           <  2:380342/148‑1 (MQ=255)                                                                             aGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCCCGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCAc                                                                     >  6:274604/1‑148 (MQ=255)                                                                                cTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCg                                                                   >  6:109414/1‑147 (MQ=255)                                                                                 ttCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCAGGTATCAATAGCATTCGCCATCTGGTACATGATTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTGTCACCGc                                                                  >  6:352035/1‑147 (MQ=255)                                                                                   ccccGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAg                                                               >  2:9966/1‑148 (MQ=255)                                                                                        cccagccccagcCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCaa                                                          <  4:192328/148‑1 (MQ=255)                                                                                         ccagccccagccCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAAt                                                         >  7:184029/1‑148 (MQ=255)                                                                                            gccccagcccGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCAGCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGc                                                      <  3:5463/148‑1 (MQ=255)                                                                                              cccagcccGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGtt                                                                         <  7:257910/127‑1 (MQ=255)                                                                                              cccagcccGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGtt                                                                         >  8:257910/1‑127 (MQ=255)                                                                                              cccagcccGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCAt                                                    >  7:98154/1‑148 (MQ=255)                                                                                                      ggCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACg                                            >  2:14692/1‑148 (MQ=255)                                                                                                        cAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGtt                                          <  4:18164/148‑1 (MQ=255)                                                                                                        cAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGtt                                          <  1:14692/148‑1 (MQ=255)                                                                                                            gCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGc                                       >  4:362953/1‑147 (MQ=255)                                                                                                            gCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGc                                       >  7:344840/1‑147 (MQ=255)                                                                                                             cACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGTATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCg                                      <  4:66823/147‑1 (MQ=255)                                                                                                             cACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGTATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCg                                      <  4:61852/147‑1 (MQ=255)                                                                                                             cACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCa                                                                                       >  5:52467/1‑98 (MQ=255)                                                                                                             cACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCa                                                                                       <  6:52467/98‑1 (MQ=255)                                                                                                                ggTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGAcgc                                  >  1:74701/1‑148 (MQ=255)                                                                                                                  tCCGTTCTTATCGTAACTGGTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGTTGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACgcgc                                <  8:344840/148‑1 (MQ=255)                                                                                                                    cGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCa                               <  8:184029/147‑1 (MQ=255)                                                                                                                       tCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCCCGTTGGCGACGCGCacgac                           >  1:34410/1‑148 (MQ=255)                                                                                                                            aCGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGc                       >  2:145578/1‑147 (MQ=255)                                                                                                                            aCGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCa                      <  2:74701/148‑1 (MQ=255)                                                                                                                                  aTGCTTGTATGAATAGCTTTCGCCATCTGGTTCATGTTTTCGATGATGCGCGAGTTCATTGCCGCCGTGTGGTGGCTACCGCGTTCCAGTTTCACAGCTAGCGCAATGGAATACTCCACGTTGGCGAGGATAACGACGGGGaaaaac                 <  8:313444/147‑5 (MQ=255)                                                                                                                                     cTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGca             <  5:109414/148‑1 (MQ=255)                                                                                                                                      tGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGcat            <  6:177713/148‑1 (MQ=255)                                                                                                                                         tATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCCCGACGGGCAACACCGGcatcat         >  4:212775/1‑148 (MQ=255)                                                                                                                                         tATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGcatcat         >  1:33253/1‑148 (MQ=255)                                                                                                                                          aTCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGcatcat         >  6:13366/1‑147 (MQ=255)                                                                                                                                           tCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGcatcat         >  1:277539/1‑146 (MQ=255)                                                                                                                                           tCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCatcatc        <  1:9966/147‑1 (MQ=255)                                                                                                                                            cAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcag      >  7:267444/1‑148 (MQ=255)                                                                                                                                              aTAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAAt                                                         <  3:172454/95‑1 (MQ=255)                                                                                                                                              aTAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAAt                                                         >  4:172454/1‑95 (MQ=255)                                                                                                                                              aTAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcagtt    <  3:285015/148‑1 (MQ=255)                                                                                                                                               ttgCATTCGCCATCTGTTTCATGTTTTCGATGTTGCTCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcagttc   <  6:140376/146‑1 (MQ=255)                                                                                                                                               tAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcagttc   <  7:213089/148‑1 (MQ=255)                                                                                                                                               tAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcagttc   <  8:98154/148‑1 (MQ=255)                                                                                                                                                 gCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCAtcagttca  <  5:341309/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CAGCGACGCAGACGGACAAACGTACGCGCGCTGGTTACCCCTTCACCGGTTGGCGTGGTGATGGTCATGGTGGTCCAGCCTTCCCCGCCCAGCCCCAGCCCGGCAATGCACGGTCCGTTCTTAACGAAAATGCTGGTATCAATAGCATTCGCCATCTGGTTCATGTTTTCGATGTTGCGCGAGTGCATTGCCGCCGTGTGGTGGCAACCGCCTTCCAGTTTCACCGCTAGCGCAATGGCATCCGCCACGTTGGCGACGCGCACGACGGGCAACACCGGCATCATCAGTTCA  >  NC_000913/2570379‑2570669

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