Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,385,206	T→A	5.7%	P456P (CCT→CCA)	yjeM →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,385,206	0	T	A	5.7%	286.5 / 5.9	87	P456P (CCT→CCA)	yjeM	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base A (3/2);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385058‑4385353                                                                                                                                                     |                                                                                                                                                    tCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCt                                                                                                                                                     >  1:27983/1‑149 (MQ=255)    tCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                  >  8:146968/1‑149 (MQ=255)      gcgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGt                                                                                                                                                <  8:138865/149‑1 (MQ=255)       cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  7:213152/149‑1 (MQ=255)          tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTg                                                                                                                                              >  4:73523/1‑147 (MQ=255)              tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  7:175813/1‑149 (MQ=255)              tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  7:84642/1‑149 (MQ=255)              tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGACGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaa                                                                                                                                            >  1:143086/1‑145 (MQ=255)              tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGACGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaa                                                                                                                                            <  2:143086/145‑1 (MQ=255)                  ccGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgag                                                                                                                                     <  1:64576/148‑1 (MQ=255)                    gTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACt                                                                                                                                  >  2:92454/1‑149 (MQ=255)                       cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                <  1:171118/148‑1 (MQ=255)                       cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                <  2:27983/148‑1 (MQ=255)                        tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  7:71814/1‑149 (MQ=255)                          tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  6:204674/1‑147 (MQ=255)                            aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGc                                                                                                                          <  8:48199/149‑1 (MQ=255)                                  gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 <  1:12042/136‑1 (MQ=255)                                  gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  2:12042/1‑136 (MQ=255)                                  gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAg                                                                                                                           >  3:130356/1‑142 (MQ=255)                                  gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAg                                                                                                                           <  4:130356/142‑1 (MQ=255)                                  gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtg                                                                                                                    <  5:55241/149‑1 (MQ=255)                                     aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCCACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  7:218804/1‑148 (MQ=255)                                     aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  8:30423/1‑148 (MQ=255)                                     aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  1:207937/1‑148 (MQ=255)                                      ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGGTGTCGTACTGGTGGTGACATTTGCGAACGTCTGCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGga                                                                                                                <  4:151397/149‑1 (MQ=255)                                       gATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgat                                                                                                               <  8:2018/149‑1 (MQ=255)                                         tctcGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATAATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatga                                                                                                             >  2:132498/1‑149 (MQ=255)                                          ctcGACAGACCGTTTGTGATTTTCAAAACGCATTAGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCAATCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatgat                                                                                                            >  7:81974/1‑149 (MQ=255)                                                  aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  2:86112/106‑1 (MQ=255)                                                  aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               >  1:86112/1‑106 (MQ=255)                                                  aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  4:134138/1‑149 (MQ=255)                                                  aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  7:4115/1‑149 (MQ=255)                                                   ccGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGccc                                                                                                   >  8:140551/1‑149 (MQ=255)                                                        tgtgATTTTAATAACGCATTTGTCGGCAATTATTGCGACAGTTGTTGTCGTACTTGTGGTGACATTTGGGAACGTCTTCACCATCATTCATCCTGTGGTTGAAGCCTGAGACTGGGACAGCACTTTGTGTATGATTGGCGTCCCTGtc                                                                                               <  8:81974/148‑1 (MQ=255)                                                          tgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGGTGTCGTTCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttc                                                                                            <  1:237748/149‑1 (MQ=255)                                                           gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGga                                                                                                                                      <  5:152583/106‑1 (MQ=255)                                                           gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGga                                                                                                                                      >  6:152583/1‑106 (MQ=255)                                                              tttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttcttct                                                                                        >  1:127320/1‑149 (MQ=255)                                                               ttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                              <  8:131538/142‑1 (MQ=255)                                                               ttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                              >  7:131538/1‑142 (MQ=255)                                                                tCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTtctc                                                                                       <  8:134450/148‑1 (MQ=255)                                                                            tGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGAt                                                                           <  3:73523/148‑1 (MQ=255)                                                                                ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCCCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGa                                                                      >  2:44434/1‑149 (MQ=255)                                                                                 gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                     >  4:101878/1‑149 (MQ=255)                                                                                  cAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAtt                                                                    <  8:4115/149‑1 (MQ=255)                                                                                    aTGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTa                                                                  >  4:166833/1‑149 (MQ=255)                                                                                     tGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAc                                                                 <  6:191429/149‑1 (MQ=255)                                                                                       aTTGCGACAGTGGTTGTCGTTCTGGTGGTGACATTTGCGAACGTCTTCTCCATGATTGAACCTGTGGTTGAAGCCGGAGTCTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTTGCGATGGCGATTTAcct                                                               <  5:197561/149‑2 (MQ=255)                                                                                       aTTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCa                                                               <  7:146968/149‑1 (MQ=255)                                                                                        ttGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAg                                                              <  4:113048/149‑1 (MQ=255)                                                                                           cGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGaa                                                            <  8:84642/148‑1 (MQ=255)                                                                                              cAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAt                                                        <  3:134138/149‑1 (MQ=255)                                                                                                gTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATAATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      <  1:132498/149‑1 (MQ=255)                                                                                                 tGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      <  2:207937/148‑1 (MQ=255)                                                                                                   gTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    <  1:98983/148‑1 (MQ=255)                                                                                                       tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCaga                                               <  1:92454/149‑1 (MQ=255)                                                                                                        cGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCagag                                              >  6:235189/1‑149 (MQ=255)                                                                                                         gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                     <  4:8408/141‑1 (MQ=255)                                                                                                         gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                     >  3:8408/1‑141 (MQ=255)                                                                                                         gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGt                                             >  6:115357/1‑149 (MQ=255)                                                                                                         gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGt                                             <  4:164605/149‑1 (MQ=255)                                                                                                           aCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcagca                                                 <  4:69183/143‑1 (MQ=255)                                                                                                           aCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcagca                                                 >  3:69183/1‑143 (MQ=255)                                                                                                             tggtggTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGc                                          >  3:56443/1‑148 (MQ=255)                                                                                                                    gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGTCAAAAAATc                                  >  4:34351/1‑149 (MQ=255)                                                                                                                    gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATc                                  >  6:159197/1‑149 (MQ=255)                                                                                                                        tttGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGGCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCa                              >  8:63196/1‑149 (MQ=255)                                                                                                                        attGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGCGACAGCACAGTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGCTGGCGATTTACCAGAACTATTGCAGCAGAGAGGCAAAAAATCCg                                >  3:50663/2‑147 (MQ=255)                                                                                                                             gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                         >  6:28483/1‑149 (MQ=255)                                                                                                                                       aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  8:67904/1‑144 (MQ=255)                                                                                                                                       aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  7:67904/144‑1 (MQ=255)                                                                                                                                       aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAAt                >  3:233746/1‑148 (MQ=255)                                                                                                                                        ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAtt                                                       >  2:49048/1‑108 (MQ=255)                                                                                                                                        ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAtt                                                       <  1:49048/108‑1 (MQ=255)                                                                                                                                        ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATaa              <  3:166833/149‑1 (MQ=255)                                                                                                                                         catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             >  8:128758/1‑149 (MQ=255)                                                                                                                                         catcatTCAACATGTGGTTGAAGCCGGAGACTGGGACAGCAAATTGTGGATGTTTGGCGGCCCTGTCTTGTTCTCTATGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             <  4:50663/149‑1 (MQ=255)                                                                                                                                            gcttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  2:147985/3‑139 (MQ=255)                                                                                                                                            gcttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  1:147985/137‑1 (MQ=255)                                                                                                                                             cttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa         >  3:215337/2‑149 (MQ=255)                                                                                                                                             atTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGACGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa         >  1:248510/1‑149 (MQ=255)                                                                                                                                              ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  1:122196/137‑1 (MQ=255)                                                                                                                                              ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  2:122196/1‑137 (MQ=255)                                                                                                                                              ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt        <  7:228657/149‑1 (MQ=255)                                                                                                                                                 aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc     >  1:71570/1‑149 (MQ=255)                                                                                                                                                    cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGCTGGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  <  8:218804/149‑1 (MQ=255)                                                                                                                                                    cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    >  6:56541/1‑99 (MQ=255)                                                                                                                                                    cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    <  5:56541/99‑1 (MQ=255)                                                                                                                                                    cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  <  7:30423/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385058‑4385353

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