Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,686,021	0	A	C	45.1%	135.9 / 15.0	91	G1051G (GGT→GGG)	bcsC	cellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base A (30/20);  new base C (0/41);  total (30/61)
Fisher's exact test for biased strand distribution p-value = 6.88e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTC  >  NC_000913/3685878‑3686165                                                                                                                                                 |                                                                                                                                                    cgcgCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCt                                                                                                                                              <  2:553494/151‑1 (MQ=255)  gcgcCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCt                                                                                                                                              >  1:1434732/1‑150 (MQ=255)     ccGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCa                                                                                                                                          <  2:204629/151‑1 (MQ=255)      cGTGTAGCCGAATCCCTGGCTGCTGCCCCCGTCCTTGGATTGCCGCGCAGCTTCTTCCTGACAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGTCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTcccg                                                                                                                                          <  1:1555785/150‑2 (MQ=255)       gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        >  2:918896/1‑151 (MQ=255)       gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        >  2:362482/1‑151 (MQ=255)       gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        >  2:555296/1‑151 (MQ=255)       gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACAACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        >  1:88699/1‑151 (MQ=255)        tgtAGCCAAAGCCCTGGCTGCTGCCGCCAACGTTGGTTTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGGCCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        <  2:913325/150‑1 (MQ=255)            gCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                   >  1:172478/1‑151 (MQ=255)             ccGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                   <  2:1355204/150‑1 (MQ=255)             ccGAAGCCCTGACTGCTGCCCCCATCGTGGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCTACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                   <  1:379758/150‑1 (MQ=255)                aaGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑ttt                                                                                                                               <  2:796690/150‑1 (MQ=255)                  gCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑c                                                                                                                            >  1:77632/1‑150 (MQ=255)                  gCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑c                                                                                                                            >  1:1196476/1‑150 (MQ=255)                   cccTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑c                                                                                                                            <  1:406895/149‑1 (MQ=255)                          tgctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTc                                                                                                                  <  1:280238/151‑1 (MQ=255)                           gctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGc                                                                                                                                               <  2:116838/124‑1 (MQ=255)                           gctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGc                                                                                                                                               >  1:116804/1‑124 (MQ=255)                                     tCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCa                                                                                                       >  2:1365774/1‑151 (MQ=255)                                      ctttGGTTTGCCGCGCAGTTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGTTTTTGGTGCGTG‑AATGCGCCCCCGAGCCAGACGCCCCCAGCCCCCACGCCCAT‑TTTT‑CCGTG‑CGCCCCCGCCACATCAc                                                                                                      <  1:955502/149‑1 (MQ=255)                                      cGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCa                                                                                                       >  2:1061377/1‑150 (MQ=255)                                         tGGTTTGCCGCCCGCCTTTTCCCTGCCGGTCGGTCGGGATAAAATCCATCAGGGGAAAACGGGGCATGTTTTTGGTGCGTG‑AATGCGCCCACGAGCCAGACGCCCCCAGCCCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCgg                                                                                                   <  2:1280475/151‑1 (MQ=255)                                         tGGATGGCCGCGCAGCTTCTTCCTGCCAGTCGGTGGGGATCAGATTCATCAGCGGAAAACGCGGCATGTTTTTGGGGCGTG‑AATGCGCCCACGAGCCAGCCGCCCCCAGCCCCCCCGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCgg                                                                                                   <  1:1380794/151‑1 (MQ=255)                                            tttcccgcccggtttctccCTGCCAGTCGGTGGGGATAAGATTCATCAGCGGATAGCGCGGCATTGTTTTGGGGCGTA‑AATGCGACCACGAGCCAGACGCCCCCAGCCCCCCCGACCGG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTa                                                                                                 <  2:88728/137‑1 (MQ=255)                                            ttcgCCGCGCATTTTCTCCCGGCCGGTCGGCGGGGATCAGATTCATCAGCGAATAACGCGGCATGTTTTTGGTGCGTG‑ATTGCGCCCACGACCCAGCCCCCCCCAGCCCCCCCGCCCAG‑TTTT‑CGGTG‑CGCTCCCGCCACATCACCGGTAt                                                                                                <  1:969266/148‑1 (MQ=255)                                             ttGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATg                                                                                               >  1:68026/1‑151 (MQ=255)                                              tgcgccgcattttctccCCGCCAGTCGGTAGGGATCAGATTCATCAGCGGATACCGCGGCATGGTTTTGGGGCGGG‑AAGGCGCCCACGACCCAGCCGCCCCCTCCCCCCACGACCAT‑TTTT‑CCGGG‑CGCCCCCGCCACATCACCGGTATgg                                                                                              <  1:1674530/140‑1 (MQ=255)                                              tGCCGCCCAGCTTTTCCCTGCGGGTCGGGCGGGATCAAATTCAACAGCGGAAAGCGCGCCGTGGTTTTGGGGCGGG‑AATGCGACCACGGGCCAGCCCCCCCCAGCTCCCGCGACCAT‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATgg                                                                                              <  1:40166/151‑1 (MQ=255)                                                   cgcattttctccCTCCCAGTCGGTCGGAATCAGATTCATCAGCGAATACCGCGGCATGGTTTTGGTGGGTG‑AATGCGCCCCCGACCCAGACCCCCCCAGCCCCCCCACCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAAc                                                                                         <  1:1548749/145‑1 (MQ=255)                                                       gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:1285119/1‑151 (MQ=255)                                                       gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:322222/1‑151 (MQ=255)                                                       gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:1363349/1‑151 (MQ=255)                                                       gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:665905/1‑151 (MQ=255)                                                        tttcttcCTGCCGGTCGGTGGTGACGAGATCCGGGGGGGGATAGCGAGGCATGTTTTTGGTGGGTA‑AAGGCAACCACGAGCCAGACGCACCCGGCACCCACAACAAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGAc                                                                                      <  1:813147/148‑1 (MQ=255)                                                           ttcccTGCCGGTGGGTCGGGATCAAATTCACCGGGGGAAACCGCGGCAGGTTTTGGGGGCGTG‑AATGCGCCCACGCCCCGGCCGCCCCCGGCCCCCCCGACCTT‑TTTT‑CGGTC‑CCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTa                                                                                 <  1:348902/148‑1 (MQ=255)                                                            tccccccgaGTGGGTGGGGATAAGATTCACCAGGGGAAAACGCGGCATTTTTTTGGGGCGAA‑AATGCGACCACAAGCCAGGCGCACCCAGCTCCCCCACCCAT‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTAt                                                                                <  2:442512/143‑1 (MQ=255)                                                               cggccggtggggggggATAAGATCCACCGGGGAATACCGCGGCGGGGTTTGGGGGGGGG‑AAGGAGCCCACGGGCAAAACGCCCCCAGCCCCCCCG‑CCAGTTTTT‑CCGTG‑CCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                             <  2:1129137/138‑1 (MQ=255)                                                               cgcgcAGTGGGTCGGAATCAGATTCATCAGCGAATACCGCGGCTTGTTTTTGGGGCGTG‑AACGCGCCCACGAGCCAGCCGCACCCAGCCCCCACAACCTT‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                             <  1:1568246/147‑1 (MQ=255)                                                               cTGCCAGTGGGTCGGGATCAGATTCACCAGCGAAAACCGCGGCAGGTTTTTGGGGCGTA‑AAGGCGCCCCCGAGCCAGACGCCCCCGGCCCCCACGCCCAT‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                             <  2:1061395/151‑1 (MQ=255)                                                                  ccggtgggtgggGAGCAGATTCATCAGCGGAAACCGCGGCTTGTTTTGGGGGCGTG‑AATGCGCCCACGAGCCGAACCCCCCCAGCCCCCACGCCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTg                                                                           <  1:1496024/144‑1 (MQ=255)                                                                      tcgggggggATCAGATCCATCAGCGGATACCGCGGCATGTTTTGGGTGCGTG‑AATGCGCCCACGAGCCAGACGCACCCAGCCCCCACGACCAG‑TTTT‑CCGGG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGa                                                                       <  1:402986/144‑1 (MQ=255)                                                                        gggcggGACAAGATTCAGCGCCGGAAAACGCGGCATGTTTTGGGGGGGTG‑AATGCGCCCCCGAGCCAGACGCCCCCACCCCCCACGCCCAT‑TTTT‑CCGGG‑CGCGCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACt                                                                     <  1:845578/147‑1 (MQ=255)                                                                         gtcggGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGa                                                                       >  1:11348/1‑147 (MQ=255)                                                                            gggacaaGATACACCGGCGAATACCGCGGCGGGGTTTGGGTGCGTG‑AATGCGACCCCGAGCCAGACGCACCCGGCCCCCACAACCGG‑TTTT‑CCGGG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtag                                                                 <  1:1356456/144‑1 (MQ=255)                                                                            gggATAAGATTCATCAGGGGAAACCGGGGCTTGTTTTTGGTGGGGG‑AAGGCGCCCACGAGCCAGACGCCCCCAGCCCCCCCGACCAG‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtag                                                                 <  2:114913/150‑1 (MQ=255)                                                                             ggATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtagta                                                               >  1:731225/1‑151 (MQ=255)                                                                             ggACCAAATTCCCCAGGGGAAAACGGGGCATGGTTTTGGGGCGGG‑AATGCGACCACGAGCCAGACCCACCCACCTCCCCCGACCAT‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtagt                                                                <  2:1706749/150‑1 (MQ=255)                                                                               acaaaactcataaGCGGAAACGGCGTCATGTTTTGGGGGGTTG‑AGTCCGACCACGGCCAAGCCCCCCCCAGCCCCCACACCCAT‑TTTC‑CCGTC‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                             <  1:108052/144‑1 (MQ=255)                                                                               aacAGATTCACCGCGGGAAAAGGCGGCATGG‑TTTGGGGGGTGAAATGCGACCGCGCCCCAGACGCCCCCAGCCCCCACGACCGG‑TTTC‑CGGGC‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                             <  1:450279/149‑1 (MQ=255)                                                                               aacAGATACACCGGCGGAAAACGCGGCTTGGTTTTGGGGCGTG‑AAGGCGACCACGACCCAGCCCCCCCCCGCCCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                             <  2:836966/149‑1 (MQ=255)                                                                               aTCAGATTCACCAGCGAAACACGCGGCGTGTTTTTGGGGCGTG‑ACTGCGCCCACGAGCCGGCCGCCCCCGCCCCCCACGACCAT‑TTTT‑CCTGG‑CCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                             <  1:1305090/151‑1 (MQ=255)                                                                                taaGATCCACCGGCGGAAACCGCGGCATGTTTTTGGTCCGGA‑AATGCCCCCACGACCCAGCCGCCCCCCGCCCCCACGACCAT‑TTTC‑CGGCT‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                            <  1:1364377/149‑1 (MQ=255)                                                                                tCAGTTTCATCGGCGGATACCGCGGCTTGGTTTTGGTGCGGA‑AATGCGCCCACGGCCCAGCCGCCCCCAGCCCCCCCGCCCAT‑TTTT‑CCGTG‑CCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                            <  2:280687/151‑1 (MQ=255)                                                                                tCAGATTAACCAGCGAATAACGCGGCGTGTTTTTGGTGTGTG‑AAGGCGACCCCGAGCAAGACGCACCAGGCCCCCCCGCCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                            <  1:512608/151‑1 (MQ=255)                                                                                tCAGATCCACCAGCGGATACCGCGCCATGTTTTTGGTGCGTG‑AATGCGACCACGAGCAAGACCCCCCCACCTCCCCCCCCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                            <  2:1628079/151‑1 (MQ=255)                                                                                aaaGATTCACCAGCGAAAACCGGGGCTTGTTTTGGGTGGGTG‑AAGGCGCCCACGGCCCAGACGCCCCCAGCCCCCACGCCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                            <  1:1283476/149‑1 (MQ=255)                                                                                    tttcaccgggggAAAACGCGGCATGTTTTGGGTGGGTA‑AAGGCGCCCCCGAGCCAGACGCACCCGGCCCCCACACCCAT‑TTTT‑CCGGG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGccc                                                        <  2:1609024/141‑1 (MQ=255)                                                                                    tttCACCAGGGGAAAACGGGGCATGTTTTTGGTGCGTG‑AAGGCGACCCCGAGCCAGACCCCCCCAGCTCCCCCGACCAT‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAccc                                                           <  2:571026/147‑2 (MQ=255)                                                                                        ctgaGCGAAAAACGCGGCAGGTTTTTGGTGCGTA‑AATGCGCCCACGCCCCAGCCGCCCCCAGCCCCCCCGCCCGT‑TTTC‑CCGTC‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAc                                                    <  1:834159/148‑1 (MQ=255)                                                                                        cccgggggAAAACGGGGCATGTTTTTGGTGGGTG‑AAGGCGCCCACGGGCAAGACGCCCCCAGCCCCCACGACCAT‑TTTC‑CCGGG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAc                                                    <  2:1360295/145‑1 (MQ=255)                                                                                         tcagggaaaaccgcgGCGTGGTTTTGGGGGGGA‑AATGCGCCCACGACCGAGAGCCCCCCAGCCCCCACGCCCAT‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAcc                                                   <  2:1467437/140‑1 (MQ=255)                                                                                         ccaGCGGAAACCGGGGCATGGTTTTGGTGCGTA‑AATGCGACCCCGAGCCAGACGCCCCCAGCCCCCACGACCAG‑TTTT‑CCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAcc                                                   <  2:187801/150‑1 (MQ=255)                                                                                            gCGAATACCGCGCCTTGTTTTTGGTGGGTG‑GAGGCGCCCACGAGCCAGACCCCCCAACCCCCCCCGCCCA‑‑TTTC‑CCGGG‑GGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgag                                                <  1:39521/150‑1 (MQ=255)                                                                                             cGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagt                                               >  2:1628035/1‑151 (MQ=255)                                                                                              ggATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                              >  1:1728244/1‑151 (MQ=255)                                                                                               gaaaccgcgCCTTGTTTTTGGTGCGTA‑AATGCGCCCCCGCCCCAGCCGCCCCCACCCCCCCCGCCCAT‑TTTT‑CCGTG‑CCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                              <  2:489204/145‑1 (MQ=255)                                                                                               gaaaaCGCGGCATGTTTTTGGTGGGGG‑AAGGCGCCCACGAGCCAGCCGCCCCCAGCCCCCACGACCAT‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                              <  2:344814/147‑1 (MQ=255)                                                                                               gATAACGGGGCGTGGTTTGGGTGCGGA‑AATGCGACCACGGCCCAGACGCCCCCGGCCCCCCCGACCAG‑TTTT‑CCGGG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                              <  1:1353330/150‑1 (MQ=255)                                                                                               gATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtga                                             >  2:1253151/1‑151 (MQ=255)                                                                                                   aCGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAg                                         >  2:864074/1‑151 (MQ=255)                                                                                                    cgcgGTGTGTTTTTGGGGCGGA‑AAGGCCCCCCCGCGCCACACGCCCCCAGCCCCCCCGCCCTG‑TTTT‑CCGGG‑CCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGc                                        <  1:1015189/151‑1 (MQ=255)                                                                                                      cgGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCa                                      >  2:950698/1‑151 (MQ=255)                                                                                                             gTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:163423/1‑151 (MQ=255)                                                                                                             gTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:1594543/1‑151 (MQ=255)                                                                                                             gTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:1551350/1‑151 (MQ=255)                                                                                                             gTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:413793/1‑151 (MQ=255)                                                                                                              ttttGGTGCGTG‑AATGCGCCCACGACCAAGCCGCCCCCGGCCCCCACGACCAG‑TTTC‑CCGGG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                              <  2:914662/151‑1 (MQ=255)                                                                                                              tttgGGTGCGTG‑AATGGCCCCACGCGCCAGACCCACCCGGC‑CCCCCGGCCAT‑TTTTCCCGTG‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGagc                              <  1:1175468/151‑3 (MQ=255)                                                                                                              tttgGGGGCGTA‑AATGCCCCCACGGGCCCGCCCCCCCCAGCCCCCGGGACCTT‑TTTT‑CCGTC‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                              <  1:680621/151‑1 (MQ=255)                                                                                                              ttgtGGTGCGTG‑AATGCCCCCACGAGCCAGACCCCCCCGGCCCCCCCGCCCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                              <  1:1067711/148‑1 (MQ=255)                                                                                                               tttggggggggaaaTGCCCCCACGACCCAGCCCCCCCCAGCCCCCACGACCAG‑TTTT‑CCGGG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                              <  2:1455827/139‑1 (MQ=255)                                                                                                               tttGGGGCGGG‑AAGGCGCCCACGAGCCAGACCCCCCCGGCCCCCACGACCA‑‑TTTT‑CCCTGCCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCt                             <  1:1358736/151‑1 (MQ=255)                                                                                                                   gTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGt                         >  2:1032678/1‑151 (MQ=255)                                                                                                                     gCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCg                       >  1:1246753/1‑151 (MQ=255)                                                                                                                               cgcccACGACCCAGACGCCCCCAGCCCCCACGACCAG‑TTTC‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAg              <  1:470186/148‑1 (MQ=255)                                                                                                                                 cccccggcccAGACCCCCCCACCCCCCACGACCAT‑TTTC‑CCGTG‑CCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGAt            <  2:1153290/143‑1 (MQ=255)                                                                                                                                     cGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtg        >  1:730588/1‑151 (MQ=255)                                                                                                                                       agccagccCCCCCCACCCCCCCCGCCCAT‑TTTC‑CCGTG‑CCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtgt       <  2:89030/150‑1 (MQ=255)                                                                                                                                           cgACCCACCCCGCCCCCAGAACCAT‑TTTC‑CCGCG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc  <  1:362366/150‑1 (MQ=255)                                                                                                                                            gccGCCCCCGGCCCCCCCGACCGT‑TTTT‑CGGTC‑CGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc  <  2:1707002/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTG‑AATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAG‑TTTT‑CCGTG‑CGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTC  >  NC_000913/3685878‑3686165

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