Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	836,745	C→T	6.4%	A132V (GCA→GTA)	ybiC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	836,745	0	C	T	6.4%	244.6 / 4.6	78	A132V (GCA→GTA)	ybiC	putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/40);  new base T (3/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATT  >  NC_000913/836598‑836891                                                                                                                                                    |                                                                                                                                                   gcGATGGCGATCGCGCAGTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCa                                                                                                                                                   <  1:1132199/148‑1 (MQ=255)   gATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAgc                                                                                                                                                 >  2:1334100/1‑149 (MQ=255)      ggCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAgcggc                                                                                                                                              >  1:1304684/1‑149 (MQ=255)      ggCGATCGCGCATTTGGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAgcggc                                                                                                                                              <  2:1202945/149‑1 (MQ=255)         gATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCggg                                                                                                                                           >  1:1440577/1‑149 (MQ=255)             gcgcATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGttt                                                                                                                                       >  1:1005635/1‑149 (MQ=255)               gctgatggtctggtcGCGGCTCATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGt                                                                                                                                     <  1:24301/144‑1 (MQ=255)               gcATTTGGTCAGGGCGAGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGt                                                                                                                                     <  1:1076775/149‑1 (MQ=255)                 aTTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGtct                                                                                                                                   >  2:179263/1‑149 (MQ=255)                     ggTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTAt                                                                                                                               <  1:1213679/149‑1 (MQ=255)                     ggTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTAt                                                                                                                               <  1:1593322/149‑1 (MQ=255)                      gTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCCCGGTATTGCCGCCGTGGCGCTACATACCTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATc                                                                                                                              >  1:1494533/1‑149 (MQ=255)                         aGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCAc                                                                                                                           >  2:1107325/1‑149 (MQ=255)                            tCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACtt                                                                                                                         >  2:288963/1‑148 (MQ=255)                            tCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACtt                                                                                                                         <  1:288963/148‑1 (MQ=255)                            gcgcgGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACttt                                                                                                                        <  1:445088/148‑1 (MQ=255)                              gcgGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGt                                                                                                                      >  2:59544/1‑149 (MQ=255)                              gcgGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGt                                                                                                                      <  1:641523/149‑1 (MQ=255)                              gcgGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCCCCGTGGCGCTAAATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGt                                                                                                                      <  2:98950/149‑1 (MQ=255)                                ggCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTa                                                                                                                    >  1:980031/1‑149 (MQ=255)                                ggCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTa                                                                                                                    <  2:578453/149‑1 (MQ=255)                                 gCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCCTCAGCACGGTATTGCCGCCGTGGCGCGACATACCTCGCATAATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGATTGTCTCTATCCACGTTGTTAg                                                                                                                   >  1:263075/1‑149 (MQ=255)                                   acaTGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCg                                                                                                                 >  1:194216/1‑149 (MQ=255)                                       gAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCCCGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                                                                                             >  2:80722/1‑149 (MQ=255)                                        caGCGAGGGCGCTGGGGATTGAGAAAGCGCATCAGCAAGGTATAGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGGAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTc                                                                                                            <  2:605766/148‑1 (MQ=255)                                         aGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCg                                                                                                           >  1:382907/1‑149 (MQ=255)                                           cGATGGCGCTGGGGATTGAGAAAGCGCATGAGCACGGTATTGCCGCCGTGGCGGGACATAACTGGCATCATCTCGGCCGTATGGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGt                                                                                                         <  1:765347/149‑1 (MQ=255)                                           cGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTATTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGt                                                                                                         >  1:1217609/1‑149 (MQ=255)                                            taTGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTa                                                                                                        <  2:47947/148‑1 (MQ=255)                                              gggCGCTGGGGATTGAGACAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGAATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTAtt                                                                                                      <  2:115696/148‑1 (MQ=255)                                               ggCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGTGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTc                                                                                                     <  2:1592529/149‑1 (MQ=255)                                               ggCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTc                                                                                                     <  2:1212239/149‑1 (MQ=255)                                               ggCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTc                                                                                                     <  1:1334100/149‑1 (MQ=255)                                                   cTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAAt                                                                                                 >  2:958962/1‑149 (MQ=255)                                                    tGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATg                                                                                                <  2:896336/149‑1 (MQ=255)                                                      gggATTGAGAAAGTGCATCAGAACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGGATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTATGTTAGCGTGGTCGGTATTCCAATGGt                                                                                              >  2:264036/1‑149 (MQ=255)                                                             agaAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCCTATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGACg                                                                                       >  2:805849/1‑149 (MQ=255)                                                              gaAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTcgcgaccg                                                                                      >  1:894067/1‑145 (MQ=255)                                                              gaAAGCGCATCAGCAAGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGCCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACATTGTTAGCGTGGTCGGTATTCCAATGATCGCGCCGt                                                                                      >  1:716539/1‑149 (MQ=255)                                                                  gcgcATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCa                                                                                  >  1:1533982/1‑149 (MQ=255)                                                                         aGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTcgc                                                                           <  2:980031/149‑1 (MQ=255)                                                                              ggTATTGCCGCCGTGGAGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTTGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAg                                                                      <  2:174457/149‑1 (MQ=255)                                                                                    gccgccGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTGTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACCGCCGCtt                                                                >  2:781048/1‑149 (MQ=255)                                                                                     ccgccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCttt                                                               >  1:1538627/1‑149 (MQ=255)                                                                                       gccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCaa                                                                                                  <  2:642619/112‑1 (MQ=255)                                                                                       gccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCaa                                                                                                  >  1:642619/1‑112 (MQ=255)                                                                                       gccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTgg                                                             <  2:1434321/149‑1 (MQ=255)                                                                                        ccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGc                                                            >  1:340554/1‑149 (MQ=255)                                                                                          ggggCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCAc                                                          <  1:1295788/147‑1 (MQ=255)                                                                                          gTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCAc                                                          <  2:1416763/149‑1 (MQ=255)                                                                                              cgcTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAAt                                                      <  1:958962/149‑1 (MQ=255)                                                                                                 tACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCg                                                   >  2:1068662/1‑149 (MQ=255)                                                                                                  aCATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGt                                                  <  1:887784/149‑1 (MQ=255)                                                                                                   cATACCTCGCATCATATCGGCCGTATCGGTTACTGGGCGGACCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                                 <  1:89377/149‑1 (MQ=255)                                                                                                     tAACTCGCATCATATCGGCCGTATCGGTTCCTGGGCGGAGCAGTGGGCAGCGGCGGGGTGTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCt                                               <  2:804838/149‑1 (MQ=255)                                                                                                     tAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCt                                               <  2:893540/149‑1 (MQ=255)                                                                                                        cTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCtgtg                                            <  2:1585282/149‑1 (MQ=255)                                                                                                         tCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCtgtgt                                           >  1:1501108/1‑149 (MQ=255)                                                                                                           gCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTgg                                         >  2:579930/1‑149 (MQ=255)                                                                                                              tcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTTTTAGCGTGGTCGGTATTCAAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATACGTTCTGTGTGGttt                                      >  1:1017613/1‑149 (MQ=255)                                                                                                               catATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGtttt                                     <  2:40359/149‑1 (MQ=255)                                                                                                                  ctCGGCCGTATCGGTCACTGGGCGGAGCAGGGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGGCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCACTCCGTTCTGTGTGGTTTTccc                                  <  2:810421/148‑1 (MQ=255)                                                                                                                   tCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCt                                 <  1:179263/149‑1 (MQ=255)                                                                                                                    cGGCCGTCTCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTCTCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTc                                <  1:1569111/149‑1 (MQ=255)                                                                                                                    cGGCCGTATCGGTTCCTGGGCGGAGCAGTGTGCCGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTc                                <  1:318520/149‑1 (MQ=255)                                                                                                                          tATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAg                          <  1:910655/149‑1 (MQ=255)                                                                                                                                 tACTGGGCGGAGCAGTGTGTAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAAtct                   >  2:207828/1‑147 (MQ=255)                                                                                                                                  aCTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTc                  <  2:1533982/149‑1 (MQ=255)                                                                                                                                      ggCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGc              >  2:1408971/1‑149 (MQ=255)                                                                                                                                         ggAGCAGGGTGCAGCGGCGGGGTTTGGCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGt           <  2:1167787/149‑1 (MQ=255)                                                                                                                                          gAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGtt          >  1:1184417/1‑149 (MQ=255)                                                                                                                                             cAGTGTGTAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGa                         >  1:930552/1‑131 (MQ=255)                                                                                                                                             cAGTGTGTAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGa                         <  2:930552/131‑1 (MQ=255)                                                                                                                                             cAGTGTGTAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGa                         <  2:936515/131‑1 (MQ=255)                                                                                                                                             cAGTGTGTAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGa                         >  1:936515/1‑131 (MQ=255)                                                                                                                                               gtgtGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTg     >  2:1434035/1‑149 (MQ=255)                                                                                                                                                 gtgCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGAt   <  2:382907/149‑1 (MQ=255)                                                                                                                                                  tgCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGAtt  >  1:2948/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATT  >  NC_000913/836598‑836891

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: ‑