Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,566,493	G→C	8.1%	A131G (GCC→GGC)	eutA ←	reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,566,493	0	G	C	8.1%	254.6 / 11.8	74	A131G (GCC→GGC)	eutA	reactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base C (4/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTT  >  NC_000913/2566347‑2566633                                                                                                                                                   |                                                                                                                                             ggCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGc                                                                                                                                             >  3:151844/1‑148 (MQ=255)     accaccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATc                                                                                                                                         >  3:131220/1‑148 (MQ=255)         ccGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  4:151844/148‑1 (MQ=255)            aCGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACAGCCGATGTCGATATTCAGTACCAGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGAtt                                                                                                                                  >  4:581790/1‑148 (MQ=255)              gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                <  3:581790/148‑1 (MQ=255)                                gCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              <  4:131220/148‑1 (MQ=255)                                 cTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  3:507171/124‑1 (MQ=255)                                 cTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     >  4:507171/1‑124 (MQ=255)                                      ttttCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  4:361546/1‑148 (MQ=255)                                         tCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        <  2:367104/145‑1 (MQ=255)                                         tCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  1:367104/1‑145 (MQ=255)                                                tCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                                                  <  4:966028/144‑1 (MQ=255)                                                tCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                                                  >  3:966028/1‑144 (MQ=255)                                                   aaCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCg                                                                                                                    <  3:441624/123‑1 (MQ=255)                                                   aaCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCg                                                                                                                    >  4:441624/1‑123 (MQ=255)                                                      aGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACgg                                                                                                                                     <  2:824277/103‑1 (MQ=255)                                                      aGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGATCGCCGGTCACgg                                                                                                                                     >  1:824277/1‑103 (MQ=255)                                                         gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                          >  1:228520/1‑95 (MQ=255)                                                         gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGAt                                                                                                                                          <  2:228520/95‑1 (MQ=255)                                                         gCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACt                                                                                     >  3:1064957/1‑148 (MQ=255)                                                          cGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    <  4:1064957/148‑1 (MQ=255)                                                              gTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     >  1:176026/1‑95 (MQ=255)                                                              gTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACgg                                                                                                                                     <  2:176026/95‑1 (MQ=255)                                                                   cgGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgcgc                                                                           >  4:45131/1‑148 (MQ=255)                                                                          aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAAt                                                                                                 >  4:479205/1‑119 (MQ=255)                                                                          aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAAt                                                                                                 <  3:479205/119‑1 (MQ=255)                                                                                    tCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAg                                                          >  1:143683/1‑148 (MQ=255)                                                                                     cGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGc                                                         >  2:519674/1‑148 (MQ=255)                                                                                       atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgc                                                                             >  2:841069/1‑126 (MQ=255)                                                                                       atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgc                                                                             <  1:841069/126‑1 (MQ=255)                                                                                       atatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAt                                                       >  1:387403/1‑148 (MQ=255)                                                                                        tatTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAtt                                                      >  4:372405/1‑148 (MQ=255)                                                                                          ttCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGc                                                    <  1:519674/148‑1 (MQ=255)                                                                                           tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                <  4:232720/71‑1 (MQ=255)                                                                                           tCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCg                                                                                                                                >  3:232720/1‑71 (MQ=255)                                                                                                aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTgg                                                                                                             <  2:679243/85‑1 (MQ=255)                                                                                                aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTgg                                                                                                             >  1:679243/1‑85 (MQ=255)                                                                                                aCCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGttt                                              <  2:387403/148‑1 (MQ=255)                                                                                                 cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGtttt                                             <  3:361546/148‑1 (MQ=255)                                                                                                 cccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGtttt                                             <  2:143683/148‑1 (MQ=255)                                                                                                  ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCggg                                                              >  2:1054212/1‑130 (MQ=255)                                                                                                  ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCggg                                                              <  1:1054212/130‑1 (MQ=255)                                                                                                     acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    <  4:1029730/105‑1 (MQ=255)                                                                                                     acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTg                                                                                    >  3:1029730/1‑105 (MQ=255)                                                                                                           ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAgg                                                                                                                             >  1:1017598/1‑58 (MQ=255)                                                                                                           ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAgg                                                                                                                             <  2:1017598/58‑1 (MQ=255)                                                                                                           ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              <  2:821735/73‑1 (MQ=255)                                                                                                           ccGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTg                                                                                                              >  1:821735/1‑73 (MQ=255)                                                                                                             gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCTGGCg                                                            <  2:662367/121‑1 (MQ=255)                                                                                                             gttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCg                                                            >  1:662367/1‑121 (MQ=255)                                                                                                              ttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagag                                                                                <  1:197332/100‑1 (MQ=255)                                                                                                              ttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagag                                                                                >  2:197332/1‑100 (MQ=255)                                                                                                                   cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                                           >  4:592104/1‑132 (MQ=255)                                                                                                                   cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                                           <  3:592104/132‑1 (MQ=255)                                                                                                                   cAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAGCGGtgat                           >  2:77911/1‑148 (MQ=255)                                                                                                                     gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgccg                                                                >  3:151594/1‑109 (MQ=255)                                                                                                                     gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgccg                                                                <  4:151594/109‑1 (MQ=255)                                                                                                                     gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatga                         <  3:372405/148‑1 (MQ=255)                                                                                                                      agaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTg                                                     <  1:731652/117‑1 (MQ=255)                                                                                                                      aaaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTg                                                     >  2:731652/1‑119 (MQ=255)                                                                                                                       aaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        <  2:709790/67‑1 (MQ=255)                                                                                                                       aaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACg                                                                                                        >  1:709790/1‑67 (MQ=255)                                                                                                                            tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGc                  <  3:45131/148‑1 (MQ=255)                                                                                                                            tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGCGCGAGCATTGCGGGTTTTCGCGCTTTCa                                   >  2:294420/1‑131 (MQ=255)                                                                                                                            tttGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTAGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCa                                   <  1:294420/131‑1 (MQ=255)                                                                                                                                     cGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATc         >  1:49221/1‑148 (MQ=255)                                                                                                                                     cGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATc         >  4:993567/1‑148 (MQ=255)                                                                                                                                       gCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCaa       <  1:77911/148‑1 (MQ=255)                                                                                                                                         tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgatga                      <  3:522667/131‑1 (MQ=255)                                                                                                                                         tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgatga                      >  4:522667/1‑131 (MQ=255)                                                                                                                                         tCCGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACg     >  4:992314/1‑148 (MQ=255)                                                                                                                                          ccGTGATCGCCGGTCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGc    >  3:803444/1‑148 (MQ=255)                                                                                                                                          ccGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGc    <  3:993567/148‑1 (MQ=255)                                                                                                                                            gTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCtt  >  4:277872/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                             GGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGATGATGGCACCAGAATCAACGCTT  >  NC_000913/2566347‑2566633

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