Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,093,463	A→G	12.0%	S270S (AGT→AGC)	rhaD ←	rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,463	0	A	G	12.0%	204.1 / 16.2	75	S270S (AGT→AGC)	rhaD	rhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (6/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTT  >  NC_000913/4093316‑4093605                                                                                                                                                    |                                                                                                                                               gTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCAc                                                                                                                                               >  2:323415/1‑149 (MQ=255)   aGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                             >  3:407925/1‑149 (MQ=255)    gTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTgg                                                                                                                                            <  1:358847/149‑1 (MQ=255)     ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           >  1:439416/1‑149 (MQ=255)     ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           >  2:493438/1‑149 (MQ=255)     ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           <  4:4381/149‑1 (MQ=255)     ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           >  1:376028/1‑149 (MQ=255)     ttATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           >  3:473671/1‑149 (MQ=255)       aTGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGa                                                                                                                                         >  2:260170/1‑149 (MQ=255)        tGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTAGGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGTGCGCTGGCgct                                                                                                                                        <  4:577524/149‑3 (MQ=255)          cgcACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCATCGCACTGGCGAGTg                                                                                                                                      <  1:448298/149‑1 (MQ=255)               cGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTa                                                                                                                                 >  2:84238/1‑149 (MQ=255)                     aTTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGCGTGGCGTAACGCCg                                                                                                                           >  3:68557/1‑149 (MQ=255)                        tCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGaaa                                                                                                                        <  2:119875/149‑1 (MQ=255)                            cACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCt                                                                                                                    <  2:148938/149‑1 (MQ=255)                             aCTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                   <  2:182264/149‑1 (MQ=255)                                  ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                              >  1:52474/1‑149 (MQ=255)                                  ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                              >  1:398021/1‑149 (MQ=255)                                  ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                              >  1:418006/1‑149 (MQ=255)                                   tCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAg                                                                                                             <  4:97893/149‑1 (MQ=255)                                          cTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATc                                                                                                      <  4:68557/149‑1 (MQ=255)                                            ggCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCaa                                                                                                    <  1:260170/149‑1 (MQ=255)                                                              cGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                  <  3:551203/149‑1 (MQ=255)                                                               gACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTc                                                                                 <  2:398021/149‑1 (MQ=255)                                                                  gCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                              >  1:335745/1‑149 (MQ=255)                                                                     gATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTc                                                                           <  3:331362/149‑1 (MQ=255)                                                                      aTACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCa                                                                          <  1:557855/149‑1 (MQ=255)                                                                                    cGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGaa                                                            <  1:323415/149‑1 (MQ=255)                                                                                     gTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             <  2:424455/147‑1 (MQ=255)                                                                                     gTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             >  1:424455/1‑147 (MQ=255)                                                                                         gggTTTACCGGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGACAGCGCACTGGCGAGTGGCTTAACGTCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                       <  2:209559/149‑1 (MQ=255)                                                                                            tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCtt                                                    <  2:376028/149‑1 (MQ=255)                                                                                              tACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCa                                                  <  2:52474/149‑1 (MQ=255)                                                                                                ccTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                             >  3:270110/1‑120 (MQ=255)                                                                                                ccTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                             <  4:270110/120‑1 (MQ=255)                                                                                                ccTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACt                                                >  3:328937/1‑149 (MQ=255)                                                                                                 cTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                               >  4:220026/1‑149 (MQ=255)                                                                                                 cTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                               >  1:476366/1‑149 (MQ=255)                                                                                                       gACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACtt                                         <  3:220026/149‑1 (MQ=255)                                                                                                         cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                       <  4:121280/149‑1 (MQ=255)                                                                                                         cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                       <  1:84238/149‑1 (MQ=255)                                                                                                           aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTGTCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                     >  1:320297/1‑149 (MQ=255)                                                                                                           aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                     >  2:281824/1‑149 (MQ=255)                                                                                                            aGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                         <  3:399118/112‑1 (MQ=255)                                                                                                            aGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                         >  4:399118/1‑112 (MQ=255)                                                                                                                tttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAtt                                >  4:108660/1‑149 (MQ=255)                                                                                                                 ttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAttt                               <  2:369614/149‑1 (MQ=255)                                                                                                                 ttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAttt                               <  1:132920/149‑1 (MQ=255)                                                                                                                   gCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAttttt                             <  2:418006/149‑1 (MQ=255)                                                                                                                    cGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTc                            <  3:289578/149‑1 (MQ=255)                                                                                                                     ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCt                           <  3:362631/149‑1 (MQ=255)                                                                                                                            gTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                       <  1:26390/130‑1 (MQ=255)                                                                                                                            gTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgt                                       >  2:26390/1‑130 (MQ=255)                                                                                                                             tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                   <  4:407925/149‑1 (MQ=255)                                                                                                                             tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                   <  4:240565/149‑1 (MQ=255)                                                                                                                                    cTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa            >  3:496720/1‑149 (MQ=255)                                                                                                                                     ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc           >  2:70211/1‑149 (MQ=255)                                                                                                                                     ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc           >  4:202557/1‑149 (MQ=255)                                                                                                                                       cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGGTGGTATGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCg         >  2:15440/5‑149 (MQ=255)                                                                                                                                       cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa            >  2:30107/5‑146 (MQ=255)                                                                                                                                       cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa            <  1:30107/142‑1 (MQ=255)                                                                                                                                       cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCg         >  3:292838/5‑149 (MQ=255)                                                                                                                                       cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCg         >  1:331696/5‑149 (MQ=255)                                                                                                                                        actcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGa        >  3:8591/4‑149 (MQ=255)                                                                                                                                         agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa       <  2:145958/149‑1 (MQ=255)                                                                                                                                          tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATCAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa            <  1:412159/142‑1 (MQ=255)                                                                                                                                          tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa            >  2:412159/2‑143 (MQ=255)                                                                                                                                          gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg      <  2:335745/149‑1 (MQ=255)                                                                                                                                          gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg      <  1:493438/149‑1 (MQ=255)                                                                                                                                          gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAg      >  3:372385/1‑149 (MQ=255)                                                                                                                                              cagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                >  3:452232/1‑87 (MQ=255)                                                                                                                                              cagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                <  4:452232/87‑1 (MQ=255)                                                                                                                                              cagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGttt  >  2:401185/1‑149 (MQ=255)                                                                                                                                                gcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa       <  2:298532/142‑1 (MQ=255)                                                                                                                                                gcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa       >  1:298532/1‑142 (MQ=255)                                                                                                                                                    |                                                                                                                                               GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTT  >  NC_000913/4093316‑4093605

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