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

GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907454‑907751                                                                                                                                                       |                                                                                                                                                      gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTAcc                                                                                                                                                     >  2:895918/1‑151 (MQ=255)     cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGc                                                                                                                                                   <  1:220287/149‑1 (MQ=255)      gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTc                                                                                                                                                 >  2:865915/1‑150 (MQ=255)      gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCa                                                                                                                                                <  2:1024885/151‑1 (MQ=255)        ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCAc                                                                                                                                               <  1:912090/150‑1 (MQ=255)          tATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGc                                                                                                                                            >  2:1006250/1‑151 (MQ=255)          tATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGGGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGc                                                                                                                                            >  1:471278/1‑151 (MQ=255)           atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGcc                                                                                                                                           >  2:661322/1‑151 (MQ=255)                   ccAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATg                                                                                                                                   >  2:989458/1‑151 (MQ=255)                        agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                               >  1:1092016/1‑150 (MQ=255)                          atggtgtTTCAGCCCCTGAATCAACGCCCGGGCGGTGGTGGGGGGTAC‑GGAAGAAATCACATGAAAAAACCACCCCGAACCCCGGGGGGGGGGGGGGGGGAGGGGAATCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCaaa                                                                                                                             <  2:72379/147‑1 (MQ=255)                           tCGTGTTTCAGCCCCTAAATCACCCCCCTGGGGGTGGTGGGCGGTAC‑GAAAGTAATCACATTGAGAACCCCCACCAAAATCCAGCGGGTGGGGAGGGGGAGGGGGATCAAC‑CAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCaa                                                                                                                              <  1:817619/148‑1 (MQ=255)                              tttttCAGCCCCTAAATCAACGCCCGGGGGGGGGGGGGGGGTAG‑GGAAGTAATCAAATTGAAAACCCCCACCAAAACCCGGGGGGGGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCaa                                                                                                                              <  2:445661/143‑1 (MQ=255)                              ttttccacccccTAAACCACCCCCCGGGGGGGGGGGGGGGGTAG‑GAAAGTAACCAAATGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGt                                                                                                                         <  2:127922/142‑1 (MQ=255)                              ttttccAGCCCCTAAATCAACGCCCTGGCGGGGGGGGGCGGTAC‑GGAAGTAACCACATTGAAAAACCCCACCAAACCCCAGGGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGt                                                                                                                         <  2:993068/145‑1 (MQ=255)                              tgtTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGtt                                                                                                                        >  2:402626/1‑151 (MQ=255)                                      ccccTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                                                 >  1:603252/1‑150 (MQ=255)                                           gAATCAACGCCCGGGCGGGGGGGGGGGGGCC‑GGAAGTAACCACATTGAAAACCCCCACCAAAACCCAGCGGGTGGGGGGGGGGATGGGGATCACC‑CGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                             <  2:893963/149‑1 (MQ=255)                                            aaTCAGCTCCTTGGGGGTGG‑GGGGGGTGCGGGAAGTAATATCAGTAAAAACCCCCCCCAAACCCCGGCGGGTGGGG‑GGGTGATGGGGGTCCACGCGGCGCCGTTAACCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAg                                                                                                           <  1:893385/150‑1 (MQ=255)                                                  cGCCCGGGGGGGGGTGGGCGGTAG‑GGAATTAATCACATAGAAAAACCACCCCAAAACCCGGGGGGTGGGGGGGGGGAGGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                      <  2:130511/149‑1 (MQ=255)                                                   gCCCTGGCGGTGGAGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                     >  1:812496/1‑149 (MQ=255)                                                      cTGGCGGTGGTGGGGGAGCC‑AGAATTAAACTCGGTGAAAACCCCCCCCGAAACCCGGCGGGGGGGGGGGGGGGTGGGGATAACC‑CAGCGCCGCTCGCCCCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGt                                                                                                 <  2:436275/150‑1 (MQ=255)                                                           ggtggtGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           >  2:746115/1‑151 (MQ=255)                                                           ggggggggggggTAG‑GGAAGAAACCACTTGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCACC‑GAGCGGCGTTTCCCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           <  1:416699/145‑1 (MQ=255)                                                              ggtggGCGGAGC‑GAAAGAAATCACATGAAGAACCCCCCCCAAAACCCGGGGGGTGGGGGGGGGGTTGGGTACCACC‑CAGCGGCGCTCCCCCCTCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCt                                                                                        <  1:576754/151‑1 (MQ=255)                                                                         c‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                             >  1:135997/1‑151 (MQ=255)                                                                           ggAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                            >  1:786914/1‑151 (MQ=255)                                                                           gaaattAATCACTTTAAGAACCCCCCCCGAATCCCGGGGGGTGGGGAGGGGGATGGGAATCACC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAg                                                                                  <  2:288096/140‑1 (MQ=255)                                                                            gAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                           >  2:483752/1‑151 (MQ=255)                                                                             aaGTAATCACATTGAAAACCCACACAAAAACCGGGGGGGGGGGGGGGGGGAGGGGAACCACC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                          <  1:1028103/151‑1 (MQ=255)                                                                                tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                       >  1:921687/1‑151 (MQ=255)                                                                                  aTCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGc                                                                     >  2:393194/1‑151 (MQ=255)                                                                                     acaTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTgg                                                                  >  1:886891/1‑151 (MQ=255)                                                                                     acaTTAAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGAATAAAC‑CGGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                                    <  1:1092164/149‑1 (MQ=255)                                                                                           agaAACCACCCAAAAACCCCGGGGGGGGGGGGGGGGGGGGGGTTCCCC‑CGGCGGCCCTCCCCCCTCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                            <  1:224199/151‑1 (MQ=255)                                                                                           agaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                            >  2:874324/1‑151 (MQ=255)                                                                                             aacccACCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGTATCACC‑CGGGGGCGCTCCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatga                                                          <  1:834934/148‑1 (MQ=255)                                                                                                ccACGCCGAAATCCGGGGGGGGGGGGGGGGGAGGGGGATCAGC‑CAGGGGCGCTTCCCTCTAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATc                                                        <  2:958292/150‑1 (MQ=255)                                                                                                 cacaCCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                       >  1:129133/1‑150 (MQ=255)                                                                                                    acCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaac                                                    >  1:254716/1‑150 (MQ=255)                                                                                                     ccGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                                   >  1:1068468/1‑150 (MQ=255)                                                                                                     ccGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                                   >  1:552535/1‑150 (MQ=255)                                                                                                        aaccccgggggggggggggggggggggggACCACC‑CGGGGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                                   <  1:87780/138‑1 (MQ=255)                                                                                                        aaacccgggggggggggggggggAGGGGGATCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                                 <  2:576561/140‑1 (MQ=255)                                                                                                        aaaCCCAGGGGGGGGGGGGGGGGAGGGGGACCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                                 <  2:223663/149‑1 (MQ=255)                                                                                                         actCCAGGGGGGGGGGGGGGGGAGGGGAATCAAC‑CGGCGGCTCTTCCCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                              <  2:218426/149‑1 (MQ=255)                                                                                                              gggggggggggggggggAGGGGAAACAAC‑CAGCGGCCCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                          <  1:499216/147‑1 (MQ=255)                                                                                                               ggggggggggggggggAGGGGGATCACC‑CGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAg                                         <  2:302933/148‑1 (MQ=255)                                                                                                                gggggggggggggggATGGGGATCAAC‑CGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                       <  1:549194/150‑1 (MQ=255)                                                                                                                   ggggggggggggATGGGAATAAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGt                                    <  2:322665/149‑1 (MQ=255)                                                                                                                     ggggAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGAACCACTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCGACCAGATAGTg                                   >  1:1058964/1‑150 (MQ=255)                                                                                                                       ggAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgc                                  >  2:380122/1‑149 (MQ=255)                                                                                                                          gggggagggggaTCACC‑CGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCaa                             <  1:310754/151‑1 (MQ=255)                                                                                                                           gggggggggaaccAAC‑CAGCGGCCCGTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAg                            <  1:708395/139‑1 (MQ=255)                                                                                                                            gggagggggatCAAC‑CGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGt                           <  2:993884/151‑1 (MQ=255)                                                                                                                               atggggatCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                        >  2:840099/1‑151 (MQ=255)                                                                                                                                 ggggatCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                      >  2:519347/1‑151 (MQ=255)                                                                                                                                 gggaatCAAC‑CAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                        <  2:558788/149‑1 (MQ=255)                                                                                                                                  gggaaccacccggcggcgCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                         <  1:648723/136‑1 (MQ=255)                                                                                                                                   ggatCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                      >  2:780895/1‑149 (MQ=255)                                                                                                                                     atCAAC‑CAGCGGCTTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgc                   <  2:756711/150‑1 (MQ=255)                                                                                                                                      tCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc                 >  2:26520/1‑151 (MQ=255)                                                                                                                                       cAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc                 >  2:631445/1‑150 (MQ=255)                                                                                                                                        aaC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                >  2:339974/1‑150 (MQ=255)                                                                                                                                        aaC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                >  2:529846/1‑150 (MQ=255)                                                                                                                                        aaC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                >  2:1010052/1‑150 (MQ=255)                                                                                                                                           ccgggggccCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAAt              <  1:907171/145‑1 (MQ=255)                                                                                                                                                 gcccTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg       <  2:467583/148‑1 (MQ=255)                                                                                                                                                      tACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCg  >  2:326090/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTAC‑GGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAAC‑CAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907454‑907751

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