Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)7→8	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	98.7%	277.1 / ‑4.6	77	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (41/35);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGG  >  NC_000913/1597158‑1597449                                                                                                                                                        |                                                                                                                                                 aGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTaacca‑a‑‑cc                                                                                                                                               >  1:189112/1‑151 (MQ=255) aGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTaacca‑a‑‑cc                                                                                                                                               <  1:702297/151‑1 (MQ=255) aGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTaacca‑a‑‑cc                                                                                                                                               <  1:1236293/151‑1 (MQ=255)  gTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑ccc                                                                                                                                              >  2:1347904/1‑151 (MQ=255)    ccATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑ccccc                                                                                                                                            >  1:1180026/1‑151 (MQ=255)     cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑cccccc                                                                                                                                           >  2:1454680/1‑151 (MQ=255)     cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑cccccc                                                                                                                                           >  1:1099761/1‑151 (MQ=255)      attattGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑ccccccc                                                                                                                                          <  1:376202/151‑1 (MQ=255)         attGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑cccccccctc                                                                                                                                       <  1:140027/151‑4 (MQ=255)             gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTcgcg                                                                                                                                     >  2:1460955/1‑150 (MQ=255)              aaaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGa                                                                                                                                    >  2:138821/1‑150 (MQ=255)              aaaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGa                                                                                                                                    >  1:947262/1‑150 (MQ=255)               aaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAAt                                                                                                                                  >  2:194095/1‑151 (MQ=255)                  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAg                                                                                                                                <  2:189014/150‑1 (MQ=255)                  atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGt                                                                                                                               >  2:581096/1‑151 (MQ=255)                      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGcc                                                                                                                           >  1:216101/1‑151 (MQ=255)                      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGcc                                                                                                                           >  1:888502/1‑151 (MQ=255)                      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGcc                                                                                                                           <  1:499202/151‑1 (MQ=255)                            ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCa                                                                                                                     <  2:999119/151‑1 (MQ=255)                             cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAg                                                                                                                    <  2:1118234/151‑1 (MQ=255)                              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                                                   >  2:408943/1‑151 (MQ=255)                                 aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                                                                                >  2:877920/1‑151 (MQ=255)                                 aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                                                                                >  1:612455/1‑151 (MQ=255)                                       atatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGCCGCCTCCCCAGATCCTTGAtt                                                                                                          >  1:1357419/1‑151 (MQ=255)                                            cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                                                     <  1:1382908/151‑1 (MQ=255)                                            cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                                                     <  2:22278/151‑1 (MQ=255)                                               tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                                                  >  2:1286758/1‑151 (MQ=255)                                               tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                                                  >  1:721920/1‑151 (MQ=255)                                                ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                                                  <  1:30653/150‑1 (MQ=255)                                                 taGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                                                                 >  2:181209/1‑150 (MQ=255)                                                 taGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                                                                 >  1:552527/1‑150 (MQ=255)                                                  agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                                                                                <  1:919520/150‑1 (MQ=255)                                                    agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                                                             >  2:1185833/1‑151 (MQ=255)                                                      tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                                                                           >  2:422512/1‑151 (MQ=255)                                                         gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                                                                        <  1:697853/151‑1 (MQ=255)                                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                                       >  2:627498/1‑151 (MQ=255)                                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                                       <  2:25876/151‑1 (MQ=255)                                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                                                                   >  1:492815/1‑151 (MQ=255)                                                              ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                                                                   <  2:1169958/151‑1 (MQ=255)                                                                  gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAt                                                                                <  1:1068012/150‑1 (MQ=255)                                                                     gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                            <  1:88052/151‑1 (MQ=255)                                                                     gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                            >  2:385084/1‑151 (MQ=255)                                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                                                          <  1:100903/151‑1 (MQ=255)                                                                       cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                                                          >  2:910520/1‑151 (MQ=255)                                                                        agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                                                                         >  1:738058/1‑151 (MQ=255)                                                                          agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                                       >  2:673969/1‑151 (MQ=255)                                                                            tCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                     >  2:1016656/1‑151 (MQ=255)                                                                             cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                                       >  1:1280628/1‑148 (MQ=255)                                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                     <  1:1092205/148‑1 (MQ=255)                                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                     <  2:522376/148‑1 (MQ=255)                                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                     <  1:148670/148‑1 (MQ=255)                                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                                                    <  2:1276218/149‑1 (MQ=255)                                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                                                                <  1:483165/151‑1 (MQ=255)                                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                                                                <  1:617391/151‑1 (MQ=255)                                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCACGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                                                                  >  1:1137201/1‑149 (MQ=255)                                                                                             cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc                                                     <  1:1454787/150‑1 (MQ=255)                                                                                             cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCCAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGTGGAATCATTATTTTTTATTATCCACTCGCCAcc                                                    >  1:956790/1‑151 (MQ=255)                                                                                                 gTTCCGATTCCCGTTGCATTGGCCCCGGTAATCGTATTTTCAGCCGGTACCCACC‑‑CCCCCCCTCGCGAATGGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCtt                                                 <  2:978265/151‑1 (MQ=255)                                                                                                  ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCtttt                                               >  2:729813/1‑151 (MQ=255)                                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTc                                              <  2:1217602/150‑1 (MQ=255)                                                                                                     cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCaa                                            <  2:104381/151‑1 (MQ=255)                                                                                                       aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                                          <  1:271053/151‑1 (MQ=255)                                                                                                         aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGAc                                        <  1:909551/151‑1 (MQ=255)                                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACaa                                      <  1:378585/151‑1 (MQ=255)                                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCACGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACa                                       >  2:823140/1‑150 (MQ=255)                                                                                                               ggCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTg                                  <  1:761619/151‑1 (MQ=255)                                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                 >  2:851445/1‑151 (MQ=255)                                                                                                                 cATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGc                                <  1:1348011/151‑1 (MQ=255)                                                                                                                    tGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCa                              <  1:1197621/150‑1 (MQ=255)                                                                                                                       cACCGATAATCGTATTTTCAGCC‑GTAACCA‑ACCCCCCCCCCCGCGAATATTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTTTTTTTTATTAGCCACTCGCCACCCTTTTCACCGACAACCTGGCCAt                             >  1:989970/1‑149 (MQ=255)                                                                                                                               aTCGTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTTCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGcc                  >  1:92852/1‑151 (MQ=255)                                                                                                                                  gTATTTTCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGa               <  1:194117/151‑1 (MQ=255)                                                                                                                                     ttttCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGAt              >  1:734584/1‑149 (MQ=255)                                                                                                                                     ttttCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGAt              >  2:1467729/1‑149 (MQ=255)                                                                                                                                      tttCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAg            <  1:729892/150‑1 (MQ=255)                                                                                                                                      tttCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGc           >  2:1401797/1‑151 (MQ=255)                                                                                                                                        tCAGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCt         >  1:497201/1‑151 (MQ=255)                                                                                                                                          aGCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTa        >  1:463944/1‑150 (MQ=255)                                                                                                                                           gCC‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTaaa      <  2:761967/151‑1 (MQ=255)                                                                                                                                             c‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCGCCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAAt    <  2:1139620/151‑1 (MQ=255)                                                                                                                                             c‑GTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAAt    >  2:113287/1‑151 (MQ=255)                                                                                                                                             c‑GTAACCA‑AC‑CCCCCCCTCGCGAATAATCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTaaaa     >  1:1149938/1‑150 (MQ=255)                                                                                                                                               gTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAAt    >  1:1093850/1‑150 (MQ=255)                                                                                                                                               gTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg   <  2:287076/151‑1 (MQ=255)                                                                                                                                               gTAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg   <  2:888414/151‑1 (MQ=255)                                                                                                                                                tAACCA‑AC‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATgg  >  1:276860/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                 AGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCC‑GTAACCA‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGG  >  NC_000913/1597158‑1597449

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