Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	44.0%	123.3 / 19.5	75	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (35/7);  new base C (0/33);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 9.30e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGG  >  NC_000913/907455‑907747                                                                                                                                                     |                                                                                                                                                   tAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCg                                                                                                                                                >  2:1034594/1‑151 (MQ=255)  aGCGCGGTATCATCCGCCAGCAGATCGTGTTGCAGCCCCTGAATCAACGCCCTGGCGGGGGGGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGc                                                                                                                                               >  2:231372/1‑151 (MQ=255)   gcgcGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCt                                                                                                                                              <  1:324505/151‑1 (MQ=255)    cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTc                                                                                                                                             >  1:718542/1‑151 (MQ=255)     gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCa                                                                                                                                            <  2:1102627/151‑1 (MQ=255)      cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCAc                                                                                                                                           >  2:1459798/1‑151 (MQ=255)          atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGcc                                                                                                                                       >  1:373654/1‑151 (MQ=255)           tcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGc                                                                                                                                        <  2:489472/149‑1 (MQ=255)                  ccAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATg                                                                                                                               >  1:447716/1‑151 (MQ=255)                    agcagATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGtt                                                                                                                             >  1:527824/1‑151 (MQ=255)                        gttcttgtTTCACCCCCTAAATCACCGCCCGGGGGGTGGGGGGGGGTAGGAAAGTAACCACATTAAAAACCCACCCCAAACCCCGGGGGGGGGGGGGGGGGAGGGGGAACAACCGGGGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCaaa                                                                                                                         <  2:627884/146‑1 (MQ=255)                        gATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCaa                                                                                                                          >  1:1350856/1‑150 (MQ=255)                            gtgtTTCAGCCCCTGAATCAACGCCCTGGCGGTGGCGGGCGGTACGGGAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTAACCGCCCTAAAATGATCAAACTGt                                                                                                                     >  2:37107/1‑151 (MQ=255)                             tttttCAGCCCCTAAATCAACGCCCTGGGGGGGGGGGGGGGTAGGAAAAAAATCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGAGGGGGAGGGGGATCACCCGGGGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGt                                                                                                                     <  1:773246/148‑1 (MQ=255)                               tttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                                  >  1:234582/1‑151 (MQ=255)                                ttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                                 >  1:643426/1‑151 (MQ=255)                                tccacccccTAAACCACCGCCGGGGGGGTGGGGGGGGGTAGGGAAGAAACCCAATTAAAAACCCCCACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAACCAACCGGCGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                                  <  1:105779/145‑1 (MQ=255)                                    gCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                                              >  1:1387675/1‑150 (MQ=255)                                     ccccTGAACCACCGCCCGGGCGGTGGTGGGGGGTAGGGAAATAATCACATTGAAAAACCACCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAACCGGGGGC‑CCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                                 <  2:569100/146‑1 (MQ=255)                                          gAATCAACGCCCGGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAACCCACCCCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GTTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                        <  1:1022705/150‑1 (MQ=255)                                              cAACGCCCGGGGGGTGGGGGGCGGTCCGGAATTAAGAAAATAGAGAACCCCCCCGAAAACCCGGGGGGGGGGGAGGGGGTGGGGAATCACCCGGGGGC‑CCTTC‑CCGCAATCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCa                                                                                                     <  1:1204858/149‑1 (MQ=255)                                                    ccTGGCGGGGGAGGGCGGCACGGAAGTAATCACATTGAGAAACCCCACCGAAATCCAGCGGGGGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCCCACCCCCATAAAATGTTCCAACTGTTGCTGATAACTGAACACctct                                                                                                >  2:990237/1‑148 (MQ=255)                                                       gggggtggtGGGCGGTAGGGAAGTAATCAAATTGAGAACCCCCCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCACCCGGCGGC‑GCTTC‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGt                                                                                             <  2:777403/145‑1 (MQ=255)                                                       gggggtgggGGGCGGTACGAAAGTAAACAATTTGAAAACCCACACCAAATCCCGGGGGGTGGGGGGGGGGAGGGGGATCAACCGGCGGC‑GCTTC‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                           <  1:899225/147‑1 (MQ=255)                                                        ggggggggggggggTAGGGAAGTAATCACATTGAGAAACCACCCCAAAACCCGGGGGGTGGGGAGGGGGAGGGGGATCAACCAGCGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcgg                                                                                         <  2:117115/143‑1 (MQ=255)                                                          ggggggggggggTAGGAAAGTAATCACATTGAAACCCCACACAGAAACCGGGGGGGGGGGGGGGGGGGGGGGGATCAACCAGCGGC‑GCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                       <  1:579218/145‑1 (MQ=255)                                                           gtggtgGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                      >  2:1169026/1‑151 (MQ=255)                                                           gggggggggggTACGAAAGTAATCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCACCCGGCGGC‑GCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                      <  1:661860/146‑1 (MQ=255)                                                             ggtggGCGGTAGGGAAAAAACCACATTGAGAACCCACCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCACCCAGCGGC‑GCTTC‑CCGCCAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                      <  2:133028/149‑1 (MQ=255)                                                                ggggggtggggAAGTAACCAAATGGAAAACCCCCACCGAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGGC‑GTTCC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                  <  1:521867/141‑1 (MQ=255)                                                                  ggggTACGGAAGGAACCACATGGAAAACCCACCCCGAACCCGAGGGGGGGGGGGGGGGGAGGGGGATCACCCGGCGGC‑TCTTC‑CCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCtt                                                                                   <  2:508303/145‑1 (MQ=255)                                                                   cGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGAGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                                                >  2:1026856/1‑149 (MQ=255)                                                                        cGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCTCCTCTGTTCCGGCGGCTTCGAAGAt                                                                            >  2:1335126/1‑148 (MQ=255)                                                                         ggAAGTAATCACATTGAGACACCACACGGAAACCCAGGGGGGGGGGAGGGGGTGGGGGATCAACCAGCGCC‑GCATA‑CCCCCCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                            <  2:1338001/147‑1 (MQ=255)                                                                         gaaaGTAACCACATGGAGAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCACCCAGCGGC‑CTTCC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                          <  1:526992/147‑1 (MQ=255)                                                                            aGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                     >  1:1056065/1‑151 (MQ=255)                                                                             gTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                    >  1:721875/1‑151 (MQ=255)                                                                              tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                   >  2:318712/1‑151 (MQ=255)                                                                              tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                   >  2:1052784/1‑151 (MQ=255)                                                                               aaccacaTGGAGAACCCACCCCAAACCCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCGGCGGC‑GCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                                  <  1:424586/148‑1 (MQ=255)                                                                               aacatcaTGGAGAAACCACCCAGAAATCCTCGGGGTGGGGAGGGGAATGGGTATCAACCAGCGGC‑‑CTTACCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                                  <  1:1250494/146‑1 (MQ=255)                                                                                   acaTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTgg                                                              >  1:1297495/1‑151 (MQ=255)                                                                                     atggaaaaccccccccGAAATCCAGGGGGGGGGGAGGGGGATGGGGACCAACCAGCGGC‑GCTTA‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                                <  1:1189903/133‑1 (MQ=255)                                                                                      ttGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                           >  1:54436/1‑151 (MQ=255)                                                                                      taaagaACCCCCACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCACCCAGCGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                            <  2:1190438/147‑1 (MQ=255)                                                                                          gaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                       >  1:1046821/1‑151 (MQ=255)                                                                                          aaacccccacCAAAACCCGGGGGGGGGGGGGGGGGATGGGGATCACCCAGCGGC‑GTTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                          <  1:1299074/141‑1 (MQ=255)                                                                                           aacccACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGACCCACCCGGCGGC‑GCTTC‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatga                                                      <  2:977752/148‑1 (MQ=255)                                                                                           aaaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                       >  2:1450758/1‑150 (MQ=255)                                                                                            aaCCACCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                        <  2:976235/148‑1 (MQ=255)                                                                                                 cccagAAACCCAGGGGGTGGGGGGGGGGATGGGGATCAACCGGGGGC‑CCTTC‑CCGCACACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                   <  1:849415/144‑1 (MQ=255)                                                                                                   ccGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                               >  1:480560/1‑150 (MQ=255)                                                                                                     gAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                            >  2:99135/1‑151 (MQ=255)                                                                                                      aaacccgggggggggggggggggAGGGGGATCAACCAGCGGC‑GCTTC‑CCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                             <  1:250101/140‑1 (MQ=255)                                                                                                      aaaCCCAGGGGGGGGGGAGGGGGAGGGGGACCAACCGGCGGC‑GCTTC‑CCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                           <  2:1362114/151‑1 (MQ=255)                                                                                                       aaccgaGCGGGGGGGGGGGGGGGGGGGGATCACCCGGGGGC‑GTTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                          <  2:323678/146‑1 (MQ=255)                                                                                                         tCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                        >  1:467803/1‑151 (MQ=255)                                                                                                          ccAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                       >  1:1103408/1‑151 (MQ=255)                                                                                                               ggggggggggggggATGGGAACCAACCAGCGGC‑GCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                  <  1:474313/151‑1 (MQ=255)                                                                                                                gggggggggggggAGGGGGACCACCCGGCGGC‑GCTTC‑CCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg                                 <  2:925812/148‑1 (MQ=255)                                                                                                                gggggggggggggAGGGGAATCACCCAGCGGC‑GCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg                                 <  1:29072/148‑1 (MQ=255)                                                                                                                ggTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                   >  2:110189/1‑149 (MQ=255)                                                                                                                 ggggggggggggggggGAATAACCCAGGGGC‑TCTTC‑CCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGt                                <  2:284507/149‑1 (MQ=255)                                                                                                                     ggggggggAGGGGAACCAACCGGCGGC‑TCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcag                            <  2:447663/148‑1 (MQ=255)                                                                                                                        gggggagggggaTCAACCAGCGGC‑GCTTC‑CCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagca                          <  1:934504/150‑1 (MQ=255)                                                                                                                           ggatggggatCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                      >  2:514439/1‑151 (MQ=255)                                                                                                                           ggagggggatCAACCGGCGGC‑TCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                      <  1:659132/151‑1 (MQ=255)                                                                                                                               ggggatCACCCAGCGGC‑GCTTC‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                  <  1:346357/151‑1 (MQ=255)                                                                                                                               ggggaaAAACCGGCCGCGGCTCA‑CCTCTAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                   <  1:1159023/151‑1 (MQ=255)                                                                                                                                    tCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc             >  1:472989/1‑151 (MQ=255)                                                                                                                                      aaCCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa            >  2:411505/1‑150 (MQ=255)                                                                                                                                      aaCCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa            >  1:963412/1‑150 (MQ=255)                                                                                                                                            cggc‑gCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCg     >  2:148306/1‑151 (MQ=255)                                                                                                                                               c‑gcTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg   <  1:1014942/150‑1 (MQ=255)                                                                                                                                               c‑gcTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtgg  >  2:1251724/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGC‑GCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGG  >  NC_000913/907455‑907747

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