Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,526	0	C	A	45.5%	111.3 / 53.8	100	G594C (GGC→TGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (38/16);  new base A (0/45);  total (38/62)
Fisher's exact test for biased strand distribution p-value = 6.36e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTA  >  NC_000913/364381‑364674                                                                                                                                                   |                                                                                                                                                      tCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑t                                                                                                                                                 >  1:373855/1‑151 (MQ=255)  cGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGCCGTATCGCCAAAA‑t                                                                                                                                                 >  1:572009/1‑150 (MQ=255)        tCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑CCAccgc                                                                                                                                           <  2:199861/150‑1 (MQ=255)           cTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTa                                                                                                                                       >  1:1060571/1‑151 (MQ=255)                   gTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGAc                                                                                                                                >  2:399755/1‑150 (MQ=255)                     ttGCCCGGAAAACCGGAACGGGAAAAACTGCTGCGGGTGTTTTGCTCCCGTCAGCGCTGGATGCGGCGTGCGGTGGGC‑AAAGACCAGCCCTTTCATACAGACCGGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCAc                                                                                                                             <  1:485377/151‑1 (MQ=255)                            gATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTg                                                                                                                       >  2:1402167/1‑150 (MQ=255)                            gATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTg                                                                                                                       >  2:269129/1‑150 (MQ=255)                            gATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTg                                                                                                                       >  2:334189/1‑150 (MQ=255)                            gATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGCCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGGTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGtcg                                                                                                                       >  2:543453/1‑148 (MQ=255)                                    gAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtca                                                                                                              >  1:1590116/1‑151 (MQ=255)                                       cTGAAAAATCTGCTGCTGGTGTTTTGCTTCCGTCCGCGGGGGATGCGGGGGGGGGCGGGA‑AAAGCCCAGCCCGTCCATAAAGAACGGGGGATGGTTCGGCGTCCCCCCAAAA‑TCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTtcatca                                                                                                           <  2:1378326/151‑1 (MQ=255)                                         ggAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGCATCGCCAACA‑TCACCGCCGTAAGCCGACCACGGGTTTCCGTTTTCATCata                                                                                                         >  1:877199/1‑151 (MQ=255)                                           aaaacCTGCTGCTGGTGTTTTGCTCCCGCAAGCGTGGGATGGGGCGTGCGGTGGAC‑AAGGACCAGACCGTTCATACACAACTGGATATCGTTCGGCGTATCTCCAAAA‑ACGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAtt                                                                                                       <  2:1363538/151‑1 (MQ=255)                                                ctgctgctgGTGTTTTGCTCCCGTCAGCGCGGGATGCGGCGTGCGGTCGGC‑AAAGACCAGCCCGTTCATACAGACCTGGCGATCGTTCGGCTTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                  <  1:269040/151‑1 (MQ=255)                                                ctgcgggtgGTTTTTTGCTTCCGTCAGCGCTGGAGGCGGCGTGCGGTCGAC‑AAAGCCCAGCCCGTTCATACAGAACTGGCGTTCGTTCGGCGTATCGCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                  <  1:1339546/151‑1 (MQ=255)                                                cggttggtgtttttttGCTCCTGTCAGCGCAGGAGGGGGGGTGCGG‑CGGCAAAAGCCCAACCTGTACATAAAAAACGGGCGATGGTTGGGTGTTTCCCAAAAA‑TCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                  <  1:842315/139‑1 (MQ=255)                                                cgcctgctgTTTTTTTGCTTCCGCCAGCGCTGGAGGCGGCGTGGGGTCGGC‑AAAGACCAGCCCGTTCATAAAGAACGGGCGCTCTTTCGGCGTTTCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATc                                                                                                  <  1:824913/148‑1 (MQ=255)                                                 gacgactgGTGTTTTGCTCCCGTCCGGGCTGGTGGCGGGGTGGGGTCGGA‑AAAGCCCAGGCCGTTAATACAGAGCGGGCGTTCGTCAGGGGTACCCCAAAAA‑CCACCGCCGTAAGCCGACCGCGGGTTGCCGTTTTCATCATATTTAATc                                                                                                  <  1:458239/145‑1 (MQ=255)                                                      ctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGa                                                                                             >  1:1280071/1‑150 (MQ=255)                                                        ggTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGa                                                                                             >  2:1041095/1‑148 (MQ=255)                                                         gtgtTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAC‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGa                                                                                         >  2:990532/1‑151 (MQ=255)                                                         gtgtTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGa                                                                                         >  2:1204502/1‑151 (MQ=255)                                                          ttttttGCTCCCGACGGCGCGGGGTGGGGCGTGCGGTGGGA‑AAAGCCAAGCCCTTTCATACAGACCGGGCGATCTTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGAt                                                                                        <  1:1316865/149‑1 (MQ=255)                                                           gtTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCTTACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGAt                                                                                        >  2:859696/1‑150 (MQ=255)                                                            ttttGCTTCCGTCAGCGCGGGATGGGGGGTGCGGTCGGA‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATc                                                                                       <  2:6852/150‑1 (MQ=255)                                                            ttttGCTCCCGTCAGCGGTGGATGCGGCGTGGGGTCGGC‑AAAGACCAGCCCTTTCATAAAGAACGGGCGATCGTTCGGGTTTCCGCCAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATc                                                                                       <  1:504068/150‑1 (MQ=255)                                                               tGCTTCCGTCAGGGCTGGAGGGGGCGTGGGGTCGGA‑AAAGACCAGACCGTTAATACAAAACGGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAcc                                                                                   <  2:973726/151‑1 (MQ=255)                                                                tctcccGTCAGCGCGGGATGGGGCGTGGGGTCGGA‑AAAGCCCAGCCCGTTCAAAAAGACCGGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  2:296288/147‑1 (MQ=255)                                                                gCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAcc                                                                                   >  2:1609835/1‑150 (MQ=255)                                                                cccccggccAGCGCGGGATGCGGGGGGCGTGCGGC‑AAAGACCAGCCCGTTCATACAGACCGGGCGATTGTTCGGGGTTCCCCAAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:750129/143‑1 (MQ=255)                                                                ccTTCGGTCAGCGCTGGATGGGGCGTGCGGTCGGC‑AAAGACCAGCCCGTTCATCAAGAGCGGGCGATCGTTCGGCTTTTCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:1447589/150‑1 (MQ=255)                                                                ccTTCCGTCGGCGTGGGATGCGGCGGCCGGTCGCA‑AAAGACCATCCCTTCCAAACAGACCGGGCGCTGTTTCGG‑GTTTCCCCAAAATTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:660997/150‑1 (MQ=255)                                                                 ctcccGTCGGGGCGGGATGCGGCGTGGGGTCGGA‑AAAGACCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGacccaccc                                                                                  <  2:52881/147‑1 (MQ=255)                                                                 ctcccGTCAGCGGTGGATGGGGCGTGGGGTCGGA‑AAAGACAAGACCGTTCAGAAAGACCGGGCGATGGTTGGGGGTATCGCAAAAA‑ACACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:461665/147‑1 (MQ=255)                                                                 ctcccGTCAGCGCTGGATGCGGCGTGCGGTCGGA‑AAAGCCCACCCCGTTCATACAGAACTGGCGCTCGTTGGGGGTTCCCCCAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:789947/147‑1 (MQ=255)                                                                 cTTCCGTCAGCGCTGGATGCGGGGTGCGGTCGGC‑AAAGACCAGACCGTTCAAACAGAACTGGCGATCTTTCGGCGTACCCCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  2:1196933/150‑1 (MQ=255)                                                                 cTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTACCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:105115/150‑1 (MQ=255)                                                                 cTTCCGTCAGCGCTGGAGGCGGCGTGCGGTGGGC‑AAAGACAAGACCGTTCATACAGACCGGGCGATCGTTCGGGTTACCCCCAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCAccc                                                                                  <  1:583948/150‑1 (MQ=255)                                                                     cttCAGCGCGGGGGGGGGGGGGCGGTCGGA‑AAAGCCCAGCCCGTTCATACAGACCTGGCGGTCGTTGGGAAAATCGCAAAAC‑CCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTcc                                                                             <  2:1169498/149‑1 (MQ=255)                                                                     cGTCAGCGCGGGATGCGGCGTGCGGTCGCA‑AAAGACCAGCCCGTTCATACAGAACTGGCGATCGTTCGGCGTATCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTcc                                                                             <  2:121933/151‑1 (MQ=255)                                                                          gcgcTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGAAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGa                                                                         >  2:861956/1‑150 (MQ=255)                                                                           cgttGGATGCGGCGTGCGGTCGGC‑AAAGCCCAACCCGTTCATAAAGAACGGGCGTTGGTTCGGCGTACCGCAAAAA‑TCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  2:1046178/147‑1 (MQ=255)                                                                           cggtgagcggggcgTGCGGCCGGA‑AAAGCCCGGCCCGTCCAAAGAAAACGGGCGATGTTTGGGCGTATCACAAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCCCCCAGTCCCAGAc                                                                        <  1:1122212/140‑1 (MQ=255)                                                                           cggggggtggggcgTGCGGTGGGC‑AAAGACCGGCCCGTTCATACAGACCGGGCGATCGTGCGGGGTTTCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACg                                                                       <  1:62699/141‑1 (MQ=255)                                                                           cgcggggcggggcgTGGGGGGGGC‑AAAGCCGAGACCGTCCAAACAAACCGGGGGCTGGTTCGGCGTATCGCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  2:957497/140‑1 (MQ=255)                                                                           cgcgggagggggggggCGGTCGGA‑AAAGCCAAGCCCGTTAAAACAGACCGGGCGATCGTTCGGCGTATCCCAAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  1:1135717/135‑1 (MQ=255)                                                                           cgcgGGAGGGGGGGTGGGGCCGGA‑AAAGCCAAGCCCGTACATACAAAACGGGCGATGGTGCGGGGTTTCCCAAAAA‑TAACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  1:398102/150‑1 (MQ=255)                                                                           cgcgGGAGGGGGGGTGCGGTCGGA‑AAAGACCAGACCGTTCATAAAGACCGGGCGTTCGTTCGGCGTCTCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  1:766087/150‑1 (MQ=255)                                                                           cgcTGGGTGGGGCGTGCGGGCGGA‑AAAGCCCGGCCCGTTCATACAAAACTGGCGATCGTTCGGCGTACCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  2:505925/150‑1 (MQ=255)                                                                           cgcTGGATGGGGCGTGCGGTCGGA‑AAAGACCAGCCGGTTAATACAGAACGGGCGATCGTTCGGCGTATCGCAAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                        <  1:590075/150‑1 (MQ=255)                                                                              tatgtGCGGCGTGCGTTCGGA‑AAAGCCCAGCCCGTTCATACAGATCTGGCATTCTTTCGGCTTATCGCCAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAg                                                                    <  1:82812/147‑1 (MQ=255)                                                                               ggATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAg                                                                    >  2:1330060/1‑150 (MQ=255)                                                                                    cggcggGCGGCCGGA‑AAAGCCCAGCCCGTTCATACAAAACGGGCGTTCGTTGGGCGTACCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                              <  2:1454634/151‑1 (MQ=255)                                                                                     ggcgTGGGGGCGGA‑AAAGCCCAGACCGTTCATACAGAACGGGCGATCGTTGGGCTTATCCCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                              <  2:176419/150‑1 (MQ=255)                                                                                           cGGTCGGA‑AAAGACCAGCCCGTTCATAAAGAACTGGCGATCGTTGGGCGTACCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAAc                                                       <  2:1008904/151‑1 (MQ=255)                                                                                            ggTCGGA‑AAAGACCAGACCTTTCATAAAAACCTGGCGATCGTTCGGCGTATCGCAAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAAc                                                       <  2:1559011/150‑1 (MQ=255)                                                                                             gTCGGA‑AAAGCCCAGCCCGTTCATACAAACCGGGCGATCGTTCGGCGTATCGCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                                      <  2:913366/150‑1 (MQ=255)                                                                                             gTCGAC‑AAAGCCCAGACCGTCCATACAGACCTGGCGCTCGTTCGGCTTATCGCCAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACg                                                      <  1:700691/150‑1 (MQ=255)                                                                                              cgggC‑AAAGCCCAGCCCGTTCATAAAGATCTGGCGATCTTTCGGCTTATCGCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                    <  2:821178/149‑1 (MQ=255)                                                                                                 gC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                 >  2:452634/1‑151 (MQ=255)                                                                                                 gC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                 >  2:1725039/1‑151 (MQ=255)                                                                                                 gC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACAAAGCCGCCCTGTAAACGGGGAt                                                 >  2:861620/1‑151 (MQ=255)                                                                                                   aaaaGCCCAGCCCGTTCATACAGAACGGGCGATGGTGGGGCTTATCGAAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAt                                                 <  1:269332/149‑1 (MQ=255)                                                                                                    aaagccacgcccGTTAAAACAGAACGGGAGACCGTCGGGCGTAGCGCGAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATAc                                               <  1:1254404/141‑1 (MQ=255)                                                                                                    aaaGCCCAGCCCTTTCATACAGAACGGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATAc                                               <  2:327134/151‑1 (MQ=255)                                                                                                    aaaGACCAGACCTTTCATACAGACCGGGCGATCGTTCGGCGTTTCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATAc                                               <  1:583733/151‑1 (MQ=255)                                                                                                      agaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                                >  2:98080/1‑148 (MQ=255)                                                                                                      agaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                                >  1:331755/1‑148 (MQ=255)                                                                                                      agaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                                >  2:893234/1‑148 (MQ=255)                                                                                                      agaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATa                                                >  1:1472325/1‑148 (MQ=255)                                                                                                       gaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGa                                            >  2:1043908/1‑151 (MQ=255)                                                                                                       gaccagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGa                                            >  2:1082027/1‑151 (MQ=255)                                                                                                        accagaccGTTCATACAGACCGGGCGATGGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTg                                             <  1:1154271/149‑1 (MQ=255)                                                                                                        accagaccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTg                                             <  1:508820/149‑1 (MQ=255)                                                                                                              ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                     >  2:1459848/1‑151 (MQ=255)                                                                                                              ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                     >  2:1462568/1‑151 (MQ=255)                                                                                                              ccGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                     >  2:1646088/1‑151 (MQ=255)                                                                                                                 ttCATACAAAACTGGCGATCGTTCGGCGTACCCCAAAAT‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGcc                                   <  1:231599/150‑1 (MQ=255)                                                                                                                    aTACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                >  2:25926/1‑150 (MQ=255)                                                                                                                    aTACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                >  1:1417761/1‑150 (MQ=255)                                                                                                                     acaaaaccgggCGATCGTTGGGCGTATCCCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCa                              <  2:1663187/142‑1 (MQ=255)                                                                                                                         gAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          >  2:996880/1‑151 (MQ=255)                                                                                                                         gAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          >  2:1626802/1‑151 (MQ=255)                                                                                                                         gAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          >  1:1137956/1‑151 (MQ=255)                                                                                                                         gAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATACTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          >  2:1091490/1‑151 (MQ=255)                                                                                                                         aacctgcgaatggTTCGGCGTATCCCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          <  1:54160/139‑1 (MQ=255)                                                                                                                         aaccgggcgtgggtgcGGCGTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          <  1:515014/136‑1 (MQ=255)                                                                                                                         aaccgggCGGTCGTGCGGCGTATACCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCAGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          <  1:269964/146‑1 (MQ=255)                                                                                                                         aaaCGGGCGATCGTGCGGCTTATCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAt                          <  2:85863/150‑1 (MQ=255)                                                                                                                             tGGCGATGGTTCGGGGTATCCCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                        <  2:1665407/149‑1 (MQ=255)                                                                                                                              ggCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAg                      >  2:204059/1‑150 (MQ=255)                                                                                                                               gggATCTTTGGGCGTATCGCAAAAA‑TCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                     <  2:1351583/148‑1 (MQ=255)                                                                                                                               gCGATCGTTCGGGGTATCCCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                     <  2:20132/150‑1 (MQ=255)                                                                                                                               gCGATCGTGGGGCTTATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                     <  1:614531/150‑1 (MQ=255)                                                                                                                               gCGACCGTTCGGCGTATCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                     <  2:57755/150‑1 (MQ=255)                                                                                                                                cgatggtgggggttATCGCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                    <  2:1731963/137‑1 (MQ=255)                                                                                                                                cGATCGTGCGGCTTACCCCAAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                    <  1:185943/150‑1 (MQ=255)                                                                                                                                       tgggCGTATCCCAAAAA‑CCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGcc            <  2:1445597/149‑1 (MQ=255)                                                                                                                                                 gCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTa  >  2:1270849/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGC‑AAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAA‑TCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTA  >  NC_000913/364381‑364674

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