Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,676	G→T	R87L (CGC→CTC)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,676	0	G	T	97.9%	316.4 / NA	94	R87L (CGC→CTC)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (48/44);  minor base C (0/2);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  NC_000913/3966533‑3966821                                                                                                                                                |                                                                                                                                                  gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGttt                                                                                                                                             <  2:643064/150‑1 (MQ=255) gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGttt                                                                                                                                             <  1:46476/150‑1 (MQ=255) gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGtt                                                                                                                                              <  1:928072/149‑1 (MQ=255) gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTc                                                                                                                                            >  1:1262027/1‑151 (MQ=255)  aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTc                                                                                                                                            >  1:1293733/1‑150 (MQ=255)  aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCa                                                                                                                                           <  1:1018060/151‑1 (MQ=255)  aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCa                                                                                                                                           <  1:1171974/151‑1 (MQ=255)      agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCt                                                                                                                                       >  2:561226/1‑151 (MQ=255)          cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGc                                                                                                                                   <  1:1294355/151‑1 (MQ=255)            cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCAc                                                                                                                                 <  2:168299/151‑1 (MQ=255)                gagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGt                                                                                                                             <  1:1240608/151‑1 (MQ=255)                  gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGa                                                                                                                           <  1:886136/151‑1 (MQ=255)                    ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATa                                                                                                                         <  2:718118/151‑1 (MQ=255)                       gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                      <  2:1100617/151‑1 (MQ=255)                        aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                     >  2:466678/1‑151 (MQ=255)                         aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                    >  2:1143255/1‑151 (MQ=255)                          gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   <  1:538621/151‑1 (MQ=255)                          gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                                                                   >  1:985024/1‑151 (MQ=255)                           atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtctc                                                                                                                  <  1:1152742/151‑1 (MQ=255)                            tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATctct                                                                                                                 >  1:1145701/1‑151 (MQ=255)                              tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                               >  2:428401/1‑151 (MQ=255)                                tttGGTGATGGCGTACTGGAGATATTGCAGGATGGCTTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATaaaataactggt                                                                                                              >  2:919192/1‑139 (MQ=255)                                 ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                             >  1:862182/1‑150 (MQ=255)                                 ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                             >  1:699455/1‑150 (MQ=255)                                  tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                             >  1:427763/1‑149 (MQ=255)                                    gTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                             >  1:809113/1‑147 (MQ=255)                                      gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt                                                                                                        <  2:91575/150‑1 (MQ=255)                                      gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt                                                                                                        <  2:1356339/150‑1 (MQ=255)                                        tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                     >  1:1021003/1‑151 (MQ=255)                                        tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                     >  1:1131224/1‑151 (MQ=255)                                         ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                    >  2:1067886/1‑151 (MQ=255)                                         ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                                                                                                     <  2:331154/150‑1 (MQ=255)                                           cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                                                                                                  <  1:389218/151‑1 (MQ=255)                                           cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                                                                                                  >  2:380868/1‑151 (MQ=255)                                            gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                                                                                                 <  1:53906/151‑1 (MQ=255)                                               cTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                               >  2:395815/1‑150 (MQ=255)                                                 ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                                                                            <  1:27511/151‑1 (MQ=255)                                                  gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                             >  2:540159/1‑149 (MQ=255)                                                  gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                           <  1:44186/151‑1 (MQ=255)                                                  gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                           <  1:404955/151‑1 (MQ=255)                                                      tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg                                                                                       >  2:1002410/1‑151 (MQ=255)                                                        ttGCAGGATGGTTTTGGTTTCCCCCTTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                      <  2:409586/150‑1 (MQ=255)                                                          gCAGGATGGATTTGTTTTCCTCCGTCCCGCAGACAGCTCCTCCCTCGCCGTTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCCCCGCACTGGGGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                                                                                    <  2:776621/150‑1 (MQ=255)                                                          gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                                                                                    >  1:1036657/1‑150 (MQ=255)                                                           cAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAAc                                                                                  >  2:752653/1‑151 (MQ=255)                                                            aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                 <  1:289882/151‑1 (MQ=255)                                                                 ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                            <  1:377509/151‑1 (MQ=255)                                                                  gatTTGGTTTCCTCCTTCCCGCAGACAGCTCCTCCCTCGCCGTTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACGGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                           <  2:131238/151‑1 (MQ=255)                                                                       ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                      <  1:669858/151‑1 (MQ=255)                                                                        gTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                     <  1:242832/151‑1 (MQ=255)                                                                            cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                  <  1:38647/150‑1 (MQ=255)                                                                             ctccgctccgCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                 >  1:763846/1‑150 (MQ=255)                                                                             ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                 >  1:891365/1‑150 (MQ=255)                                                                              tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                >  1:639286/1‑150 (MQ=255)                                                                               ccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                               <  1:100644/150‑1 (MQ=255)                                                                                cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                               >  1:398085/1‑149 (MQ=255)                                                                                 gTTCCGCGGACACCTCCTACCTCGCCGGTCCTGATGACATCTACGTCTCCCCGAGCAAAATCCCCCGTTTCACCCACCGCACGGGTAATCCCATCTCTGGAAAGATTCGGCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  2:593693/151‑1 (MQ=255)                                                                                 gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            >  1:505427/1‑151 (MQ=255)                                                                                  tcccGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACTTCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATCCCATCTCGGGAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  2:958359/148‑1 (MQ=255)                                                                                   tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                           >  2:1162173/1‑150 (MQ=255)                                                                                    ccGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCTTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                         <  2:427746/151‑1 (MQ=255)                                                                                     cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCTCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                        >  1:1287123/1‑151 (MQ=255)                                                                                     cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                        >  1:708470/1‑151 (MQ=255)                                                                                      gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                       >  1:890094/1‑151 (MQ=255)                                                                                        aGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                     <  2:437944/151‑1 (MQ=255)                                                                                          aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                   >  1:468450/1‑151 (MQ=255)                                                                                           cacccccacccccGCCGGCCCTGATGACATCTACGTTCCCCCCACCCAAACCCTCCTTTTCACCCTCCGCACGGGTGATCCCATTTCGGGTAGGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                  <  2:108789/139‑1 (MQ=255)                                                                                            aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                  >  1:735517/1‑150 (MQ=255)                                                                                               ccctccctcGCCGGTCCTGATGACATCTACGTTCCCCCTAGCCAAATCCTCCTTTTCAACCCCCGCACGGGTGATACCATCTCTGGAAAGATCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                              <  2:549059/150‑1 (MQ=255)                                                                                                ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                             >  2:66871/1‑151 (MQ=255)                                                                                                   cccTCGCCGGTCCTGATGACATCTACGTTCCCCCTAGCAAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                            <  2:1106674/148‑1 (MQ=255)                                                                                                   aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                            <  1:290758/149‑1 (MQ=255)                                                                                                    ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                         <  1:877996/151‑1 (MQ=255)                                                                                                        gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                     <  1:449332/151‑1 (MQ=255)                                                                                                          cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                   >  1:294731/1‑151 (MQ=255)                                                                                                              ccTGATGACATCTACGTTCCCCCTAGCCAAATCCCCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                <  2:599349/150‑1 (MQ=255)                                                                                                                tgatgaCATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                             >  2:1302966/1‑151 (MQ=255)                                                                                                                 gatgaCATCTACGTTTCCCCTAGCCAAACCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                             <  1:27454/150‑1 (MQ=255)                                                                                                                    gaCATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                         >  1:338870/1‑151 (MQ=255)                                                                                                                      cATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                       >  2:430997/1‑151 (MQ=255)                                                                                                                       aTCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                       >  1:1232462/1‑150 (MQ=255)                                                                                                                       aTCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                       >  1:1337718/1‑150 (MQ=255)                                                                                                                           aCGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                     <  1:373907/148‑1 (MQ=255)                                                                                                                              tttCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa                <  2:1232190/150‑1 (MQ=255)                                                                                                                               ttCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc              >  1:141970/1‑151 (MQ=255)                                                                                                                                tCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc             >  2:580196/1‑151 (MQ=255)                                                                                                                                   ccTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          >  2:1064674/1‑151 (MQ=255)                                                                                                                                   ccTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          >  1:421657/1‑151 (MQ=255)                                                                                                                                    cTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt         >  2:621977/1‑151 (MQ=255)                                                                                                                                         cAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga    <  2:103064/151‑1 (MQ=255)                                                                                                                                          aaaTCCTCCGTTTCACCCCCCGCACGGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa   <  2:608770/151‑1 (MQ=255)                                                                                                                                          aaaTCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga    >  1:387447/1‑150 (MQ=255)                                                                                                                                           aaTCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc  >  1:1299819/1‑151 (MQ=255)                                                                                                                                           aaTCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc  >  1:598858/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                  GAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  NC_000913/3966533‑3966821

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