Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	809,978	T→G	5.1%	L212R (CTG→CGG)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	809,978	0	T	G	5.1%	249.9 / 4.9	79	L212R (CTG→CGG)	bioB	biotin synthase
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base G (2/2);  total (42/37)
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.49e-01

ATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCC  >  NC_000913/809831‑810124                                                                                                                                                    |                                                                                                                                                   atatCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGAAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGc                                                                                                                                                  >  2:1642135/1‑150 (MQ=255)   atcatcACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCg                                                                                                                                                >  2:128668/1‑150 (MQ=255)      atcaCCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGAcg                                                                                                                                             >  1:13935/1‑150 (MQ=255)        caccacACGCACTTATCAGGACCGCCTCGATACGCTGG‑AAAAGTCCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGAcgccg                                                                                                                                          >  2:1726641/1‑150 (MQ=255)         accacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGAcgccg                                                                                                                                          <  2:948505/150‑1 (MQ=255)         accacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGAcgccg                                                                                                                                          >  1:1727059/1‑150 (MQ=255)         accacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGAcgccg                                                                                                                                          >  2:730134/1‑150 (MQ=255)         accacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCGGCCGAcgccg                                                                                                                                          >  2:361030/1‑150 (MQ=255)         accacACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCGGCCGAcgccg                                                                                                                                          >  2:368043/1‑150 (MQ=255)            acacGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGccgccg                                                                                                                                       >  1:1098350/1‑150 (MQ=255)                  aCTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGc                                                                                                                                 <  2:150715/150‑1 (MQ=255)                     tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGACGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTg                                                                                                                              >  2:1418871/1‑150 (MQ=255)                      aTCAGTAACGCCTCTATACGCTGTAAAAAGTGCGGGATGCCGGTATCAAAGTCTGTTCTGGCGGCTTTGTGGGCTTAGGCTAAACGGTAAAAGATCGCGCCGTATTATTGCTGCAACTGGCAAACCTGCCTACGCCGCCGGAAAGCTTGc                                                                                                                             <  2:13935/150‑1 (MQ=255)                       tCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGcc                                                                                                                            <  2:1098350/150‑1 (MQ=255)                       tCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCGGCCGACGCCGCCGGAAAGCGTGcc                                                                                                                            <  1:368043/150‑1 (MQ=255)                       tCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCGGCCGACGCCGCCGGAAAGCGTGcc                                                                                                                            <  1:361030/150‑1 (MQ=255)                         aGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCaa                                                                                                                          >  1:1751289/1‑150 (MQ=255)                           gAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATc                                                                                                                        <  1:730134/150‑1 (MQ=255)                            aaCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCa                                                                                                                       <  2:21393/150‑1 (MQ=255)                               gCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACa                                                                                                                    <  1:472252/150‑1 (MQ=255)                                 cTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATg                                                                                                                  >  2:1883643/1‑150 (MQ=255)                                         gCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtgaa                                                                                                          >  2:904729/1‑150 (MQ=255)                                         gCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtgaa                                                                                                          >  2:836649/1‑150 (MQ=255)                                         gCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtgaa                                                                                                          <  1:902362/150‑1 (MQ=255)                                            ggAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtgaaggt                                                                                                       >  1:1469277/1‑150 (MQ=255)                                               aaaaGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtgaaggtgaa                                                                                                    <  2:1751289/150‑1 (MQ=255)                                                 aaGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAg                                                                                                  <  2:1960516/150‑1 (MQ=255)                                                 aaGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAg                                                                                                  <  1:1732217/150‑1 (MQ=255)                                                        cgATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGcc                                                                                                                            >  2:1068193/1‑117 (MQ=255)                                                        cgATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGcc                                                                                                                            <  1:1068193/117‑1 (MQ=255)                                                            gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCtt                                                                                       >  1:1782656/1‑150 (MQ=255)                                                             ccGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTg                                                                                      >  2:87551/1‑150 (MQ=255)                                                             ccGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGCCGAAACGGTAAAAGATCGCCCCGGATTATTGCTGCAACTGGCAAACCTCCCGACGCCGCCGGAAAGCGTGCCCATCAACATCCTCGTGAAGGTGAAAGGCACGCCGCTTg                                                                                      >  2:1245464/1‑150 (MQ=255)                                                              cGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGc                                                                                     >  2:643399/1‑150 (MQ=255)                                                               gggATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGcc                                                                                    <  1:836649/150‑1 (MQ=255)                                                               gggATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGcc                                                                                    <  1:1726641/150‑1 (MQ=255)                                                                ggATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCg                                                                                   <  1:1642135/150‑1 (MQ=255)                                                                    cAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATaa                                                                               >  1:475002/1‑150 (MQ=255)                                                                     aaaGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAAc                                                                              <  1:1491874/150‑1 (MQ=255)                                                                     aaaGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATAGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAAc                                                                              >  1:88686/1‑150 (MQ=255)                                                                         tCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACgatg                                                                          >  1:283304/1‑150 (MQ=255)                                                                           tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACgatgat                                                                        >  1:230809/1‑150 (MQ=255)                                                                              tCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTc                                                                     <  1:1147515/150‑1 (MQ=255)                                                                                tGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGa                                                                   <  1:1557281/150‑1 (MQ=255)                                                                                 ggcggcATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtg                                                                                                            <  1:1730927/108‑1 (MQ=255)                                                                                 ggcggcATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTggtg                                                                                                            >  2:1730927/1‑108 (MQ=255)                                                                                            gggCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGAttt                                                       >  1:231005/1‑150 (MQ=255)                                                                                            gggCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGAttt                                                       >  2:841483/1‑150 (MQ=255)                                                                                            gggCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGAttt                                                       >  2:1831968/1‑150 (MQ=255)                                                                                              gCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTa                                                     <  1:904729/150‑1 (MQ=255)                                                                                              gCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTa                                                     <  1:1418871/150‑1 (MQ=255)                                                                                                 tAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTc                                                  >  2:1845728/1‑150 (MQ=255)                                                                                                    gCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCa                                               <  2:283304/150‑1 (MQ=255)                                                                                                    gCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCa                                               >  1:700731/1‑150 (MQ=255)                                                                                                     cGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCAc                                              >  1:1851592/1‑150 (MQ=255)                                                                                                          cGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTg                                         >  1:1585513/1‑150 (MQ=255)                                                                                                          cGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTg                                         >  2:1126665/1‑150 (MQ=255)                                                                                                            gTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGc                                                <  2:416328/141‑1 (MQ=255)                                                                                                            gTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGc                                                >  1:416328/1‑141 (MQ=255)                                                                                                             tAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTTAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCgg                                      <  2:412427/150‑1 (MQ=255)                                                                                                             tAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCgg                                      <  2:1795502/150‑1 (MQ=255)                                                                                                               aaaGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAACGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTc                                    >  1:1657756/1‑150 (MQ=255)                                                                                                                 aGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTcgc                                  >  1:1338777/1‑150 (MQ=255)                                                                                                                    tCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTcgcgcg                               <  2:1782656/150‑1 (MQ=255)                                                                                                                     cgcgCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCgg                              <  2:1851592/150‑1 (MQ=255)                                                                                                                         ccGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCa                          <  1:1730700/150‑1 (MQ=255)                                                                                                                         ccGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCa                          <  2:552323/150‑1 (MQ=255)                                                                                                                         ccGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCa                          >  2:1493620/1‑150 (MQ=255)                                                                                                                           ggATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCatg                        <  2:1727059/150‑1 (MQ=255)                                                                                                                            gATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCatga                       <  2:475002/150‑1 (MQ=255)                                                                                                                                attGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTACGGTCGCGCGGATCATGATGcc                   >  1:1023362/1‑150 (MQ=255)                                                                                                                                   gctgcAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAAc                >  1:1255794/1‑150 (MQ=255)                                                                                                                                   gctgcAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAAc                >  2:1122636/1‑150 (MQ=255)                                                                                                                                         aCTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTa          <  1:128668/150‑1 (MQ=255)                                                                                                                                         aCTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTa          <  1:87551/150‑1 (MQ=255)                                                                                                                                          cTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTAc         >  1:1632154/1‑150 (MQ=255)                                                                                                                                              cAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTgc     >  1:740714/1‑150 (MQ=255)                                                                                                                                               aaaCCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTgcg    <  2:26044/150‑1 (MQ=255)                                                                                                                                                 aCCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGcc  <  1:1831968/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ATATCATCACCACACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCC  >  NC_000913/809831‑810124

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