Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,406,861	T→A	intergenic (‑220/‑700)	lrhA ← / → alaA	DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,406,861	0	T	A	100.0%	146.8 / NA	40	intergenic (‑220/‑700)	lrhA/alaA	DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (17/23);  total (17/23)

CGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAGGGTCATTGCGAATTATCTGATGCAAAGTTATGTTATGTTAGGCAAAGTAATCTTCTA  >  NC_000913/2406621‑2407105                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     cGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCAGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCActct                                                                                                                                                                                                                                            <  1:220376/251‑1 (MQ=255)       agttGCACTTATCAGATATTATCTCGTACTGGCGGCTCATACTGAGCTGGCTAAAGCTGTGCACACACAAACAAGCAATTCCAGTTGTAATGTGACTCCCTGGCAGTTGCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTGGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTACGATACTCACTCTTTTCTa                                                                                                                                                                                                                                      <  2:423755/249‑1 (MQ=255)           gCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTAttt                                                                                                                                                                                                                                   <  1:340093/250‑1 (MQ=255)            cACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTAt                                                                                                                                                                                                                                 <  1:237645/251‑1 (MQ=255)            cACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGAGAACAAACAAACAAGCAATTCCATATGTAATGTGCCTACCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTACGTTATAATGCCTGAAGTAGAACACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTAt                                                                                                                                                                                                                                 >  1:280605/1‑251 (MQ=255)                tATCATATATTATCCCTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACCCACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGAc                                                                                                                                                                                                                              <  2:194008/250‑1 (MQ=255)                tATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGAc                                                                                                                                                                                                                              >  2:378652/1‑250 (MQ=255)                 aTCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACAt                                                                                                                                                                                                                            <  1:292013/251‑1 (MQ=255)                   cATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCGTTTCTATTTATGACATgc                                                                                                                                                                                                                          <  2:422121/251‑1 (MQ=255)                        attatCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCgt                                                                                                                                                                                                                      <  1:378652/250‑1 (MQ=255)                           gccaaaacggactggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGttt                                                                                                                                                                                                                  <  1:416073/238‑1 (MQ=255)                            tCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTg                                                                                                                                                                                                                 >  2:357300/1‑251 (MQ=255)                                                    gCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTg                                                                                                                                                                                         >  1:318987/1‑251 (MQ=255)                                                                        cacaAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATc                                                                                                                                                                     <  2:332427/251‑1 (MQ=255)                                                                                   aaTTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATAc                                                                                                                                                          >  1:8039/1‑251 (MQ=255)                                                                                    aTTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTATATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATAcg                                                                                                                                                         <  2:403006/251‑1 (MQ=255)                                                                                      tCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGGTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATAcgcg                                                                                                                                                       >  1:224091/1‑251 (MQ=255)                                                                                                atgggCCTCCCTGGCAGTTTCATCCGGATACCGGACGAAAGTAAAAATGCATAGGAGTTGCACTAAAAAAGCGACTCAGATTTTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGTTCGCATATCTATTAAGATACTCACTCTTTTCTATTTAGGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCGCGTGCAGCATCGCTGGTCATACGCGAACACGTAcc                                                                                                                                             <  1:223535/247‑1 (MQ=255)                                                                                                  gtgCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCa                                                                                                                                            >  1:105457/1‑250 (MQ=255)                                                                                                        cccTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCa                                                                                                                                      >  2:81488/1‑250 (MQ=255)                                                                                                            ggCAGGTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGt                                                                                                                                 <  2:122949/251‑1 (MQ=255)                                                                                                                         ggAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTCAACAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGCTCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGGGAGGCGGCATCa                                                                                                                     >  2:59818/1‑250 (MQ=255)                                                                                                                            aCCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATAc                                                                                                                                                          <  1:210169/210‑1 (MQ=255)                                                                                                                            aCCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATAc                                                                                                                                                          >  2:210169/1‑210 (MQ=255)                                                                                                                                  aCGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAggg                                                                                                            >  1:104988/1‑250 (MQ=255)                                                                                                                                    gaaagaaaAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATc                                                                                                         >  1:102132/1‑251 (MQ=255)                                                                                                                                     aaaGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATcc                                                                                                        <  1:81488/251‑1 (MQ=255)                                                                                                                                                gCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATcc                                                                                             >  2:393847/1‑251 (MQ=255)                                                                                                                                                                aaaaaaGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCTCTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCa                                                                              <  1:31422/250‑1 (MQ=255)                                                                                                                                                                aaaaaaGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGNATCCGTCATTTATCCGAGCATTTTACCCCa                                                                              >  2:226461/1‑250 (MQ=255)                                                                                                                                                                  aaaaGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCaaaa                                                                           <  2:14751/251‑1 (MQ=255)                                                                                                                                                                      gagACTGACAATGTTCCGTTATAATGCCTGAAGTAGAGCACCGAATAAATATTCAGGGATCGCATATCTGTTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCAATATCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTCCGGCGGCATTAAGAGCAGGGATACGTCATTCAACCGAGCATGTTACCCCAAAAACCt                                                                       >  2:478600/3‑251 (MQ=255)                                                                                                                                                                       cGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCtt                                                                      >  1:450997/1‑251 (MQ=255)                                                                                                                                                                            cacaTTGTTCCGCTATAATGCCTGAGGTAGATCACAGAATATATCTTCAGGGATCGCAGATCTATTAAGATACTCACTCTGTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTTCTTGTTTCACTCGCGTGCAGCCTCGCTGTTCATACGCGAACACGTACCAACAGCAATGGTGAGTCTGCAACAAGAGCAGTGATACGTCATTTATCCGAGCATTTTACCCCAAAAACCTtttatc                                                                 <  2:48663/251‑2 (MQ=255)                                                                                                                                                                                    tCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAgg                                                          <  1:357300/250‑1 (MQ=255)                                                                                                                                                                                     ccGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAggg                                                         <  1:252166/250‑1 (MQ=255)                                                                                                                                                                                         tagaaTGCCTGAAGTAGATCACAGAATATATCTTCATGGAGCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATTTAAATTTGAGTCCTTTTTCCACTCGCGTGCAGCATCGCTGGTCATACGCGAACACGTTCCAACAGCAATGGGGAGGCGGCATCAAGAGCATGTATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATCAGGATCAt                                                     <  2:174220/247‑1 (MQ=255)                                                                                                                                                                                                        aGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAGGGTCATTGCGAATTATCTGAt                                      <  1:345920/250‑1 (MQ=255)                                                                                                                                                                                                            cacaGAATATATCTTCAGGGATCGCATATCTATTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAGGGTCATTGCAAATTATCTGATGCaa                                  <  2:11445/250‑1 (MQ=255)                                                                                                                                                                                                                                           aTTAAGATACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTCTAAGGGTCATTGCGAATTATCTGATGCAANGTTATGTTATGTTAGGCAAAGTAAtcttctt  >  2:226042/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAACCTTTTATTTATAAGGGTCATTGCGAATTATCTGATGCAAAGTTATGTTATGTTAGGCAAAGTAATCTTCTA  >  NC_000913/2406621‑2407105

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