Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,635,457	T→C	5.1%	intergenic (‑147/‑64)	rob ← / → creA	right oriC‑binding transcriptional activator, AraC family/putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,457	0	T	C	5.1%	237.7 / 2.8	78	intergenic (‑147/‑64)	rob/creA	right oriC‑binding transcriptional activator, AraC family/putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base C (2/2);  total (40/38)
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.19e-01

CTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGA  >  NC_000913/4635319‑4635603                                                                                                                                           |                                                                                                                                                   cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTa                                                                                                                                                   <  6:182530/140‑1 (MQ=255) cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTa                                                                                                                                                   >  5:182530/1‑140 (MQ=255) cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           <  3:168916/148‑1 (MQ=255)       cTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTAtt                                                                                                                                     >  5:71514/1‑148 (MQ=255)         ttCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTAc                                                                                                                                  >  4:60088/1‑149 (MQ=255)              caacGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGa                                                                                                                             >  7:212452/1‑149 (MQ=255)                 cGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACg                                                                                                                           >  3:209824/1‑148 (MQ=255)                 cGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGa                                                                                                                          <  1:122736/149‑1 (MQ=255)                  gAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           >  3:202143/1‑131 (MQ=255)                  gAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           <  4:202143/131‑1 (MQ=255)                        cTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   >  4:110372/1‑149 (MQ=255)                            cATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  5:198771/149‑1 (MQ=255)                            cATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  2:103923/149‑1 (MQ=255)                             aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  6:46041/148‑1 (MQ=255)                             aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               >  5:46041/1‑148 (MQ=255)                             aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCg                                                                                                              >  5:201607/1‑149 (MQ=255)                              tCAGGTAATTGGATAATATCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCgt                                                                                                             <  5:102331/149‑1 (MQ=255)                              tCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAAt                                                                                                                                                 <  1:102498/113‑1 (MQ=255)                              tCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAAt                                                                                                                                                 >  2:102498/1‑113 (MQ=255)                                 ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          >  2:134090/1‑117 (MQ=255)                                 ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          >  2:69891/1‑117 (MQ=255)                                 ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          <  1:69891/117‑1 (MQ=255)                                 ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          <  1:134090/117‑1 (MQ=255)                                  gTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  4:8144/149‑1 (MQ=255)                                  gTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTATTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         >  5:210785/1‑149 (MQ=255)                                   tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  3:60088/148‑1 (MQ=255)                                   tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  5:190275/148‑1 (MQ=255)                                   tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAAt                                                                                                        <  2:30343/149‑1 (MQ=255)                                             ataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATaaaaa                                                                                               <  6:201607/148‑1 (MQ=255)                                               aGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            >  7:110940/1‑149 (MQ=255)                                                gCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            <  8:234872/148‑1 (MQ=255)                                                gCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACAGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            >  1:135837/1‑148 (MQ=255)                                                 cATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGt                                                                                           >  5:63364/1‑148 (MQ=255)                                                  aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         <  8:212452/149‑1 (MQ=255)                                                  aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         >  8:230506/1‑149 (MQ=255)                                                      tttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACaa                                                                                      <  3:110372/148‑1 (MQ=255)                                                      tttGCTGTTTTAGCATTCAGTGTTTTTATCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACaa                                                                                      <  2:206661/148‑1 (MQ=255)                                                       ttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAata                                                                                    >  1:182270/1‑149 (MQ=255)                                                          cTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGa                                                                                 <  7:121095/149‑1 (MQ=255)                                                                tAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATAc                                                                            <  6:71514/148‑1 (MQ=255)                                                                 aGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTGGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTaa                                                                                                                     <  2:98888/106‑1 (MQ=255)                                                                 aGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTGGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTaa                                                                                                                     >  1:98888/1‑106 (MQ=255)                                                                     ttCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAt                                                                      >  7:99008/1‑149 (MQ=255)                                                                      tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   >  5:65639/1‑103 (MQ=255)                                                                      tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   <  6:65639/103‑1 (MQ=255)                                                                        aGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAttt                                                                    >  5:170729/1‑148 (MQ=255)                                                                            tttttttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATc                                                                <  8:99008/148‑1 (MQ=255)                                                                                tttCTTAGTAGAGTATGTTTTAGGGCTGGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATGTGATCCTGTc                                                           >  4:55335/1‑149 (MQ=255)                                                                                 ttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCt                                                          >  3:161509/1‑149 (MQ=255)                                                                                    ttAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCttt                                                        >  2:24145/1‑148 (MQ=255)                                                                                       gtagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTATCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGc                                                    <  1:113037/149‑1 (MQ=255)                                                                                        tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                   >  2:170590/1‑149 (MQ=255)                                                                                        tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                   <  6:42816/149‑1 (MQ=255)                                                                                          gagTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCt                                                  >  4:17795/1‑148 (MQ=255)                                                                                             tATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGAt                                               >  3:63063/1‑148 (MQ=255)                                                                                                gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         >  8:216900/1‑103 (MQ=255)                                                                                                gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         <  7:216900/103‑1 (MQ=255)                                                                                                gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAAt                                            >  1:90171/1‑148 (MQ=255)                                                                                                 ttttAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGc                                          >  3:233217/1‑149 (MQ=255)                                                                                                  tttAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCt                                         >  6:78225/1‑149 (MQ=255)                                                                                                    tAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTg                                        <  1:24145/148‑1 (MQ=255)                                                                                                     aGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgg                                       <  4:63063/148‑1 (MQ=255)                                                                                                      gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgggg                                     <  3:17795/149‑1 (MQ=255)                                                                                                      gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      <  1:170590/148‑1 (MQ=255)                                                                                                      gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATGTGATCCTGTCTTTAAGCCTGATAATGCTgggg                                     >  4:48727/1‑149 (MQ=255)                                                                                                       ggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGc                                    <  4:161509/149‑1 (MQ=255)                                                                                                            gACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCAtt                                <  4:209824/148‑1 (MQ=255)                                                                                                             aCGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTgg                              >  5:100883/1‑149 (MQ=255)                                                                                                             aCGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGATAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTgg                              >  4:34604/1‑149 (MQ=255)                                                                                                               ggAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGATGGGGCCATTGGc                             >  2:80178/1‑148 (MQ=255)                                                                                                                gAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTc                           <  2:208215/149‑1 (MQ=255)                                                                                                                   aTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATg                        <  7:230506/149‑1 (MQ=255)                                                                                                                    tAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTACTTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGc                       <  5:78225/149‑1 (MQ=255)                                                                                                                     aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      <  7:154923/133‑1 (MQ=255)                                                                                                                     aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      >  8:154923/1‑133 (MQ=255)                                                                                                                      aaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCag                     >  2:45732/1‑149 (MQ=255)                                                                                                                                      gTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGt     >  3:230214/1‑149 (MQ=255)                                                                                                                                        tCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTc    >  5:183544/1‑148 (MQ=255)                                                                                                                                          ttAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGa  >  3:105278/1‑148 (MQ=255)                                                                                                                                           |                                                                                                                                                   CTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGA  >  NC_000913/4635319‑4635603

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