Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,192,132	T→A	6.0%	intergenic (‑193/‑76)	glgS ← / → yqiJ	motility and biofilm regulator/DUF1449 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,192,132	0	T	A	6.0%	239.2 / 5.0	67	intergenic (‑193/‑76)	glgS/yqiJ	motility and biofilm regulator/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ACCTTCTGCGTGCATTCTGGCAAAACTACGCGCCAGGAAATCGAAATTATTTAAAGAATTAAGACTATGATCCATAAGCACATCCTCCGGTTATCGCATACCTTTTAAGGCCATATGCTTTTTCTGATTTCAGTATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATCGGCCTGTTAGAAATTTTTGCGCTTATCTGCGGTCACATGCTCTCCGGTGCTCTTGATGCACATCTTGACCATTATGATTCTATAACCACAGGTCATATTAGC  >  NC_000913/3191865‑3192372                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 aCCTTCTGCGTGCATTCTGGCAAAACTACGCGCCAGGAAATCGAAATTATTTAAAGAATTAAGACTATGATCCATAAGCACATCCTCCGGTTATCGCATACCTTTTAAGGCCATATGCTTTTTCTGATTTCAGTATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAac                                                                                                                                                                                                                                              <  2:402111/272‑1 (MQ=255)                                                                                                                              gATTTCAGTATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAacac                                                                                                                                                                                                                                            <  8:456517/149‑1 (MQ=255)                                                                                                                                tttCAGTATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTAATAATAACACAg                                                                                                                                                                                                                                          <  7:246042/149‑1 (MQ=255)                                                                                                                                      tATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAAt                                                                                                                                                                                                                                    <  5:28438/149‑1 (MQ=255)                                                                                                                                      tATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTAATAATAAcaca                                                                                                                                                                                                                                           >  6:297183/1‑142 (MQ=255)                                                                                                                                      tATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTAATAATAAcaca                                                                                                                                                                                                                                           <  5:297183/142‑1 (MQ=255)                                                                                                                                           cccAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTaa                                                                                                                                                                                                                                               <  8:227568/133‑1 (MQ=255)                                                                                                                                           cccAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTaa                                                                                                                                                                                                                                               >  7:227568/1‑133 (MQ=255)                                                                                                                                           cccAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAg                                                                                                                                                                                                                               >  5:50751/1‑149 (MQ=255)                                                                                                                                           cccAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAg                                                                                                                                                                                                                               >  8:126942/1‑149 (MQ=255)                                                                                                                                                         tGCCTTTAAACATAACGTTCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTa                                                                                                                                                                                                                 <  7:205789/149‑1 (MQ=255)                                                                                                                                                             tttAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACaaa                                                                                                                                                                                                             >  7:325209/1‑149 (MQ=255)                                                                                                                                                             tttAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACaaa                                                                                                                                                                                                             <  5:38005/149‑1 (MQ=255)                                                                                                                                                             tttAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACaaa                                                                                                                                                                                                             <  5:293707/149‑1 (MQ=255)                                                                                                                                                             tttAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACaaa                                                                                                                                                                                                             >  5:276934/1‑149 (MQ=255)                                                                                                                                                                aaaCATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAg                                                                                                                                                                                                                               >  7:31570/1‑128 (MQ=255)                                                                                                                                                                aaaCATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAg                                                                                                                                                                                                                               <  8:31570/128‑1 (MQ=255)                                                                                                                                                                aaaCATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTa                                                                                                                                                                                                          >  8:79970/1‑149 (MQ=255)                                                                                                                                                                    aTAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAata                                                                                                                                                                                                      >  7:45875/1‑149 (MQ=255)                                                                                                                                                                          tGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGt                                                                                                                                                                                                <  8:357565/149‑1 (MQ=255)                                                                                                                                                                                    tatTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATgg                                                                                                                                                                                      <  8:51064/149‑1 (MQ=255)                                                                                                                                                                                           cccGATCACACTTTGTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTgg                                                                                                                                                                               >  7:35643/1‑149 (MQ=255)                                                                                                                                                                                            ccGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAAc                                                                                                                                                                                                                      >  6:428872/1‑109 (MQ=255)                                                                                                                                                                                            ccGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAAc                                                                                                                                                                                                                      <  5:428872/109‑1 (MQ=255)                                                                                                                                                                                            ccGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGa                                                                                                                                                                              <  6:91767/149‑1 (MQ=255)                                                                                                                                                                                              gATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATaaa                                                                                                                                                                                                                                     >  8:436180/1‑92 (MQ=255)                                                                                                                                                                                              gATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATaaa                                                                                                                                                                                                                                     <  7:436180/92‑1 (MQ=255)                                                                                                                                                                                                 cacaCTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGat                                                                                                                                                                         <  8:45875/149‑1 (MQ=255)                                                                                                                                                                                                 cacaCTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGat                                                                                                                                                                         <  6:417665/149‑1 (MQ=255)                                                                                                                                                                                                   cacTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATgg                                                                                                                                                                                      >  4:515251/1‑134 (MQ=255)                                                                                                                                                                                                       ttttAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACAt                                                                                                                                                                                  <  2:515256/134‑1 (MQ=255)                                                                                                                                                                                                                 aaaTAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTTGc                                                                                                                                                         >  7:390889/1‑149 (MQ=255)                                                                                                                                                                                                                 aaaTAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGc                                                                                                                                                         <  8:325209/149‑1 (MQ=255)                                                                                                                                                                                                                  aaTAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGcc                                                                                                                                                        <  7:126942/149‑1 (MQ=255)                                                                                                                                                                                                                       aaaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGc                                                                                                                                                                                                  <  8:188073/102‑1 (MQ=255)                                                                                                                                                                                                                       aaaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGa                                                                                                                                                                   >  6:144088/1‑133 (MQ=255)                                                                                                                                                                                                                       aaaaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGa                                                                                                                                                                   <  5:144088/133‑1 (MQ=255)                                                                                                                                                                                                                       aaaaCCATATTAAAAGACTTTCTCTTACAGAAATTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAATGTAAATATTGc                                                                                                                                                                                                  >  7:188073/1‑102 (MQ=255)                                                                                                                                                                                                                         aaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAg                                                                                                                                                                             >  6:245302/1‑121 (MQ=255)                                                                                                                                                                                                                         aaCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAg                                                                                                                                                                             <  5:245302/121‑1 (MQ=255)                                                                                                                                                                                                                          aCCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATaa                                                                                                                                                <  6:50751/149‑1 (MQ=255)                                                                                                                                                                                                                             atatTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAacac                                                                                                                                             >  7:417355/1‑149 (MQ=255)                                                                                                                                                                                                                                   aaaGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTa                                                                                                                                       >  6:194671/1‑149 (MQ=255)                                                                                                                                                                                                                                   aaaGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTa                                                                                                                                       >  5:97519/1‑149 (MQ=255)                                                                                                                                                                                                                                         tttCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCtttt                                                                                                                                 >  5:24633/1‑149 (MQ=255)                                                                                                                                                                                                                                          ttCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATCGGCCTGTTAGAAAATTTTGCGCTTATCTGCGGTCACATGCTCTCCGGTGCTCTTGATGCACATCTTGACCATTATGCTTCTATCACCCAAGGTCATATTAGc  >  1:54131/1‑275 (MQ=255)                                                                                                                                                                                                                                           tctcTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAAt                                                                                                                                                                                        >  6:186911/1‑92 (MQ=255)                                                                                                                                                                                                                                           tctcTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAAt                                                                                                                                                                                        <  5:186911/92‑1 (MQ=255)                                                                                                                                                                                                                                             tctTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCt                                                                                                                             >  5:342511/1‑149 (MQ=255)                                                                                                                                                                                                                                             tctTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCt                                                                                                                             >  6:171436/1‑149 (MQ=255)                                                                                                                                                                                                                                               ttCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTAt                                                                                                                           >  7:55863/1‑149 (MQ=255)                                                                                                                                                                                                                                                   aGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTcc                                                                                                                       >  5:245089/1‑149 (MQ=255)                                                                                                                                                                                                                                                     aaCTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCtt                                                                                                                     <  6:276934/149‑1 (MQ=255)                                                                                                                                                                                                                                                      aCTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCttt                                                                                                                    <  6:24633/149‑1 (MQ=255)                                                                                                                                                                                                                                                      aCTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCttt                                                                                                                    <  8:35643/149‑1 (MQ=255)                                                                                                                                                                                                                                                        ttCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATg                                                                                                                                                                                       <  7:213484/80‑1 (MQ=255)                                                                                                                                                                                                                                                        ttCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATg                                                                                                                                                                                       >  8:213484/1‑80 (MQ=255)                                                                                                                                                                                                                                                         tctcTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTa                                                                                                                 <  6:342511/149‑1 (MQ=255)                                                                                                                                                                                                                                                           tctAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGAtt                                                                                                                                                                 >  7:469375/1‑99 (MQ=255)                                                                                                                                                                                                                                                           tctAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGAtt                                                                                                                                                                 <  8:469375/99‑1 (MQ=255)                                                                                                                                                                                                                                                                  tgtgTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATc                                                                                                        >  6:27554/1‑149 (MQ=255)                                                                                                                                                                                                                                                                  tgtgTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATc                                                                                                        <  8:417355/149‑1 (MQ=255)                                                                                                                                                                                                                                                                   gtgtTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATCg                                                                                                       <  7:79970/149‑1 (MQ=255)                                                                                                                                                                                                                                                                    tgtTATTATTAACACAGATCAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATCGGCCTGTTAGAAATTTTTGCGCTTATCTGCGGTCACATGCTCTCCGGTGCTCTTGATTCACAt                                         >  2:501562/1‑210 (MQ=255)                                                                                                                                                                                                                                                                    tgtTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTAtt                                                                                                                                                            >  3:20130/1‑95 (MQ=255)                                                                                                                                                                                                                                                                    tgtTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTa                                                                                                                                       >  2:268831/1‑116 (MQ=255)                                                                                                                                                                                                                                                                     gtTAATAATAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATNACACCCCTTATCTTTTCGCTATTNNCTTTGTACNTTTGATCGGCCTGTAAGAAATTTTTGCGCTTNTc                                                                             >  3:35436/1‑173 (MQ=255)                                                                                                                                                                                                                                                                       tattattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTGTTTCCTTTGTACTGTTGATCGGCCt                                                                                                   <  8:390889/149‑1 (MQ=255)                                                                                                                                                                                                                                                                       tattattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTCCTGTTGATCGGCCt                                                                                                   <  6:97519/149‑1 (MQ=255)                                                                                                                                                                                                                                                                        attattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGcc                                                                                                                                                        <  1:20130/95‑1 (MQ=255)                                                                                                                                                                                                                                                                        attattaACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCtt                                                                                                                                   <  4:268827/116‑1 (MQ=255)                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 ACCTTCTGCGTGCATTCTGGCAAAACTACGCGCCAGGAAATCGAAATTATTTAAAGAATTAAGACTATGATCCATAAGCACATCCTCCGGTTATCGCATACCTTTTAAGGCCATATGCTTTTTCTGATTTCAGTATTCCCCAATCAGTGTAGTGCCTTTAAACATAACGTGCGTAAATATATTGTCCCCGATCACACTTTTTAGTGAAAAATAAAAAACCATATTAAAAGACTTTCTCTTCCAGAACTTCTCTAACCTGTGTTATTATTAACACAGATAAATGTAAGCAAGGAAACTATTACAAAGTAAATATTGCGTCACTAAATGGACATTGGAGTGATATATGATTTTATTCGCCGACTATAACACCCCTTATCTTTTCGCTATTTCCTTTGTACTGTTGATCGGCCTGTTAGAAATTTTTGCGCTTATCTGCGGTCACATGCTCTCCGGTGCTCTTGATGCACATCTTGACCATTATGATTCTATAACCACAGGTCATATTAGC  >  NC_000913/3191865‑3192372

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